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Abstract. We continue the development, by reduction to a first order system 
for the conormal gradient, of a priori estimates and solvabihty for boundary 
value problems of Dirichlet, regularity, Neumann type for divergence form second 
order, complex, elliptic systems. We work here on the unit ball and more generally 
its bi-Lipschitz images, assuming a Carleson condition as introduced by Dahlberg 
which measures the discrepancy of the coefficients to their boundary trace near 
the boundary. We sharpen our estimates by proving a general result concerning 
a priori almost everywhere non-tangential convergence at the boundary. Also, 
compactness of the boundary yields more solvability results using Fredholm theory. 
Comparison between classes of solutions and uniqueness issues are discussed. As 
a consequence, we are able to solve a long standing regularity problem for real 
equations, which may not be true on the upper half-space, justifying a posteriori 
a separate work on bounded domains. 
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1. Introduction and main results 

We refer to where we took up this study, for a comprehensive historical account 
of the theory of boundary value problems for second order equations of divergence 
form. Before we come to our work here, let us connect more deeply to even earlier 
references going back to the seminal work of Stein and Weiss [39] that paved the 
way for the development of Hardy spaces on the Euclidean space in several 
dimensions. Their key discovery was to look at the system of differential equations in 
the upper-half space satisfied by the gradient F = {dfU, Vxu) of a harmonic function 
u on the upper half-space, to which they gave the name of conjugate system or M. 
Riesz system. The system of differential equations is in fact a generalized Cauchy- 
Riemann system which can be put into a vector-valued ODE form. They did not 
exploit this ODE structure but used instead subharmonicity properties of \F\p for 
p > to define the (harmonic) Hardy spaces as the space of those conjugate 
systems satisfying 

sup / \F{t, x)f dx < oo 

and to prove that the elements in this space have boundary values 

F{t,x) -> F{0,x) 

in the norm and almost everywhere non-tangentially. Further, they proved that 
elements in can be obtained as Poisson integrals of their boundary traces. In 
other words, there is a one-to-one correspondence between and its trace space 
"H^. By using Riesz transforms, the trace space "H^ is in one-to-one correspondence 
with the space defined by taking the first component of trace elements. As they 
pointed out, it was nothing new for j» > 1 as we get L^, but for p < 1 it gave a new 
space. Over the years, this last space turned out to have many characterizations, 
including the ones with Littlewood-Paley functionals of Fefferman and Stein |25] 
and the atomic ones of Coifman [18] and Latter [32] , and is now part of a rich and 
well understood family of spaces. 

In our earlier work with Mcintosh [S], and in [1], we wrote down the Cauchy- 
Riemann equations corresponding to the second order equation and the key point 
was a further algebraic transformation that transformed this system to a vector- 
valued ODE. In some sense, we were going back in time since elliptic equations with 
non-smooth coefficients have been developed by other methods since then (see [28]). 
In this respect, it is no surprise in view of the above discussion that we denote our 
trace spaces by "H. They are in a sense generalized Hardy spaces, and this notation 
was used as well in our earlier work with Hofmann [5]. We shall use again such 
notation and terminology here. What today allows the methods of Hardy spaces 
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to be applicable in the case of non-smooth coefficients, are the quadratic estimates 
related to the solution of the Kato conjecture for square roots. These are a starting 
point of the analysis. Indeed, the quadratic estimates are equivalent to the fact 
that two Hardy spaces split the function space topologically, as it is the case for the 
classical upper and lower Hardy spaces in complex analysis, essentially from the F. 
and M. Riesz theorem on the boundedness of the Hilbert transform. So in a sense 
everything looks like the case of harmonic functions (for p = 2 at this time). But this 
is not the case. The difference is in the last step, taking only one component of the 
trace of a conjugate system. This may or may not be a one-to-one correspondence, 
which translates to well- or ill-posedness for the boundary value problems of the 
original second order equation. 

See also [11] for a different generalization of Stein-Weiss conjugate systems of 
harmonic functions. There conjugate differential forms on Lipschitz domains were 
constructed by inverting a generalized double layer potential equation on the bound- 
ary. 

Let us introduce some notation in order to state our results. Our system of 
equations is of the form 



where di = < i < n and the matrix of coefficients is A = (^"^ (x))°^^q G 

Loo{^'t C{C^^~^'^^"^)), n,m > 1. We emphasize that the methods used here work 
equally well for systems (m > 2) as for equations (m = 1). For the time being, 
n = 0^+" := {x G R^+" ; |x| < 1} for the unit ball in Ri+" (see the end of 
the introduction for more general Lipschitz domains). The coefficient matrix A is 
assumed to satisfy the strict accretivity condition 



for some k > 0, uniformly for a.e. r G (0, 1) and u G C^(0^+"; C") where we use 
polar coordinates x = rx, r > 0,a; G 5", and dx is the standard (non-normalized) 
surface measure on S*" = dO^^"". The optimal k is denoted ka. This ellipticity 
condition is natural when viewing A as a perturbation of its boundary trace. See 
below. 

The boundary value problems we consider are to find u G P'(0^^"; C™) solving 
in distribution sense, with appropriate interior estimates of VxM and Dirichlet 
data in L2, or Neumann data in L2, or regular Dirichlet data with gradient in L2. 
Note that since we shall impose distributional VxM G L^^'^, u can be identified with 
a function u G W^' °'^(0^+", C"), i.e. with a weak solution. In order to study these 
boundary value problems, our task, and this is the ffist main core of the work, is to 
obtain a priori estimates. 

As in the previous work [1] on the upper-half space R+'^'^, we reduce ([I]) to a ffist 
order system with the conormal gradient as unknown function, so the strategy and 
the scale-invariant estimates are similar. See the Road Map Section in [1] for an 
overview. Some changes will arise in the algebraic setup and in the analysis though. 
Here, the curvature of the boundary (the sphere) will play a role in the algebraic 
setup, making the unit circle slightly different from the higher dimensional spheres. 



(1) 
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In addition, owing to the fact that the boundary is compact, we may use Fredholm 
theory to obtain representations and solvability by only making assumptions on 
the coefficients near the boundary. We shall focus on this part here and give full 
details. We also mention that the whole story relies on a quadratic estimate for 
a first order bisectorial operator acting on the boundary function space. On the 
upper-half space, this estimate was already available from [H] as a consequence of 
the strategy to prove the Kato conjecture on R". We shall need to prove it on 
the sphere, essentially by localization and reduction to [13], where such estimates 
were proved for first order operators with boundary conditions. An implication of 
independent interest is the solution to the Kato square root on Lipschitz manifolds. 
This is explained in Section [HI 

As is known already for real equations (m = 1) from work of Caffarelli, Fabes and 
Kenig [16] , solvability requires a Dini square regularity condition on the coefficients 
in the transverse direction to the boundary. So it is natural to work under a condition 
of this type. We use the discrepancy function and the Carleson condition introduced 
by Dahlberg [21]. For a measurable function / on 0^+", set 

(3) r(x):=esssup|/(y)|, 

ygiy°{x) 

where W"{x.) denotes a Whitney region around x G O^"^" and 

(4) 11/11^:= sup // /*(x)2^^y for some fixed c < 1, 

where the supremum is over all geodesic balls Q C S"" of radius r{Q) < c. We 
make the standing assumption on A throughout that there exists Ai a measurable 
coefficient matrix on S*", identified with radially independent coefficients in O^"*""", 
such that S{y) := A{y) — Ai{y),y = y/|y|, satisfies the large Carleson condition 

(5) ll^'llc < oo. 

The choice of c is irrelevant. Note that this means in particular that S* vanishes 
on 5*" in a certain sense and so Ai{y) = yl(y/|y|). In fact, it can be shown as in 
[H Lemma 2.2] that if there is one such Ai, it is uniquely defined, ||v4i||oo < ll^lloo 
and > i^A- So we call Ai the boundary trace of A. It turns out that this is a 
very natural assumption with our method, implying a wealth of a priori information 
about weak solutions as stated in Theorem 11.11 Such a result applies in particular 
to all systems with radially independent coefficients since £^ = in that case. 

For a function / defined in O^"*"", its truncated modified non-tangential maximal 
function is defined as in [SU] by 



(6) iV:(/)(x):= sup (iWirx)]-' [ \f{y)\'dy)\ x e S\ 

1-T<r<l \ Jw°(rx) / 

for some fixed r < 1. Note that changing the value of r will not affect the results. 
We shall use the notation fr{x) := f{rx) for < r < l,x G S^. Our main result is 
the following. 

Theorem 1.1 (A priori representations and estimates, existence of a trace, Fa- 
tou type convergence). Consider coefficients A G Loo(0^"'"'^; £(0^^"'"'^^'")) which are 
strictly accretive in the sense of ([2]) and satisfy ([5]) with boundary trace Ai. Consider 
u G W^''^"'^ {O^^'^ ] C"*) which satisfies (QP in 0-*^+" distributional sense. 
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(i) // ||Ar°(VxM)||L2(5") < oo, then 
(a) Vx^i has limit 
1 



lim 

r->l 1 — r 



|Vxm(x) -^i(x)prfx = 



r<|x|<(l+r)/2 



for some g, G 1.2(5- C(i+")-) wtth ||^?l|L,(s^c(l+")"') < I|A^°(Vxm)||l,(S"). 

(b) r Ur belongs to C(0, 1; L2(S""; C")) and has L2 limit ui at the 
boundary with 

\\Ur - Ui\\l2{S";C"') ^ 1 - ^ 

and Ml e 1^2^(5"; C™). 

(c) Fatou type results: For almost every x G , 



\im\W°{rx)\-^ 

r— >■! 



lim|iy°(rx)|-^ 



lim|W^''(rx)r^ 

r— S>1 



W°(rx) 

dtu{y)dy = (^i)x(x), 
(AVxM)ii(y)f/y = 



Ty°(r2:) 



and if m = 1 (equations) or n = 1 (unit disk) we also have 



lim \ W°{rx) 
1 — ^1 



lim\W°{rx)\-^ 
1 — ^1 



Vx?i(y)dy = gi{x), 

W°{rx) 

{AV^u){y)dy = {Aigi){x). 



'W°(rx) 

ill) If /qi+„ |VxUp(l - |x|)(ix < cx), i/ien 

(a) r ^ Ur belongs to C(0, 1; L2{S"'] C")) and has L2 limit 

lim \\Ur — ttl||L2(S" C"M = 



for some Ui E L2{S^; C"). 

(b) We have a priori estimates 



(8) 



|A^»IIL(5")< 



|Vxm|^(1 - |x|)dx + 



Ql + n 



Ui{x)dx 



(9) ||^r||i2(5n;C'") ~ 



|Vxu|^(l - |x|)dx 



Ql + n 



Ui{x)dx 



5" 



(c) Fatou type results: For almost every x G S^, 

u{y)dy = ui{x). 



r G (0,1). 



\Ym\W°{rx)\-^ 

'■"^^ JW°{rx) 

The definition of the normal component (■)^ and tangential part (■)|| of a vector 
field will be given later. Not stated here are representation formulas giving ansatzes 
to find solutions as they use a formalism defined later. In particular, we introduce a 
notion of a pair of conjugate systems associated to a solution. We note that the non- 
tangential maximal estimate ([H]) was already proved in the R^"^" setting of [1]. Again, 
this is an a priori estimate showing that, under the assumption ||£^||c7 < 00, the 
class of weak solutions with square function estimate /qi+„ |Vx^ip(l — |x|)(ix < 00 
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is contained in the class of weak solutions with non-tangential maximal estimate 
||-^°(w)||2 < oo. The almost everywhere convergences of Whitney averages are new. 
They apply as well to the setup in [1]. 

Theorem [LT] enables us to make the following rigorous definition of well-posedness 
of the EVPs. 

Definition 1.2. Consider coefficients A E L^{0^~^^; £{0^^^"^^^)) which are strictly 
accretive in the sense of (E]). 

• By the Neumann problem with coefficients A being well-posed, we mean 
that given ip G L2(S'";C'") with Jg„ip{x)dx = 0, there is a function u G 
p^^i,ioc^Qi+„. with estimates || A^°(Vx^i) ||l2(5") < oo, unique modulo con- 
stants, solving ([1]) and having trace gi = limf,_s.i(Vx^i)r in the sense of ([7]) 
such that (AiQi)^ = ip. 

• Well-posedness of the regularity problem is defined in the same way, but 
replacing the boundary condition (AiQi)^ = (p by {gi)\\ = (p, for a given 

• By the Dirichlet problem with coefficients A being well-posed, we mean that 
given if e L2(^";C™), there is a unique function u e •'°'(Oi+"; C"^) 
with estimates /qi+„ |VxMp(l — |x|)(ix < oo, solving ([T]) and having trace 
limr_j.i Ur = (p iTi the sense of almost everywhere convergence of Whitney 
averages. 

For the Neumann and regularity problem when \\S\\c < oo, for equations (m = 1) 
or in the unit disk (n = 1) or any system for which A is strictly accretive in pointwise 
sense, the trace can also be defined in the sense of almost everywhere convergence of 
Whitney averages of VxM and the same for the conormal derivative (AVxm)x- The 
operator V5 denotes the tangential gradient. See Section [HI 

For the Dirichlet problem, the trace is defined for the almost everywhere conver- 
gence of Whitney averages. When ||£^||c < 00, Theorem II . II shows that it is the same 
as the trace in L2 sense. We remark that we modified the meaning of the boundary 
trace in the definition of the Dirichlet problem compared to jl]. This modification 
can be made there as well and the same results hold. 

We now come to our general results on these BVPs. A small Carleson condition, 
but only near the boundary, is further imposed to obtain invertibility of some oper- 
ators. The second result is on the precise relation between Dirichlet and regularity 
problems. The first and third are perturbations results for radially dependent and 
independent perturbations respectively. The last is a well-posedness result for three 
classes of radially independent coefficients. 

Theorem 1.3. Consider coefficients A G Loo(0^"^"'; £(C^^"''"''''")) which are strictly 
accretive in the sense of i^. Then there exists e > 0, such that if A satisfies the 
small Carleson condition 

(10) lim ||xr<r<i(^-^i)||c < e 

and the Neumann problem with coefficients Ai is well-posed, then the Neumann 
problem is well-posed with coefficients A. 

The corresponding perturbation result for the regularity and Dirichlet problems 
also holds. For the Neumann and regularity problems, the solution u for datum (p 
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has estimates 

I |Vxn|2rfx< ||iV:(Vxn)||2^||^||2. 

^|x|<l/2 

For the Dirichlet problem, the solution u for datum ip has estimates 
\\K{u)\\l^ sup / |Vxwr(l-|x|)dx + 

l/2<r<l Joi+" 

An ingredient of the proof is the following relation between Dirichlet and regularity 
problems, in the spirit of [3U1 Thm. 5.4]. 

Theorem 1.4. Consider coefficients A G Loo(0^"'""'; £(C^^"^"''''")) which are strictly 
accretive in the sense of Then there exists e > 0, such that if A satisfies the 
small Carleson condition ffTUj) . then the regularity problem with coefficients A is 
well-posed if and only if the Dirichlet problem with coefficients A* is well-posed. 

Theorem 1.5. Consider radially independent coefficients Ai G Loo^S"", C{C^^~^"'^"^)) 
which are strictly accretive in the sense o/ ([21). // the Neumann problem with co- 
efficients Ai is well-posed, then there exists e > such that the Neumann problem 
with coefficients A[ G Loo(S^] C(C^^~^^^'^)) is well-posed whenever \\Ai — A[\\oo < e. 
The corresponding perturbation results for the regularity and Dirichlet problems also 
hold. 

Theorem 1.6. Consider radially independent coefficients Ai G Loo(>S'"; £(C'-^'^"''"*)) 
which are strictly accretive in the sense of ([2]). The Neumann, regularity and Dirich- 
let problems with coefficients Ai are well-posed if 

(1) either Ai is Hermitean, i.e. 74^ = ^41, 

(2) or block form, i.e. {Ai)^^^ = = (Ai)||x in the normal/tangential splitting of 
C(i+n,)m ^^gg Section\^, 

(3) or Ai has Hdlder regularity CiS""; £(C(i+'^)'^)), s > 1/2. 

Proof of Theorems \1.3[ \1.5\ and \1.6i . For Theorem 11.11 the L2-limits and L2- 
estimates of solutions follow from Theorem 112.41 and Corollary 112.81 respectively. 
The non-tangential maximal estimate ([H]) is in Theorem 114.11 Almost everywhere 
convergence of averages follows from Theorems 115.21 and 115.51 

The well-posedness results in Theorem 11.61 are in Propositions 117.16^ 117.151 and 
117.171 The radially independent perturbation result in Theorem 11.51 is in Corol- 
lary 117.131 The well-posedness result for radially dependent coefficients with good 
boundary trace in Theorem 11.31 is in Proposition 117.141 □ 

Our next result is the following semigroup representation, analogous to the result 
in [3] in the upper half-space. It is interesting to note that for harmonic functions 
u, it gives a direct proof (without passing through non-tangential maximal function 
or sup —L2 estimates) that Jqi+u |Vx^ip(l — |x|)(ix < 00 implies a representation by 
Poisson kernel from its trace (also shown to exist). We have not seen this argument 
in the literature. Another interesting feature is that it points out the importance of 
well-posedness of the Dirichlet problem when dealing with more general coefficients. 

Theorem 1.7. Consider radially independent coefficients Ai G Loo(>S'"; £(0^^"*^"^™)) 
which are strictly accretive in the sense of ([2]). Assume that the Dirichlet problem 
with coefficients Ai is well-posed. Then the mapping 

Vr : L2{S''; C™) -> L2(^"; C"^) : Ui ^ Ur, 



ip{x)dx 



5" 
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where u is the solution to the Dirichlet problem with datum Ui, defines a bounded 
operator for each r G (0, 1]. The family {Vr)re{o,i] is a multiplicative Co-semigroup 
( i. e. VrVr' = Vrr' o-^d Vr — ?■ / strongly in L2 when r ^ 1) whose infinitesimal 
generator A (i.e. Vr = e^^'^^'^-^) has domain D{A) contained in 14^2^(5'"; C"). More- 
over, D{A) = 14^2^(5'"; C™") if and only if the Dirichlet problem with coefficients 
is well-posed. 

As mentioned above two classes of weak solutions compare: the one with square 
function estimates is contained in the one with non-tangential maximal control. It 
is thus interesting to examine this further. Does the opposite containment holds? 
How do well-posedness in the two classes compare? Clearly uniqueness in the larger 
class implies uniqueness in the smaller, and conversely for existence. As we shall 
see, positive answers come a posteriori to solvability. 

Definition 1.8. The Dirichlet problem with coefficients A is said to be well-posed in 
the sense of Dahlberg if, given ip G L2(S'"; C™), there is a unique weak solution u G 
py^i,iocj^Qi+„,. diVx^VxM = with estimates ||iV°(M)||2 < c>o and convergence 

of Whitney averages to (f, almost everywhere with respect to surface measure on 

This definition has the merit to be natural not only for equations but for systems as 
well. For real equations, this is equivalent to the usual one as N° can be replaced by 
the usual pointwise non-tangential maximal operator by the DeGiorgi-Nash-Moser 
estimates on weak solutions. Even in this case, observe that the control ||A'^^('u)||2 < 
oc does not enforce the almost everywhere convergence property. Thus existence of 
the limit is part of the hypothesis in Definition II. 8[ as compared to Definition II. 2[ 
A first result is the following. 

Theorem 1.9. Consider radially independent coefficients Ai G Loo(5'"; C{C^^~^"^'^)) 
which are strictly accretive in the sense of ([2]). Assume that the Dirichlet and 
regularity problems with coefficients Ai are well-posed in the sense of Definition 
Then, all weak solutions to diVxAiVxW = with ||A^°(m)||2 < 00 are given by the 
semigroup of Theorem \1.7\ In particular, the Dirichlet problem with coefficients Ai 
is well-posed in the sense of Dahlberg. 

Theorem 11.41 implies the same conclusion for the coefficients AI. The next results 
are only for real equations where the theory based on elliptic measure brings more 
information. For (complex) equations, the strict accretivity in the sense of is 
equivalent to the usual pointwise accretivity, which is the same as the strict ellipticity 
for real coefficients. 

Theorem 1.10. Consider an equation with real coefficients A G Loo{0^'^"'; C(lV''^'^)), 
which are strictly elliptic. Assume further that the small Carleson condition ( fTOl) 
holds. Then the following are equivalent. 

(i) The Dirichlet problems with coefficients A and A* are well-posed in the sense 
of Dahlberg. 

(ii) The Dirichlet problems with coefficients A and A* are well-posed in the sense 
of Definition \1.2[ 

Moreover, in this case the solutions for coefficients A (resp. A*) from a same datum 
are the same. 
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Note that, by Theorem 11.41 we can replace (ii) by (ii'): the regularity problems 
with coefficients A and A* are well-posed. When A = A*, all the problems in 
(i) and (ii') are well-posed by [30] so there is nothing to prove. For (even non- 
symmetric) real coefficients A alone, the direction from (ii') to (i) was known from 
[30] (without assuming the Carleson condition) and the converse is unknown. It 
seems that making the statement invariant under taking adjoints solves the issue. 
We mention the equivalence in [31] concerning Lp versions of this statement for self- 
adjoint constant coefficient systems on Lipschitz domains (in this case, the L2 result 
is known and used). 

Our last result is well-posedness of the regularity problem under a transversal 
square Dini condition on the coefficients, analogous to the result obtained in [21] 
for the Dirichlet problem with real and symmetric A. This partly answers Problem 
3.3.13 in [28]. 

Theorem 1.11. Consider an equation with coefficients A € Loo(0"'^'''"; £(C^"'"'^)), 
which are strictly accretive in the pointwise sense. Then there exists e > 0, such that 
if A satisfies the small Carleson condition ffTUl) and its boundary trace Ai is real 
and continuous, then the Dirichlet problem with coefficients A is well-posed in the 
sense of Definition \l.B and in the sense of Dahlberg, and the regularity problem with 
coefficients A is well-posed. In particular, this holds if A is real, continuous in 0^+" 
and the Dini square condition ^Qw\{t)^ < 00 holds, where WA{t) = sup{|y4(ra:) — 
A{x)\]x G S*", 1 — r < t}. The corresponding results hold in for the Neumann 
problem with coefficients A. 

Proofs of Theorems 11.71 and 11.91 are in Sections [TH] and proofs of Theorems 11.101 
and 11.111 are in Section (THl 

We end this introduction with a remark on the Lipschitz invariance of the above 
results. Let Vt C R^"*"" be a domain which is Lipschitz diffeomorphic to O^"*"" and let 
p : O^"*"" — be the Lipschitz diffeomorphism. Denote the boundary by S := dVt 
and the restricted boundary Lipschitz diffeomorphism by po : "S*" — S. 

Given a function -u : f2 — )■ C™, we pull it back io u := u o p : O^"^" — ?• C™. 
By the chain rule, we have VxM = P*(Vy'u), where the puUback of an m-tuple 
of vector fields /, is defined as p*(/)(x)" := p*(x)/°'(p(x)), with p* denoting the 

transpose of Jacobian matrix p. If u satisfies diVyAVyCt = in 1], with coefficients 

A e Loo(l^;/:(C(^+")"*)), then u will satisfy div^AV^u = in 0^+", where A e 
Loo(Oi+";£(C(i+")™)) are the "pulled back" coefficients defined as 

(11) A(x) := I J(p)(x)|(p(x))-ii(p(x))(p*(x))-\ X e Ol+^ 

Here J(p) is the Jacobian determinant of p. 

The Carleson condition, non-tangential maximal functions and square functions on 
Q correspond to ones on S"" under pullback, so that 1 — |x| becomes S{y) the distance 
to S. In particular, the condition for A amounts to \\S\\c < 00 with S defined from 
A. We remark that puUbacks allow to replace normal directions by oblique (but 
transverse) ones to the sphere in the Carleson condition on the coefficients: take 
p : O^"*""" — > O^"*"" to be a suitable Lipschitz diffeomorphism. 

The boundary conditions on -u on S translate in the following way to boundary 
conditions on m on S*". 

• The Dirichlet condition -u = on S is equivalent to m = (/9 on S*", where 
(/^:=^opoeL2(^";C™). 
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• The Dirichlet regularity condition V^u = on S (Vs denoting the tangential 
gradient on E) is equivalent to V^u = ip on 5"", where (p := Pol*^) ^ C 
L2(5";C"™). 

• The Neumann condition {u, AVyu) = (p on H {u being the outward unit 
normal vector field on S) with ip{y)dy = is equivalent to (n, AVxm) = ^ 
on 5" with (p{x)dx = 0, where v? := | J(po)|<^ o Po e L2(^"; C"^). 



In this way the Dirichlet /regularity/Neumann problem with coefficients A in the 
Lipschitz domain Q is equivalent to the Dirichlet /regularity/Neumann problem with 
coefficients A in the unit ball O^^", and it is straightforward to extend the results 
on O^''"" above to Lipschitz domains Q. 

The plan of the paper is as follows. In Section |2l we transform the second order 
equation ([T]) into a system of Cauchy-Riemann type equations. In Section [3], the 
Cauchy-Riemann equations are integrated to a vector-valued ODE for the conormal 
gradient of u and a second ODE is introduced to construct a vector potential. The 
infinitesimal generators Dq and Dq for these ODE with radially independent coef- 
ficients are studied in Sections H] and [6l and it is shown in Section [7] that Dq and 
Dq have bounded holomorphic functional calculi. Section [5] treats special features of 
elliptic systems in the unit disk. In Section [9] we define the natural function spaces 
X° and y° for the BVPs and we describe in Section [TU] how to construct solutions 

from the semigroups generated by |-Do| = ^/D^ and |-Do| = In Section [TTl the 

ODE with radially dependent coefficients for the conormal gradient from Section H] 
is reformulated as an integral equation involving an operator 5*^, which is shown to 
be bounded on the natural function spaces for the BVPs. In Section [12], we obtain 
representation for X°- and 3^°-solutions. These representations are further devel- 
oped in Section [T3] where we introduce the notion of a pair of conjugate systems for 
([T]), allowing to prove in Sections \\M and [T51 non-tangential maximal estimates and 
Fatou type results. Crucial for the solvability of ([T]) is the invertibility oi I — Sa- In 
Section [121 we apply Fredholm theory to show that / — S'a is invertible on the nat- 
ural spaces whenever the small Carleson condition fllOl) holds, which proves that it 
suffices to assume transversal regularity of A near the boundary only. (For BVPs on 
the unbounded half-space studied in [1], the needed compactness was not available.) 
We then study well-posedness in Section [T71 this is where we prove the perturbation 
results, the equivalence Dirichlet /regularity up to taking adjoints and obtain classes 
of radially independent coefficients for which well-posedness holds. The section [18] 
deals with uniqueness issues, on comparisons of different classes of solutions upon 
some well-posedness assumptions. We conclude the article in Section [TH] by a dis- 
cussion in the special case of real equations (m = 1) for which we obtain further 
results. 

Acknowledgements. The support of Mathematical Science Institute of the 
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greatly acknowledged. We thank Steve Hofmann for several discussions on real 
equations pertaining to this work. We also thank Alano Ancona, Carlos Kenig and 
Martin Dindos for a few exchanges. 
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2. Generalized Cauchy-Riemann system 

Following [SI El H], the starting point of our analysis is that solving for u the 
divergence form system ([T]) amounts to solving for its gradient g a system of Cauchy- 
Riemann equations. 

Proposition 2.1. Consider coefficients A G L^{0^+''; £{0'-^+''^"')). If u is a weak 
solution to diVxAVxM = m 0^+", then g := W^u G c{i+n)m-^ ^ 

solution of the generalized Cauchy-Riemann system 

'div, (Ag) = 0, 
curlxS- = 0, 



(12) 



in O"'^^"\{0} distribution sense. Conversely, if g G L]^^{0^~^^; ig q. solution 

to (in]) m 0^+"" \ {0} distribution sense, then there exists a weak solution u to 
diVxAVxW = in O^"*"", such that g = Vx^ in 0^~^" distribution sense. 

Proof. If u is given, then g := Vx^i has the desired properties and the equation 
is even satisfied in O^"^*^ distribution sense. Conversely, assume g is given and 
satisfies ( fT2i) in O^^" \ {0} distribution sense. Then the next lemma applied to both 
operators divx and curlx implies that is a removable singularity and that ( lT2l) holds 
in O^"'""' distribution sense. Thus one can define a distribution u in O^^" such that 
g = VxM, hence diVxAVxW = in O^'^". That m is a weak solution follows from 
^ G 4°=(Oi+"; C(^+")"*). □ 

Lemma 2.2. Let X be a homogeneous first order partial differential operator on 
j^i+n jYiappifig Q^-valued distributions to C^-valued distributions, k,i & Z_|_. // 
h G L^°=(0^+"; C^') and Xh = in distributional sense on 0^+^^ \ {0}, then Xh = 
in O^^"- -distributional sense. 

Proof Let G C^{0^+'';C^). We need to show that /i)rfx = 0. To 

this end, let ?7e be a smooth radial function with 77^ = on {|x| < e}, 77^ = 1 on 
{2e < |x| < 1} and \\Vr],\\^ < e'^ Then 

[ r],{X*(j),h)d^= [ (X*(r/,0),/i)rfx- [ {{X*r],)(P, h)d^ 

{{X%)<f), h)d^. 

As e — 0, the left hand side converges to Jqi+„(X*</), /i)(ix, whereas 

/ ((X*7^,)0, h)d^ <- [ \h\d^< e(«-i)/2 ( [ \h\'d^) ^'^ ^ 0. 

^ Je<|x|<2e \Je<|x|<2e / 

This proves the lemma. □ 
3. The divergence form equation as an ODE 

We introduce a convenient framework to transform the Cauchy-Riemann system 
into an ODE. 

We systematically use boldface letters x, y, . . . to denote variables in R^^" and 
indicate the variable for differential operators in R^"*"" (e.g. Vx...). We denote 
points on S*" by x, . . . and the standard (non-normalized) surface measure on S"^ 
by dx. Polar coordinates are written x = rx, with r > and x G S*". For a function 
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/ defined in O^"*"", we write fr{x) := f{rx), x G S"^, for tlie restriction to tlie spliere 
witli radius < r < 1, parametrized by S"'. 

The radial unit vector field we denote by n = n(x) := x/|x|. Vectors v G R"^+", 
we split V = + where fx := {v,n) is the normal component and f|| := v — v^n 
is the angular or tangential part of v, which is a vector orthogonal to n. Note that Vj_ 
is a scalar, but is a vector. In the plane, i.e. when n = 1, we denote the counter 
clockwise angular unit vector field by r, and we have v = v^n + (f , r)r. For an 
m-tuple of vectors v = {v°')i<a<m, "we define its normal components and tangential 
parts componentwise as 



The tangential gradient, divergence and curl on the unit sphere are denoted by 
V5, divg and curl^ respectively. The gradient acts component-wise on tuples of 
scalar functions, whereas the divergence and curl act vector-wise on tuples of vector 
fields. In polar coordinates, the R^"'"'^ differential operators are 

Vxf = {drUr)n + r^-^VsUr, 
div,/ = r-"9,(r"(/,)x) +r-Miv5(/.)||, 
curlx/ = r^^UA (<9^(r(/r)||) - Vsifr)±) + r"^curl5(/^)||. 

We use the boundary function space L2(5'";V), writing the norm || • II2, of L2 
sections of the complex vector bundle 



over S^, where C™ is identified with the trivial vector bundle and TcS*" denotes the 
complexified tangent bundle of S*". The elements of this bundle are written in vector 

1 1 



form / 



[a /?] , and we write /x := «, /y := /3 for the normal component 



a 

and tangential part. Note that V is isomorphic to the trivial vector bundle C*^^"^"-*"*, 
when identifying scalar, i.e. C™- valued, functions and m-tuples of radial vector 
fields. More precisely, the isomorphism is V 9 [« /?] H- an + /3 G C'-"'^'^"^'", for 
a G and /3 G (Tc5'^)'". 

The differential operators on 5" can be seen as unbounded operators. We use 
D(y4), R(v4), N(y4) for the domain, range and null space respectively of unbounded 
operators. Then 

and its adjoint 

-div5 : ^2(5"; (TcST) L^S""; C™), 

with domains D(Vs) = W^{S'';C"') and D(div5) = {g e L2(^"; (Tc^'^)™) ; 
divg^f G L'^{S"'; C"^)} are closed unbounded operators with closed range. The con- 
dition g G R(div5) = N(Vs')"'" is that J^„ g{x)dx = so R(div5') is of codimension m 
in L2(5"; C"). Also when n > 2, R{Vs) = N(curls), and when n = 1, ^ G Ri^s) 
if and only if j^-^{g{x).,f)dx = 0. Thus R(V5) is of codimension m in L^(S'^;C'^) 
when n = 1, and infinite codimension when n > 2. 
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Definition 3.1. In L2(5'";V), we define operators D :- 



-I 
I 



where D{D) :- 



D(diV5) 







-dws 




Write A^+/ := i(/ + N)f 



and := 
? , and 



N-f:=\{I-N)f 



A basic observation is tliat tlie two operators D and anti-commute, i.e. 

ND = -DN. 

Of fundamental importance in tliis paper are the closed orthogonal subspaces 



n := R(Z^) 



■R(div5)' 

R(V5) 



and 



N(div5) 



We consistently denote by P-^ the orthogonal projection onto "H. We remark that 
'0' 
/„ 



• divs/ii = 



and 



constants being identified to constant functions. It can be checked that 
alent to A is strictly accretive on 

^1 := {g e 0(1+")"^) ; g, G R{Vs)}, 



IS equiv- 



More precisely, the accretivity assumption on A 



(13) 

uniformly for a.e. r G (0, 1) 
rewrites 

n m „ n m „ 

j J=0 a,/3=l j=0 a=l ^ ° 

for all (7 G "Hi, a.e. r G (0, 1). In fact, as we shall see in Lemma EH] this is equivalent 
to pointwise strict accretivity when n = 1 (unit disk), but this is in general not the 
case when n > 2 except if m = 1 (equations). 

Using the notation above, we can identify "Hi with 

R(V5) 

and see that "H is a subspace of codimension m in "Hi. 

On identifying 0^^+^^"" with V, the space of coefficients Loo(0^+"; £(C(i+")'")) 
identifies with Loo(0^'^"; C(y)), so that we can split any coefficients A as 

A^^{rx) Ax||(ra:) 
A„ I (rx) A„„(rx) 



A{rx) 



with A^^(rx) G CiC""; C"), A^^^{rx) G £((T,5")'^, C"), A^^^{rx) G £(0™, (T,^'^)™) 
and A||||(rx) G £((T,^")™, (T,^")'"). Note also that A^^{rx) = {A{rx)n,n). 

With our accretivity assumption ( |T4l) . the component ylxx(r-) seen as a multipli- 
cation operator is invertible on ^2(5'"; C"*), thus as a matrix function it is invertible 
in Loo(>S'"; C™). This is the reason why strict accretivity on "Hi is needed, and not 
only on "H, so that the transformed coefficient matrix A below can be formed in the 
next result. We make the above identification for coefficients A without mention. 

We can now state the two results on which our analysis stands. Proposition 13.31 re- 
formulates this Cauchy-Riemann system (fT2|) further, by solving for the r- derivatives. 
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as the vector-valued ODE ffT7|) for the conormal gradient / defined below. This for- 
mulation is well suited for the Neumann and regularity problems. For the Dirichlet 
problem, we use instead a similar first order system formulation of the equation; see 
Proposition 13. 5[ As explained in [H Sec. 3], the vector- valued potential v appearing 
there should be thought of as containing some generalized conjugate functions as 
tangential part. In the case of the unit disk, we make this rigorous in Section [5] and 
come back to this in Section [131 The fundamental object is the following. 

Definition 3.2. The conormal gradient of a weak solution u to diVxv4Vxtt = in 
0^+"- is the section / : R+ x S*" V defined by 

(15) ft = e-("+i)*/2 KAgrh 

I {9r)\\ 

where r = e~* and g = VxM- The map gr ^ ft is called the gradient-to-conormal 
gradient map. 

Proposition 3.3. The pointwise transformation 



All ] A , , A«« All I A , , Ai 



is a self-inverse bijective transformation of the set of bounded matrices which are 
strictly accretive on Hi. 

For a pair of coefficients A G L^{0^+''; £(0(1+")"^)) and B G L^{Yi+ x £(V)) 
which are strictly accretive on Hi and such that B = A, the gradient-to-conormal 
gradient map gives a one-to-one correspondence, with inverse the conormal gradient- 
to-gradient map 

(16) ft^9r = r-'^{{Bft)3+Ut)\^ 

wheret = ln(l/r), between solutions g G L]^^{0^~^'^; Q{^+^)"^'^ to the Cauchy-Riemann 
system ffT^ in O^^^ \ {0} distribution sense, and solutions f G L2°^(R_|_; "H), with 
Ii° WftWi^'t < oo, to the equation 

(17) dtf + {DB + ^N)f = 0, 
in R+ X S*" distributional sense. 

Recall that the Ricci curvature of S*" is n — 1, so the constant is related to 
curvature. On the other hand, the exponent appearing in the correspondence 
Or ^ ft is the only exponent for which no powers of r remain in equation ( fT71) . 
It turns out that this also makes the gradient-to-conormal gradient map an L2 
isomorphism since 

r rl poo 

(18) / |(?pdx^ / hrWlr'^dr^ / WftWldt. 
Jo^+" Jo Jo 

Proof. The stated properties of the matrix transformation are straightforward to 
verify, using the observation that e*^""*"^^* Re(i?t/t, /t) = Re{Argr, gr)- See [U Prop. 
4.1] for details. 

(i) Assume first that the equations fll2p hold on O^''"" \ {0}. In polar coordinates 
X = rx, the equations diVx(A5f) = 0, curlx(5') = give 

'r-^dr{r^{Ag)^) + r''diYs{Ag)^^ = 0, 
^r{rg^^) - Vsg± = 0. 
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Next we pull back the equations to R+ x S*". Write (Ag)^ = r and 
(Ag), = A,,g, + A,,g,. Then g, = r-("+i)/M:l(A - A,,f,) and g, = r-("+i)/2/ii, 
and the equations further become 

r-"9,(r("-i)/2/x) + r-(^^'^/Mws{B,J\ + B,J\) = 0, 
9,(r(i-")/2/|i) - r-("+i)/2V5(i?xx/x + fixii/ii) = 0. 
Using product rule for dr and the chain rule —rdr = dt, this yields the equation 

(HZD- 

It remains to check that ft ^ H for almost every t > 0. This is equivalent 
to {Argr)± G R(div5) and {gr)\\ G R(V5) for a.e. r e (0, 1). To see {Argr)± € 
R(div5) amounts to seeing that Jg„{Argr)±dx = 0. We apply Gauss' theorem as 
follows. For any radial function G C^(0^~'""'; C"*), the divergence equation gives 
jQi+n{Ag,'V(f))dyi = 0. Taking, for a.e. r G (0,1), the limit as </> approaches the 
characteristic function for balls {|x| < r} shows that j[^g„{Ag)j^dx = 0. To check 
{gr)\\ G R(V5) we distinguish first n = 1. In that similar application of 

Stokes' theorem shows that Jgi{T, gr)dx = for a.e. r G (0,1). For n > 2, that 
cuTls{{gr)\\) = is a consequence of cmlj^g = and the general fact that puUbacks 
and the exterior derivative commute. Hence ft G T-L. 

(ii) Conversely, assume that equation (1171) holds and ft & H for a.e. t > 0. 
Define the corresponding function g G L]^'^ (O^'^'^ ; C'^^"'""'-'™') by the conormal gradient- 
to-gradient map and note that cmls{{gr)\\) = 0. Reversing the rewriting of the 
equations in (i) shows that div^iAg) = 0,curlx(5') = hold on O^^" \ {0}. This 
proves the proposition. □ 

Corollary 3.4. For any coefficients A G Loo{0^~^'^; C{C^^~^'^^'^)) which are strictly 
accretive in the sense of ([2]), gradients of weak solutions to ([I]) in O^^"- are in one- 
to-one correspondence with R+ x S"' distributional solutions to the equation ( |T7I) . 
belonging to L2°'^(R_|_; 7/) with estimate ||/t||2'^^ < oo. 

Proof. Combine Proposition 12.11 and Proposition 13.31 □ 

There is a second way of constructing weak solutions that we now describe. 

Proposition 3.5. Let A and B = A be as in Proposition Assume that v G 
L^°'^(R+; D{D)) with \\Dvt\\ldt < oo satisfies 

(19) dtv + {BD - '-^N)v = 

in R+ X S*" distributional sense. Then 

Ur:=r 2 (t,j)^, r = e ' G (0, 1), 

extends to a weak solution of div^AV^x^u = in O"'^'''", and Dv equals the conormal 
gradient of u. 

Proof. By definition of u in the statement, (-Df)|| = V^Wx = r'^""'"^^/^(r~^Vsn). On 
the other hand, taking the normal component of dtV + [BD — ^^N)v = gives 

dtv^ - A-l{dwsV\i + A^^^Vsv^) + av^ = 0, 

or equivalently 

(Dv)^ = -div5W|| = -A^^{dt + (t)v^ + AxyV^Ux 

= r("+i)/2^Axx9.u + A^^^r-^Vsu) = r^^+^^/^^AVxu)^. 
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These equations hold in 0^+" \ {0}. Next, applying D to ffT^ yields 

{dt + DB + ^N){Dv) = 0. 

Thus / := Dv satisfies (fTTl) and ft G R(-D) = T-i- By Corollary 13.41 there is a 
weak solution u in O^"^*^ of the divergence form equation associated to /. In par- 
ticular, /|| = r^'^^^'>^'^{r~^'W su) and /x = r'^'^^^'^^'^{AV y{a) ^. Applying the conormal 
gradient-to-gradient map, we deduce Vx^t = Vx^i in O^^" \ {0} distribution sense. 
In particular, u = u + cm O^^^ \ {0} for some constant c. As n + c is also a weak 
solution in O^"*"" to the divergence form equation with coefficients A, this provides 
us with the desired extension for u. □ 

For perturbations A of radially independent coefficients, Corollary ll2.8( i) proves a 
converse of this result, i. e. the existence of such a vector- valued potential v containing 
a given solution u to diVxAVxW = as normal component. We do not know whether 
such V can be defined for general coefficients (except in O^, see Section [5]). 

Remark 3.6. Assume that the coefficients A are defined in R^"*"" and that the ac- 
cretivity condition ([2]) or holds for a.e r G (0, oo). As in Proposition [331 there is 
also a one-to-one correspondence between solutions g G L2°'^(R^"''"\ 0^+**; L2(5'"; V)) 
to divx(v4(7) = 0, curlx^? = in the exterior of the unit ball and solutions / : R_ — t- l-i 
to the equation dtf + {DB + ^N)f = for t < in L2(R_; Ti). Also, as in Proposi- 
tion[331 L^^'^-solutions v: R. ^ ^2(^''; D{D)) to the equation dtV+{BD-^N)v = 
for t < 0, give weak solutions u to diVxAVxW = in the exterior of the unit ball. 

4. Study of the infinitesimal generator 

In this section, we study the infinitesimal generators DBq + ^^^N and BqD — ^^^N 
for the vector-valued ODEs appearing in (fTTl) and ( |T9l) for radially independent 
coefficients -Bq = G L^^S"", CiV)), strictly accretive on l-L with constant k = 
> 0. Note that strict accretivity of Ai on "Hi is needed for the construction of 
Bq = Ai as a mutiplication operator. Once we have -Bq, only strict accretivity of Bq 
on "H is needed in our analysis. This has the following consequences used often in 
this work. First, Bo : H ^ B^H is an isomorphism. Second, the map P-hBq is an 
isomorphism of "H. 

The first operator will be used to get estimates of VxM, needed for the Neumann 
and regularity problems. The second operator will be used to get estimates of the 
potential m, needed for the Dirichlet problem. 

Definition 4.1. Let cr G R. Define the unbounded linear operators 

Dq := DBq + aN and Dq := BqD - aN 

in L2(^";V), with domains D{Dq) := Bo^D{D) and D(5o) := D{D) respectively. 
Here Bq^{X) := {/ G L2 ; BqJ G X}. When more convenient, we use the notation 
Da^ := Do and Da, := Dq. 

For these two operators, we have the following intertwining and duality relations. 

Lemma 4.2. In the sense of unbounded operators, we have DqD = DDq and 
{Da-J* = DB* -aN = -N{DA[ + aN)N . 

Proof. The proof is straightforward, using the identity B^ = NAlN for the second 
statement. □ 
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Proposition 4.3. In L2 = L2{S"';V), the operator Dq is a closed unbounded oper- 
ator with dense domain. There is a topological Hodge splitting 

i.e the projections Pj^^^ and Pg^^ onto H and Bq^T-L-^ in this splitting are bounded. 
The operator Dq leaves "H invariant, and the restricted operator Dq : ^ H, with 
domain D[Dq) fl H, is closed, densely defined, infective, onto, and has a compact 
inverse. 

If (T ^ 0, then Dq : L2 ^ L2 is also infective and onto, and Do|_b-i-h^ ~ 
If a = 0, then Dq = DBq, N{Dq) = Bq^H^ and R{Dq) = H are closed and 
invariant. In particular, when n = 1, dim N{Dq) = 2m = dim(L2/ R{Dq)). 

Proof. The splitting is a consequence of the strict accretivity of Bq on "H, and it is 
clear that "H is invariant under Dq. Note that 

{iN){DBQ + aN) = {iND)BQ + za, 

where iN is unitary on L2 as well as "H, and where iND = —iDN is a self-adjoint 
operator with range "H. This shows that Dq is closed, densely defined, injective and 
onto on "H, and on L2 when cr 7^ 0, as a consequence of properties of operators such 
as {iND)BQ stated in [71 Prop. 3.3]. 

Next we show that Dq : H H has a compact inverse. Write Dq = D{P-^Bq) + 
aN. Since P-uBq is an isomorphism on "H, it suffices to prove that the inverse of 
D ■.n^n is compact. Note that D{Vs) = W^iS""-, C™) is compactly embedded in 
L2{S^; C™) by Rellich's theorem. In particular V5 : R(div5) — >■ R(V5) has compact 
inverse. Since V5 = — div^, it follows that divg : R(Vs') — )■ R(div5') has a compact 
inverse as well. This proves that the inverse of D is compact on "H. 

The remaining properties when cr 7^ and cr = are straightforward and are left 
to the reader. □ 

Proposition 4.4. In L2 = 1/2(5""; V), the operator Dq is a closed unbounded oper- 
ator with dense domain. There is a topological Hodge splitting 

L2 = Bq'H®'H^, 

i. e the projections P}^^ and P^^ onto BqT-L and H"*" in this splitting are bounded. Here 
Ti^ C D{Dq) and Dq leaves Ti^ invariant. 

If a ^ 0, then Dq : L2 —> L2 is also injective and onto, and Dq\-^± = —aN . 

If a = , then Dq = BqD, N{Dq) = "H-*- and R{Dq) = BqT-L is closed. In particular, 
the subspace BqT-L is invariant under Dq and when n = 1, dim N{Dq) = 2m = 
dim(L2//?(5o)). 

Proof. These results for Dq follow from Proposition 14. 31 by duality, using Lemma 1421 

□ 

Remark 4.5. The reader familiar with [HI H] should carefully note the following 
fundamental difference between the cases a ^ and a = 0. When a = 0, each 
of the operators Dq and Dq is of the type considered in [HI H] , and each has two 
complementary invariant subspaces. On the other hand when a 7^ 0, the operator 
Dq has in general only the invariant subspace "H, and Dq only has the invariant 
subspace "H-*-. One can define an induced operator Dq on the quotient space L2/'H'^, 
but this cannot be realized as an action in a subspace complementary to "H"*" in L2 
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in general. As a will be set to this means for us a difference in the treatment of 
n = 1 {i.e. space dimension 2) and n > 2 {i.e. space dimension 3 and higher). 

We prove here a technical lemma for later use. 

Lemma 4.6. There is a unique isomorphism 

(20) H L2/H^ -.hy-^h 
such that DqH = Dh for h e U <r] D{Do). 

Proof. When a = 0, we can take h := BqH G BqH ^ L2/'H'^ as DqH = DBqH = Dh. 

When 0" 7^ 0, we use that D : L2 — t- "H is surjective with null space "H"*". This 
defines /i for /i e "H fl D{Do). With the compact inverse of D : Ti ^ Ti, the 
equation Doh = Dh is equivalent to 

(21) PnBoh + aD-^Nh = P^h. 

This shows that ( 1201) extends to a bounded map since ||/;.||i2/wi ~ ll-Pw^lb- More- 
over, since P-hBq is an isomorphism on "H, we have also the lower bound \\h\\2 < 
||P-^i?o/i||2 ^ II^IIl2/w-l + ||-D~^/i||2, which shows that ([20]) is a semi-Fredholm op- 
erator. If /i = 0, then dni) implies h eUr\ D{Do). Therefore Doh = and (1^ 
is injective. Since the range contains the dense subspace D{D)/T-L'^, invertibility 
follows. □ 



5. Elliptic systems in the unit disk 

In dimension n = 1, i.e. for the unit disk C with boundary S^, some 
special phenomena occurs. In this section we collect these results. 

Lemma 5.1. If n = 1 and A is strictly accretive in the sense of ([2D, then A is 
pointwise strictly accretive, i.e. 

Re{A{x.)v,v) > for all v e C^*^, and a.e. x e O^. 

Proof. By scaling and continuity, it suffices to consider v = [{za) {wa)]* G C^™, 
with Wa 7^ 0, a = 1, . . . ,m. In ([2]), let 

u"(re*^) := {iky^Wae"'^^ri{e'^)e"'^ , a = 1, . . . , m, 

with a smooth function 77 : S*^ — t- R, A; G Z_|_ and tq G (0, 1). Using polar coordinates 
and letting /c — )■ 00 yields 

Re / {A{rQx)v,v) \ri{x)\'^dx > kIv]"^ / \ri{x)\'^dx, for a.e. roG(0,l). 

Taking |?7p to be an approximation to the identity at a given point x ^ now 
proves the pointwise strict accretivity in the statement. □ 

Definition 5.2. Assume that A G Loo(0^; £(C^™)) is pointwise strictly accretive. 
Given a weak solution u G 14^2^' °^(0^; C™) to divxAVx^i = 0, we say that a solution 



u G W2'^"^{0'^', C"*) to .J'V:^u = AVx^i is a conjugate of u where J :- 



-/ 

1 
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We note that since AVx^i is divergence- free, there always exists a conjugate of u, 
unique modulo constants in C"*. The notion of conjugate solution for two dimen- 
sional divergence form equations, in the scalar case m = 1, goes back to Morrey. See 
[35]. Note that when A = I, the system JV^u = VxM is the anti Cauchy-Riemann 
equations. 

Lemma 5.3. Assume that A G Loo(0^; £(C^™)) is pointwise strictly accretive. Let 
u G 1V2^''°'^(0^; C") he a weak solution to diVxAVx^ = 0. Then 



AV^u = JVy-U ^ 



(AVxM) 

(AV^u) 



-(Vxm) 

(VxM)ii 



AVjcU = fV^u ^ divxAVxM = 



where A is the conjugate coefficient defined by 

A := J'A-^J. 

We have 

^_ {d — ca^^b)^^ {d — ca~^b)^^ca~^ 

a~^b{d — ca^^b)^^ + a~^b{d — ca~^b)~^ca~^ 

When m = 1, this reduces to A = (det A)~'^A'^. 



if A 



a b 
c d 



Here, we have identified the tangential part (■)|| with its component along r. (See 
below.) 

Proof. The equivalences and implication are verified from A = J^A~^J. The explicit 
formula for A is classical if m = 1. If m > 2, the proposed formula for A can be 
checked by a straightforward computation. Note that a,b,c,d G Loo(0^; £(0™)) 
and all the entries of A as well: the inverses are pointwise multiplications. We omit 
further details. □ 

We next show that the vector- valued potential v in Proposition 13.51 contains, along 
with u as normal component, its conjugate u as tangential component. To do that, it 
is convenient to identify V with the trivial bundle C^"* by identifying the tangential 
component (3 to the tangential part /3r G (TcS'^)™. 

Given this identification, D becomes 



D 



■ -df 
& 



where 9^ denotes the tangential counter clockwise derivative of m-tuples of scalar 
functions on S^. A coefficient A G Loo(0^; £(C^™')) is thus identified with its ma- 
trix representation in the moving frame {n,f}. We remark that this identification 
commutes with the matrix J. 



Proposition 5.4. Let A G Loo(0^; £(C^™)) be pointwise strictly accretive and let 
B := A e Loo(02;/:(C2"*)). Assume that v = [u m]* G 4°^(R+; D(L))) with 
\\Dvt\\ldt < CO is a R+ x S"- distributional solution to dtv + BDv = as in 
Proposition \3.5\. so that Ur = {vt)±, r = e~*, is a weak solution to diVxAVxM = in 
O^. Then u is a conjugate to u. 

Conversely, given a weak solution u to diVxAVxW = m and a conjugate u, 
the potential vector v = \u tt] has the above properties. 
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Note that the construction of v this way is a feature of two-dimensional systems 
as compared to higher dimensions. 

Proof. Applying J* to dtv + BDv = gives dt{Jh!) + BD{J^v) = with B = J^BJ, 

since JD = DJ. A calculation shows that B = A. Applying Proposition 13.51 shows 
that Ur = {J^Vt)± is a weak solution to diVxAVxW = 0. Also we know that Dv and 
Dv are respectively equal to the conormal gradients of u and u, and since J^v = v, 
this gives the middle term in the equivalence of Lemma 15.31 Thus £t is a conjugate 
of u. The converse is immediate to check and left to the reader. □ 

We finish this section with the following simple expressions for the projections P^^ 
and P^^ of Propositions 14.31 and 14.41 when n = 1. We still make the identification 

Lemma 5.5. Let Ai G Loo(0^; £(C^™)) be pointwise strictly accretive radially in- 
dependent coefficients, and let Bq := Ai G Loo(0^; £(C^"^)) the corresponding coef- 
ficients. Then 

P^g = fio"' dx^ B,^g dx, g G L2{S'; C^™), 

andPl = B,'P'^^Bo. 

Proof. By accretivity, (f^^ Bq^ dx)~^ is a bounded operator (called the harmonic 
mean of Bq). If g & BqH, then BQ^g G "H and J^i B^^gdx = 0, hence Pbo9 = 0; 
follows. On the other hand, if g E 'H'^, then g is constant, and therefore the right 
hand side equals 



(^^^ B^' dx^ (^^^ B,' dx^ g = g. 



This proves the expression for P^^. The formula for P^^ comes from the similarity 
relation DBq = Bq^{BqD)Bo. □ 

6. Resolvent estimates 

In this section we prove that the spectra of Dq and Dq are contained in certain 
double hyperbolic regions, and we estimate the resolvents. For parameters < w < 
z/ < 7r/2 and cr G R, define closed and open hyperbolic regions in the complex plane 
by 



QO 



{x + G C ; (tan^ uj)x^ >y+ cr }, 

{x + iy e C ; (tan^ z/)x^ > y"^ + o-^}, 

{x-^iy eC ; (tanw)x > (y^ + a^Y^"^}, 

{x + iy eC ; (tan u)x > (y^ + a^Y^^}. 



When 0" = 0, we drop the subscript a in the notation for the sectorial regions. 

Proposition 6.1. On L2 = L2{S"';V), there is a constant u G (0,7r/2), depending 
only on ||i?o||oo ond the accretivity constant Hbq, such that the spectra of the operators 
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Dq and Dq are contained in the double hyperbolic region S^^^r- Moreover, there are 
resolvent bounds 

||(A - Do)'^\\l^_^l^, ||(A - -Do)""^||l2->L2 < / „ „ - 

yy"^ + a^/ tanw — \x\ 

for all X = x + iy ^ e same estimates hold for the restriction Dq : 1-L ^ %. 

Proof, (i) To prove the spectral estimates for Dq? assume that 

{DBq + aN-x- iy)u = f. 

Introduce the auxihary operator Ny := iaN ~yl, and note that \\Ny\\ = \\Ny^\\^^ = 
A/y^~-i-~a^. Multiply with A^^ and rewrite as 

(22) {NyD)Bou + i{y^ + (t^)u = Nyf + xNyU. 
Now split the function u as 

u = ui + uo en^B^^n-^, 

and note that \\u\\ ~ + ||no||. Apply the associated bounded projections P^^ to 
(1221) to get 

{NyD)BoUi + liy" + g'')u^ = P^^Nyf + xP^^NyU, 

+ t{y^ + a^)uo = P^Nyf + xP^^^NyU. 

Take the imaginary part of the inner product between the first equation and BqUi 
(using that NyD is self-adjoint), and the second equation and Uq to get 

(y^ + a^) Re(Mi, BqU^) = Im (P^iV^/, Po^i) + Im {xP^^NyU, BqUi), 

{y^ + a')\\uof = lm{P's^Nyf,Uo) + lm{xPlNyU,Uo). 

Using the strict accretivity of Bq on gives the estimate 

{y' + c')\\ur < Civ/^^T^(||/||||«|| + IxlWuf), 

for some constant Ci < oo. Thus < {\/y'^ + cr^/Ci — |x|)~"'^||/||. 

(ii) To prove a similar lower bound on Dq, assume that [BqD — aN — x — iy)u = /, 
and rewrite as 

(23) BoDNy^NyU + iNyU = f + xu. 

Write NyU = BqUi + Uq E BqH © "H^. Apply the bounded projections P^^ to ( 1^ 
to get 

BoiDN.;')BoUi + tBoUi = P^ + xP^^u, 

Q + iuo = Plj + xPl^^u. 

Recall that Bq : ^ BqH is an isomorphism and apply its inverse Bq^ : BqH — t- Ti 
to the first equation. Then take the imaginary part of the inner product between the 
first equation and BqUi (using that DN~^ is self-adjoint), and the second equation 
and Mo to get 

Re(Mi, Boui) = Im {B^^Phj, BqUi) + Im {xB^^Pl^^u, BqUi), 
||nof = Im {PbJ, Uo) + Im {xP^^u, uq). 
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Using the strict accretivity of Bq on T-L gives the estimate 

iy' + a')\\ur < C2i\x\\\u\\ + Wmy' + a'Y/^u\ 



for some constant C2 < 00. Thus HmH < {^/y'^ + o"^ jC^ — ^\\f\\- 

(iii) Using that DBq + and B^D + aN are adjoint operators, combining the 
results in (i) and (ii) shows that both operators D^ — \ and D^ — \ are onto, with 
bounded inverse, when A ^ 5^,0-- Here uj := arctan(max(Ci, C2)). The estimates on 

n follow. ' □ 



< e, dosed sets E,F C S" 



We shall also need the following off-diagonal estimates for the resolvents, both in 
L2 and in Lp for p near 2. 

Lemma 6.2. (i) There exist e, a > such that for — ■ 
and f G Lp{S"'; V) with supp f C E and t G R, 

where d{E,F) is the distance between the sets E and F. 

(ii) There exist q > 2 with j — -<e, and a > such that for closed sets E,F C S" 
and f G L2(S'"; V) with supp f <Z E and = and |t| < 1, 



[I + ^tDor'f\\LJF)<\t[ 



n{^-^)^-ad{E,F)/\t\ 



Proof. We first prove (i). The case p = 2 follows the argument in [HI Prop. 5.1]. It 
remains to prove Lp boundedness for p near 2 as, the Lg off-diagonal bounds follow 
by interpolation with the L2 off-diagonal bounds for q between p and 2. 

For f E LpH L2, we let h = {I + itDo)~^f and wish to prove \\h\\p < when p 
is near 2 and uniformly in t. To prove this, we rewrite the equation (J + itDQ)h = f 
first as (/ -|- itaN + itDB^jh = f and then in terms of a divergence form equation, 

with coefficients Ai = Bq. Write h = [{Aih)j_ /i,,]* and / = [{AJ)^ f^]\ Then 

(1 - zta){Aih)^ - ztdws{Aih)n = (Aj)^, 
(1 + zt(T)/i|| + itV shi_ = /||. 

Using the second equation to eliminate /ly in the first equation, and letting z = 
(1 + ita)~^, we obtain 



with 



[1 —itzdivs] 
L := [1 -itzdivs] Aj 



z{A,),J, + iz-z)iA^),J, 
-^(^1)1111/11 



[1 -irdiv^] Ae 



1 

-irVsj 

\tz\, and Dg the diagonal matrix with 



1 

-itz'Vs 

and Ae = D^gAiDe with tz = e*^r, r 
entries 1, e*^ in the normal/tangential splitting. We note that Ae is strictly accretive 
on "Hi with the same constants as Ai, and that \z\ < 1 and |r| < |cr|~-^. We claim 
that L is invertible from the Sobolev space (S""; C") equipped with the scaled 
norm 

i/p 

(24) 



{\u{x)\^ + \TVsu{x)\''Yl^dx] 



5" 



to its dual, with bounds independent of r, 6, for p in a neighborhood of 2. 
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To prove this, if we rescale from the sphere 5" of radius 1 to the sphere 5"^^ 
of radius 1/r, we obtain the same equation with A0,Ai,z±. unchanged, /(x), h{x) 
replaced by fijx), hijx), and rdiv5, tV s replaced by div^n^ , Vs^^ , and we want to 
show ||/i(T-) 11^^(5" ) < ||/(t-) lUpCsj"/ ) (with implicit constant uniform in t,6). Thus 
it is enough to set r = 1 and work on S*", as long as we only use estimates on 
which hold (with same constant) on S'Jy^ as well. 

Having set r = 1, we have, for 1 < p, g < oo such that 1/p + 1/g = 1, estimates 

< PelloollMlliypi = ||Ai||oo||M||iypi, 
where ||tt||u/-i •= supn m i |(it,f)| and {u,v) denotes the L2(S'"; C") pairing ex- 

P 1 1 1 1 

tended in the sense of distributions. For p = q = 2, the accretivity assumption on 
Ae yields ||Lm||^-i > Applying the extrapolation result of Sneiberg [38] to 

the complex interpolation scale {VT^ }i<p<oo, shows the existence of e > such that 

for 1^ — || < e. (Even for r 7^ 1, one can verify that the Sobolev norms given by ( 1241) 
on S'ly^ (with rVs replaced Vsj^^ ) are equivalent to the ones given by the complex 
interpolation method, with constant independent of r. Hence e depends only on the 
ellipticity constants and dimension, and is thus independent of r, 6.) Applying this 
isomorphism, we obtain the resolvent estimate 

ll^llp ~ ll^xllp+ ll^iillp ^ \\h±\\w^ + ll^illp < ll/llp- 

In the second step we used /ly = — ze'^Vs/ix + ^/n and \z\ < 1 (recall we have rescaled 
and set r = 1), and in the third step that [l -idiv^] : Lp{S''; C(i+")™) ^ W'^ is 
an isometry since [l — iVs]* : Wg — ?■ Lq{S^; Q(^+^)^"j is §0. This finishes the proof 
of «• 

To prove the inequality (ii), the above argument shows that Lh^ = z^{Ai)^^f^ 
and = —itz-^.'Vshj_. Having rescaled in the same way, the Sobolev embedding 
1/2 C for some q > 2 with \ — ^ < ^, allows us to conclude that h G Lq{S^; V) 

and since we assume |t| < 1, we have |r| ~ \t\ and obtain \\h\\q < |t[ ||/||2, the 

power coming from scaling. It suffices to interpolate again with the L2 off-diagonal 
decay, and conclude for any exponent between 2 and q. □ 

We state the following useful corollary. Here and subsequently, is defined as 
N^, replacing averages by Lp averages and M is the Hardy-Littlewood maximal 
operator. 

Corollary 6.3. For the e as above and \ ^ — \\ < e, we- have the pointwise inequalities 

N^{{I + ztDo)-'f)<M{\fn^^ 

N^{{i+ztDor'f)<M{\fn'/p 

and for some p < 2 with ~ — | < 

iv,(((/+^t5o)-v)x)<M(i/nv.. 
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Proof. We fix a Whitney region Wq = W{to, Xq) in R+ x 5". Then \Wo\-^ J^^ \{I + 
itDo)~^ f{x)\Pdtdx < M(|/|^')(xo) follows directly from the off-diagonal decay of 
LemmaOas in [SI Prop. 2.56]. Next, \Wo\~^ Jy^^ \{I+itDoy^ f{x)\Pdtdx < M(|/|p)(xo) 
follows by testing against g G Lg{WQ; V), supported in Wq with 1/p + 1/q = 1. We 
have 

/ ((J + ttDo)-^fix),g{t, x))dtdx = r\f,{I- itDl)-^gt)dt 

JWo ^ ^ Jto/co 

SO that for each fixed t, using that Z)q = DBq — aN has the same form as Do, we 
can use the Lg off-diagonal decay for each t ~ to and obtain for any M > 0, 

(25) iW^ol"^ / \{I + itDo)-^f{x)fdtdx 

Jwo 

<J2^-^''\B{xo,2Ho)\-' [ \f{x)\^dx, 

using standard computations on annuli around i?(xo,to) iii "S*". Details are left to 
the reader. 

The last estimate starts in the same way with g G L2{Wo\ V), but since we want 
to estimate the normal component of (/ -|- itDo)~^f we assume that {gt)\\ = for 
each t. The second estimate in Lemma [6.21 implies that (/ — itDQ)~^gt =: ht has Lq 
estimates with decay. Thus using Holder's inequality on J^^/'^if, ht)dt with exponent 
q on ht and dual exponent on / yields 

1/2 



(26) \\Wot^ [ \{{I + ttDo)-'f)Ax)\''dtdx) 



Wo J 

<Y,^-''' (mxo,2HoT' [ IfixWdx' 

j>2 V JB(xo, 



xo,2Ho) 



and the conclusion follows. □ 

7. Square function estimates and functional calculus 

All the remainder of this article rests on the square function estimate below. 

Theorem 7.1. Let n > 1. The operator Do = DBo + (tN, with cr G R fixed but 
arbitrary, has square function estimates 

r°° dt 

/ \\tDo{l + eDlr'f\\l-^\\f\\l forallfeRiDo). 
Jo 

The estimate < holds for all f G L2{S"', C"). The same estimates hold for Do = 
BoD - aN. 



Proof. Note that the equivalence can only hold on R(-Do) = R(-Do) which equals 
-^2(5'"; V) if cr 7^ and "H if a = 0. By standard duality arguments, the estimates > 
on R(-Do) follows from the estimates < for Dq. See pQ. Note that Dq is of type Do- 
Hence it is enough to prove 

(27) r\\tDo{l + t'Dl)-'f\\l^<\\f\\l 
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for all / G L2{S"';V), and similarly for Dq. Consider first the operator Dq. 
(i) We first reduce to 



(2^ 



[ \\WBo{l+t\DBoY)-' 
Jo 



f\\lj< 



for all / G ^2(5"; V). First note that 



2 (it 



using that -Dq has bounded inverse by Proposition 14.31 (When n = 1, write / = 
/i + /o £ T-i® Bq^I-L^ . The above estimate goes through for /i, and the contribution 
from /o is zero.) For the integral jj^, we may ignore the zero order term in Dq, using 
the idea from [6l Sec. 9]. Indeed, 

\\{I + itDor'f-{I + itDB,r^fh = W + itD^yHtaNil + itDB^r^fh < \t 

Since 2itDo{I + f^D^)-^ = {I - UDo 
subtraction yields 

»i 



^ — (/ + itDo) ^, and similarly for DBq, 

dt 



f\\tD,{I + eDl)-'f\\f-^< f\\tDB,{I + t\DB,f)-'f\\l^+ f 
Jo Jo ^0 



tdtwfwi 



(ii) Next, using a partition of unity, it suffices to show that 



(29) 



j\mB,{i+t\DB,fr'f\\fj< 



when C is a smooth cutoff that is 1 on a neighborhood of supp /. Indeed, L2-off 
diagonal estimates of tDi?o(l + t^(-D-Bo)^)~^ from Lemma 16^ and again 2itDBo[I + 
t'^{DBof)-^ = {I - itDBo)'^ - {I + itDBo)'^ show in this case that 

- OtDBoil + t'iDBorr'f\\l<t'\\f\\l- 

(iii) To prove fpUj) . we assume that / and ( are supported inside the lower hemi- 
sphere which we parametrize by O"^ using stereographic coordinates, i.e. using the 
map 



\y\ 



eo 



22/ 



|y|2 + r 



where eo G R^"*""" is a fixed unit normal vector to R" C R^^", which covers all S"", 
except the north pole eo G S*". Note that p is a conformal map with length dilation 
d~^ and Jacobian determinant dx/dy = d~^, where 

%):=(|yP + l)/2. 

Let T : R" ^ Ri+" : y ^ dyp{y) be the differential of p, and note that T*T = d-^I. 
Define adjoint rescaled puUbacks and pushforwards 

p* : L2(p(0"); V) ^ L2(0- C^^^-)-^) : [f\ /,,]* ^ [^-"(A o p) T\f, o p)]\ 



p, :L2(0";C(i+")^ 
Note that (p*)"^ = 



d"" 




L2(p(0");V): [g^ g,]' ^ [{g, 
p*. We claim that 



(d"T^7„)op-i]*. 



p*D = D, 









where -D 







-div„ 
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Indeed, the tangential part of the equation is the chain rule, and the normal com- 
ponent is the adjoint statement. We consider Dp as a self-adjoint closed unbounded 

operator in ^2(0"; C*^'^"''"^™') with domain D(Dp) 



ifnHO";C™) 



, where Hi de- 



D(div, 

notes the Sobolev 14^2^ functions vanishing at the boundary 5""^. 

Next we map coefficients Bq in p(0"') to coefficients Bp := (p^,)~^i?o(p*)^^ in O*^, 
and claim that Bp is strictly accretive on R{Dp). To see this, let g G R{Dp). Then 
curlj^^fy = and g^^ is normal on 90" (or if n = 1 we have j_^g^^dy = 0). Writing 
g = p*f and extending / by outside ^(0*^), it follows that f E Hi. (To see 
this, write = VyU with u G i/o(0"; C"'), and extend m by to an //^(R"; C")- 
function.) The assumed strict accretivity of Bq on "Hi gives 



Re/ {Bpg,g)dy = Re f {{p.)-'Bof,p*f)dy 

Jo" JO" 



= Re/ iBofJ)dx>K \f\'dx^ \g\'dy. 

J S" J S" JO" 

Thus we obtain a bisectorial operator DpBp in ^2(0"; C*^^+")™), and we observe the 
intertwining relation 

'd"" 

rf2-n 



p*DBof = Dp 



p*p,Bpp*f = DpBpp*f. 



for / supported in the lower hemisphere. In O", let := {|?/| < 1/4}. By rotational 
invarince, it is enough to consider those / = {p*)~^g with g supported on K and 
( = {p*)~^r] = 770 with ?7 G C^(R") be such that 77 = 1 on {\y\ < 1/2} and 
suppr/ C {\y\ < 3/4}. Using rjg = g and understanding 77, C as the operators of 
pointwise multiplication by 77, (, one can check the identity 

C(/ + itDB^y'f - {p*rW{I + itDpBpY^g 

= C(/ + itDB,Y\p*Y\r^g) - CipT^I + ^tDpBpr'g 

= C{I + itDBo)-\pT' Hi + ^tDpBp) -p*{I + ttDBo){pT'v) (I + ^tDpBpY^g 

= C(/ + itDBQY\p*)-Ht[7], Dp\Bp{I + itDpBpY^g. 

As in (i) above, subtracting the corresponding equation with t replaced by — t, yields 
the estimate 

\\CtDB,{I + ^{DBofr^f - {p*r^rftDpBp{I + t\DpBpfy'gh < \t\\\g\\2, 

since [?7, -Dp] is bounded. As ||/||2 ~ ||5'||2 by the support conditions, ( 129|) will follow 
from 

/ \\tDpBp{I + t'iDpBpfr'gWl- < Ml for all g G L2(0"; c(^+")'"). 
Jo 

(iv) The latter square function estimate follows from combining [T31 Thm. 2] and 
[m Prop. 3.1(iii)], the latter purely being of functional analytic content. (See P, 
Sec. 10.1] where this is pointed out.) 

(v) Consider now Do- Similarly one can reduce to prove < for BqD. On M{BoD) = 
7/-*", this is trivial. On R{BqD) = BqH we use that BqD is similar to DBq on 
R{DBq) = Ti through the isomorphism Bq : Ti ^ BqH. Thus the square function 
upper estimate for BqD follows by similarity from the one for DBq. □ 
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The square function estimates from Theorem 17.11 provide bounds on the S°^- 

holomorphic functional calculus of the operators Dq and Do, adapting the techniques 
described in p. Write 

H{Sl^„) := {holomorphic b : ^ C}, 

i/oo(^.V) ■■= e HiS:^^) ; sup{|KA)| ; A e S^J < oo}, 

vl>(^° J := {b G J ' IKA)| < min(|Ar, {XD. for some a > 0}. 

We summarize the result for the 5*° ^-holomorphic functional calculus in the fol- 
lowing corollary. The proof is a straightforward adaption of the results in [1]. 

Corollary 7.2. Assume a G R and Dq = DBq + aN. Fix cu < u < 7r/2. There is 
a unique continuous Banach algebra homomorphism 

H^{S:^^)^C{R{D,)):b^b{D,), 

with bounds ||6(£)o)/||2 < C'lsup^o^ l^(A)|)||/||2 for all f G R{Do), where C only 
depends on ||i?o||oo; f^Bo, n and a, and with the following two properties. If b & 
then 

b{D,) = ^ [ 6(A)(A - Do)-'dX G /:(/?po)), 

where 7 := dSo^„, u < 6 < u, oriented counter clockwise around S^^^^. For any 
b G Hoo{S°^) we have strong convergence 

lim ||6fcpo)/ - K^o)/||2 = 0, for each f G /?po), 

whenever bk G \l/(S'°g.), /c = 1, 2, . . are uniformly bounded, i.e. sup^j^ I^A:(A)| < 00, 
and converges pointwise to b. 

The corresponding results hold for Dq = BqD — aN replacing Dq by throughout. 

We remark that the square function estimates in Theorem 17.11 hold when ip{z) = 
z{l + z^^~^ is replaced by any G ^1/(5'°) which is non-zero on both components of 
We have 

/■"^ dt 
(30) / ||^(t/^o)/||^- ^ ll/lli for all / G Rpo). 

Jo 

A similar extension of the square function estimates holds for Dq. 
Fundamental operators in this paper are the following. 

Definition 7.3. (i) Let x^(A) and x~(A) be the characteristic functions for the right 
and left half planes. Define spectral projections Eq := x^{Dq) and Eq := x^{Dq) 
on R(-Do) and R(-Do) respectively. 

(ii) Define closed and dense defined operators A = iZ^ol '■= sgn(-Do)-Do and A = 
\Dq\ :=sgn(5o)5o onL2(^";V). Here |A| := Asgn(A) and sgn(A) := x+(A) - x~(A). 

Define operators e~*^ and e"*'^ on R(-Do) and R(-Do) respectively by applying 
Corollary O with 6(A) = e"*!^!, t > 0. 

When cr = 0, R(-Do) = Bo^^ = -BoR(-Do) are strict subspaces of L2 and it is 
convenient to extend the above operators to all L2. Using the Hodge splitting 
L2 = BqH © "H-*-, on "H-*- the operator Do = BqD is already and A = \BqD\ 
is naturally defined by 0. Using the other Hodge splitting L2 = H ® Bq^T-L-^, on 
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Bq ^"H-*- the operator Dq = DBq is already and A = |-D-Bo| is naturally defined by 
0. It follows that e~*^ and e~^^ are naturally extended to L2, by letting e~*^|^± := / 
and e-*^|^-i^x := /• 

However, for the projections the extension is more subtle. Indeed, we see for the 
functional calculus of -Do = DBq — aN that 

'6(a)/ 



biDo) = bi-aN) 



b{-a)I 



on 7i when a ^ using the definition of A^. As we are mainly interested in a = 

it is more natural for consistency of notation towards applications to divergence form 

equations to define the operators for a = by continuity a — )■ O"*". Thus set 

'6(0+)J 



biBoD) := biBoDU,n)Pi 



pO 



6(0-)/ 

where 6(0±) := lim±Ae5t^+,A->o 6(A), assuming the limits exist, 6(i?o-D|_Bo«) is the 



operator from Corollary 17.21 and 



So' 



0, 1, denote the projections from Proposi- 



tion 14.41 onto the subspaces in the Hodge splitting L2 - 
Similarly, for cr 7^ 0, we have -Do = DBo + aN so -Dq 
set 



BqH®H^. 
-- aN on B^^n^ 



For (J = 0, 



biDBo) ■.= b{DBo\n)Pj 



Bo 



pO 



biO-)I 






6(0+)/ 



where P'^^, i = 0,1, denote the projections from Proposition 14.31 onto the subspaces 



in the Hodge splitting L2 = Ti® Bq ^"H-*-. Remark that Pg^ on the left of the matrix 
is needed to obtain an element in Bq^T-L^. An elementary calculation shows that 
this extension of the functional calculus coincides with (6(/?q-D))* with 6(A) = 6(A) 
and that the extended functional calculi of -Do and -Do thus obtained are intertwined 
hjD. 

Taking 6(A) = A or Asgn(A), this provides us with the zero extension that we 
already chose so this is consistent. For the projections, this leads to the following 
definition. 

Definition 7.4. When cr = 0, extend Eq, Eq originally defined on R{BqD) = BqH 
and R{DBq) = "H respectively from Definition 17.31 to operators on all L2{S"';V), 
letting 

^E^f:=mf, forall/e-H^, 
E^f := P^^N^f, for all / G B^'U^. 

Lemma 7.5. With L2 = L2{S^; V), the spectral projections Eq and E^ are hounded, 
we have topological spectral splittings 

L2 

restricting to = EqH © EqH in the suhspace H invariant under Dq, and 

L2 



Eq L2 



Eq L2, 



E^L2 



Eq L2, 



restricting to V. = EqU © Eq Ti in the suhspace T-L invariant under Dq . We 
also have the intertwining relation 



(31) 



E^D 



DE^ 
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SO that D : EqL2 — )■ EqU. is surjective. 



If o' then in the latter splitting we have Eq = in T-L . Hence E^T-L = 
N^H^ and EqH^ = N+H^. (On the other hand, ifcr<0, then E^ = m H^.) 

Proof. When o" 7^ 0, R(-Do) = -^2 and L2 = R(-Do) by Proposition 16.11 Boundedness 
on L2 follows from Corollary I7.2[ The intertwining property is a consequence of 
Lemma 14.21 The surjectivity of D easily follows from the spectral subspaces and 
using D : L2 ^ T-L surjective and the splittings. That E^ = in "H-*- when a > 
comes from Dq = —aN in "H-*- and x^{~^^) = The case cr = follows from 
Definition 17.41 We leave further details to the reader. □ 

8. A DETOUR TO KATO'S SQUARE ROOT ON LiPSCHITZ SURFACES 

Let S be a surface in R^"*"*^, assumed to be Lipschitz diffeomorphic to S*" through 
a bilipschitz map po ^ "S*" — S. Let da denote surface measure on S. Consider, 
for n,m > 1, coefficient matrices H G Loo(S; £((TcE)™')) (with TcS denoting 
the complexified tangent bundle) and h G Loo(S; £(0"^)), assumed to be strictly 
accretive in the sense that 



for all u G W^iJl; C™) and z G C"*, and some k > 0. Then L := -divsi/Vs, with 
divs := — (Ve)* in L2(Tj;da), constructed by the method of sesquilinear forms, is a 
maximal accretive operator and hL is defined on D(L) and can be shown to be an 
CJ-sectorial operator on L2(T,;da) for some < u; < tt. Thus it has a square root 
and we have 

Theorem 8.1. The square root of the operator hL = — Mivi/Vs has domain 
D{^/hL) = W2 (T,; C"") , and estimates \\y/hLu\\2 ~ ||Vem||2- 

In particular for h = 1, we obtain a version of the Kato square root problem 
on Lipschitz surfaces E. The presence of h makes the theorem invariant under 
bilipschitz changes of variables as we shall see in the proof. 

Our Theorems 17.11 and 18.11 are inspired by [TU Thm. 7.1], and a comparison of 
these two results is in order. The main novelty in Theorems 17.11 and 18.11 is that 
these do not require the coefficients Bq or H to be pointwise strictly accretive, 
which was needed for the localization argument in [TU Thm. 7.1]. This theorem 
considered more general forms on S, and more general compact Lipschitz surfaces 
S. It is straightforward to extend our results Theorems 17.11 and 18.11 here to more 
general compact Lipschitz manifolds. On the other hand, we do not know how to 
extend our localization argument here to the case of forms, unless pointwise strict 
accretivity is assumed. 

We also mention that in his PhD thesis [36], A. Morris proved similar results on 
embedded (possibly non-compact) Riemannian manifolds with bounds on the second 
fundamental form and a lower bound on Ricci curvature. 




Proof of Theorem 




{hdiYj,HVj:u){po{x)) = {hdiysHVs{uopo)){x), 



X G 
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where/i(x) = |J(po)(a;)| ^h{po{x)) and H{x) := \J{po){x)\{po{x)) ^ H {po{x)) {po\x)) \ 
So we assume that S = 5" from now on. Let D be as in Definition 13.11 and let 

Bq := Q G Loo{S'^] ^^))- Then Bq is strictly accretive on the space l-Li from 
(USD and 

—hdivs 
HVs 

Thus by Theorem 17.11 with a = 0, we have bounded functional calculus of BqD in 
BqI-L. Following [U], we have for u G D(V5) that 



B^D 



hLu 




sgn{BoD)BoD 



sgn{BoD) 





HVsu 



so that ||v/iLM||2 ~ ||ifV5M||2 ~ IIV5MII2, using that sgn(_Bo-D) is bounded and 
invertible on BqH and that H is bounded above and below on R(V5). □ 

Remark 8.2. It is interesting to note that we apply Theorem 17.11 with cr = no 
matter what the dimension is. If n > 2, Kato's square root problem on S*" is 
not directly linked to the boundary operator appearing in ffT7|) . associated to the 

equation diVxAVxW = on O^"^", with A - ^ ^ 



H 



i.e. when one can separate 



in the equation radial derivatives from tangential derivatives. This is different from 
the case of the half space (R"" replacing S*") and emphasizes the role of curvature. 

In view of Section HI the second order operator on the boundary associated to this 
div^AVx on O^^", comes from 



(BqD - aNf 



-hL + 

-HVshdivs + o"^ 



with a = {n — l)/2. Thus, the naturally associated Kato square root is \/—hL + cr^, 
and one has 

' ' ~||V5n||2 + ^| 



-hL 



'n~l \2 



U 



Mb- 



9. Natural function spaces 



By Corollary 13. 4^ our method to study and construct solutions u to the divergence 
form equation ([T]) consists in translating this equation to the ODE ( fTTj) for the 
conormal gradient / in R4. x S*". Conormal gradients of variational solutions belong 
to L2(R+ X S*";"!^) as noted in ( ITSi) . On the other hand, the appropriate function 
spaces for / with Dirichlet /Neumann boundary data for u in L2(S'"; C™) are the 
following. 

Definition 9.1. The (truncated) modified non-tangential maximal function of / 
defined on R_|_ x S*", is 

iV,(/)(x) := sup r(^+")/l/x.<i|U,(iy(i,x)), x G S\ 

0<t<co 

where W{t,x) := G R+ x S'" ; \y — x\ < cit, Cq ""^ < s/t < cq} for some 

fixed constants cq > 1, Ci > 0. We assume that cq ~ 1 and ci << 1, so that the 
Whitney regions x) are non-degenerate for t < cq. For a function /o on O^"*"", 



we have iV°(/o) = iV,(/) where f{t,x) := fo{e 
introduction. 



X] 



which properly defines N° in the 
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The (truncated) modified Carleson norm of / in R+ x S*" is 

/ ^ \ 1/2 

||/||c:= sup -r^ esssup|/|2— — , 

\r(Q)<ro IVI JJ{0,r(Q))xQ W{t,x) ^ J 

and the sup is taken over geodesic balls Q C 5"" with volume \Q\, and with radius 
r{Q) less than some fixed constant tq << 1. For a function /o on O^"'"'^, we have 
ll/ollc = ll/llc where f{t,x) := /o(e-*x), which corresponds to ||/o||c as in (g]). 

Note that changing the parameters ci,ci does not affect the results. 

Definition 9.2. (i) For g : 0^+" ^ c^i+n)m^ ^^g^^g norms 

\\g\\l;o:= I |(7(x)|2(l-|x|)rfx, 

ikii^=:=ii^:(^7)ii2+ / i^?(x)rrfx. 

^|x|<e-i 

Let 3^° and X° be the Hilbert/Banach spaces of functions g for which the respective 
norm is finite. 

(ii) For / : R_(_ x S"" — )■ V, define norms 

\\f\\l:= / mlmm{tA)dt, 
Jo 

/oo 
midt. 

Let 3^ and A" be the Hilbert/Banach spaces of sections / for which the respective 
norm is finite. 

The gradient-to-conormal gradient map of Proposition 13.31 is an isomorphism 
y°^y and X" X. 

Lemma 9.3. There are estimates 

1 r^* ~ ds 

sup -/ ||/s||><||iV,(/)||^< / /e4-(R+x5-V). 

0<t<l/2 1^ Jt Jo ^ 

Denoting by y* the dual space of y relative to L2(R+ x S'";V), i.e. the space of 
functions f such that II2 max(t"^, l)dt < 00, we have continuous inclusions of 

Banach spaces 

y* dx d L2(R+ X 5";V) c y. 

Note that Lemma 19.31 shows that another choice of threshold than t = 1 in the 
definition of the norms for X and y would result in equivalent norms. 

Proof. The L2°'^(L2) estimates of ||A^*(/)||2 is an adaption of the corresponding result 
for R+^'\ proved in [U Lem. 5.3]. The remaining statements, except possibly that 
X C L2(R+ X S*"; V), are straightforward consequences. To verify this embedding 
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of X, we use the lower bound on ||A^*(/)||2 to estimate 



POO °o p2-'' p 

/ \\ft\\ldt = J2 / \\ft\\ldt+ / 



WftWidt 

<Y.2~'^'\mf)\\i+ / ii/,ii^rft = 11/11^ 

k=o -^^ 



□ 



The following lemma gives necessary and (different) sufficient conditions for a 
multiplication operator S to map X into y*. Write 

ll^^llcnLoo '■= ll^llc+ ||^||Loo(R+x5")- 

Lemma 9.4. For functions E : R+ x S*" — )■ C'^^"'""-'"*, define the multiplicator norm 
\\E\\* := = sup||j||_^^]^ . Then we have estimates 

||i^||Loo(R+x5") ^ ll^"!!* ^ ll^^llcnLoo- 

Proof. This is an adaption to the unit ball of [H Lem. 5.5]. As in that proof, the 
estimate ||i^||oo ^ ||^^||* follows from the L2°^ estimates in Lemma 1^751 For the second 
estimate we write 



pa 1. POO 



As in in Lem. 5.5], the first term is estimated with Whitney averaging and Car- 
leson's theorem. The second term is controlled with ||^^||oo- In total, this gives the 
bound \\Sf\\y, < mcWfh + ll^llooll/IU as desired. □ 



Remark 9.5. It has been recently proved in [12] that \\£\\* ^ ll^^llcnLoo so all of our 
results use in fact the same condition on S. 

We end this section by introducing an auxiliary subspace ys of y. 

Definition 9.6. For 6 > 0, define the norm 

POO 

11/11^,:= / mlmin{t,l)e''dt. 
Jo 

Let ys be the Hilbert spaces of sections / : R+ x S*" — V such that is finite. 

Clearly ys C y. The motivation for introducing ys is the following result. 

Proposition 9.7. Given coefficients A G Loo(0^"^'^; £(C^^"'""-'''")), which are strictly 
accretive on Tii, there is 6 > such that 



/oo P 
mle'' dt< 



WftWldt, 

/2 



for all f G 4°=(R+;?{) solving dtf + {DB + ^N)f = 0. Hence, if f e y n 
4-(R+; V) and dj + {DB + ^N)f = 0, then f e ys and < \\f\\y. 

The proof of Proposition 19.71 uses reverse Holder inequalities. 
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Theorem 9.8. Fix c > 1. There exist C < oo andp > 2 depending only on n, m, the 
ellipticity constants ||^||oo5 i^a of A and c, such that for any ball B with cB C O^"*"" 
and any weak solution to diVx(v4Vx^^) = m O"^"'"", we have 

i/p / /• \ 1/2 



B J \JcB 

Proof. This result is due to N. Meyers [31] for equations. Here, we make sure that 
the result extends to elliptic systems in the sense of Garding by giving appropri- 
ate references. We begin by noting that the usual Caccioppoli inequality for weak 
solutions 

\ 1/2 / /. \ 1/2 

for any ball B so that cB C O^^", with r its radius, holds for any system that is 
elliptic in the sense of the Garding inequality (Ej). Although not stated like this in 
[T71 Thm. 1.5, p. 46], the proof only uses Garding's inequahty. See also [lO] where it 
is done explicitly for second (and higher order) equations and it is said [TOl p. 315] 
that this applies in extenso to such systems. The constant C depends only on n, m, 
K, ll^lloo and c. Now, this combined with Poincare inequality yields 

1/2 / ^ \ i/g 



B / \JcB 



for '^^^I^P < q < 2. Finally, Gerhing's method for improvement of reverse Holder 
inequalities with increase of radii, presented in [271 Thm. 6.3], applies. □ 



Proof of Proposition \9. 7\ Corollary 13.41 shows that / is the conormal gradient of a 



weak solution to diVxAVx^i = in 0^+". By Holder's inequality and Theorem 19.8 
we have for g = Vx^i the estimate 

^(x)l^lxr^dx < / |^(x)rdx < / |^7(x)prfx 

|x|<e-i / Vj|x|<e-i / VJ|x|<e-i/2 

for < 6 < ("+^)(p-2) ^ This translates to the stated estimate for /, using the 
gradient-to-conormal gradient map from Definition 13. 2[ □ 



10. Semi-groups and radially independent coefficients 

In this section and the subsequent ones, we set a = 

In this section, fix radially independent coefficients Ai and Bq = Ai. We show how 
to obtain weak solutions of diVxAiVxM = inside and outside O^"'"'^ using the semi- 
groups associated to A and A. Later, we show all weak solutions with prescribed 
growth towards the boundary have a representation in terms of these semi-groups. 

Theorem 10.1. Let /q belong to the spectral subspace EqI-L. Then 

ft ■■= e-*Vo 

gives an Ti-valued solution to dtf + D^f = 0, in the strong sense f G C"'^(R+; L2) fl 
C°(R+; D{Do)) and in R"*" x 5*" distribution sense. (In particular f is the conormal 
gradient of a weak solution 0/ diVxAiVxM = in 0^~^".) The function f has L2 
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limit liiiii^o /i = /o o,'>T'd rapid decay ||9j/j2 < Cj,k/t^\\fo\\2, for each k > j > 0. 
Moreover, we have estimates 

If instead /q belongs to the spectral subspace E^Ti, then define ft := e*^/o for 
t < 0. Then dtf + D^f = for t < 0. (In particular f is the conormal gradient of a 
weak solution o/diVxAiVxW = in R" \ 0^+".j Limits and estimates as above hold 
forft,t<0. 

Proof, (i) The rapid decay of ft follows from the lower bound on -DqI-h from Propo- 
sition giving 

= ||A^e-*Vo||2 < ||po)'~^A^e-*Vo||2 ^ ^1 (tA)'e-*Vo||2 < t-'ll/olb- 

(ii) That / is the conormal gradient of a solution follows from Corollary 13.41 and 
it is straightforward to show that the ODE dtf + -Do/ = is satisfied in the strong 
and distribution sense. 

(iii) Next, Hc^t/Hj; < /q°° \\dtft\\2tdt, and the square function estimate WdtftWl'tdt ~ 
||/o||| follows from flHUj) . since dtft = — Ae~*'^/o. This together with the decay from 
(i) with j = 1 shows ||/o||2 ~ \\dtf\\y- 

(iv) It remains to show that II/0II2 ~ II/IIa"- For this, the decay from (i) with 
j = implies it is enough to prove ||iV*/||2 ~ ||/o||2- The proof is an adaptation of 
the results on R^"^" from [SI Prop. 2.56] as follows. 

The estimate ||A^*(/)||2 ^ II/0II2 follows from Lemma 19. 3[ Next consider the 
estimate <. We follow the argument in [3 Prop. 2.56]. By the reverse Holder 
inequalities noted in the proof of Proposition 19.71 applied to a weak solution of the 
divergence form equation with coefficients Ai associated with / = e~*'^°'/o, we can 
bound L2 averages by Lp averages for some p < 2, i.e. A^*/ < N^f in a pointwise 
sense (up to changing to constants Cq, Ci). Since 'ipi^) = e"''^' — (1 + i\)~^ E \I/(S'°^), 
it follows from Lemma 19.31 and Theorem 17.11 or more precisely (IHUj) , that 

\\N^{ij{tDo)fo)h < \\N,mDo)fo)h < II/0II2. 

FoTht := (/+^tDo)^7o wehave \\NP{h)\\2 < ||M(|/o|P)Vf < ||/o||2 by CorollaryEJ 
and the boundedness of M on L2/p. We have proved that ||A^*/||2 ^ II /oil 2- 

(v) The modifications for /o G EqT-L are straightforward, and the correspondence 
with u follows from applying the methods of Proposition 13.31 □ 

Remark 10.2. The assumption a = is used in part (iv) to pass from N^^ to 

with some p < 2. Thus, for any a G R, /o G "H andp < 2, we have || A^f (/)||2 ^ II /oil 2- 
The converse, however, is not clear because p < 2, and this shows that the value of 
a is significant. 

Theorem 10.3. Let Vq G E^L2. Then 

vt := e'^^vo 

gives a solution to dtV+D^v = 0, in the strong sense v G C"'^(R+; L2)nC°(R+; D{Do)) 
and in R"^ x S"" distributional sense. (In particular r'~'^{vt)± extends to a weak 
solution 0/ divx^i VxM = in O^"*"" as in Proposition \3.5[ ) The function v has L2 
limit limt^oVt = vq and rapid decay ||c^ft||2 < Cj,k/t^\\'Vo\\2 for each k > j > 0. 
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(When cr = 0, this estimate for j = only holds for Vq G R{Dq) r\EQL2.) Moreover, 
forp < 2, we have estimates 

\\dtv\\y + \\N:{v)h+\\N,{v,)y<\\voh. 

In dimension n = 1, we have \\v\\x ~ ll'^olh- 

// instead vq G Eq L2, then define Vt ■= c^^Vq for t < 0. Then dtv + Dqv = 
for t < 0. (In particular r~"{vt)^ satisfies diVxAiVx^i = m R" \ 0^+" as in 
Proposition \3. 51 ) Limits and estimates as above hold for Vt, t < 0. 

Proof. The proof, except for the non-tangential maximal estimates, is identical to 
that of Theorem 110.1^ using Proposition 14.41 and Corollary 16.31 When n > 2, the 
estimate of ||A^*(f±)||2 follows, using the same ip as above and reduction to ||A^*((/ + 
itDo)~^^o)±)||2; from Corollary 16.31 and the maximal theorem. When = 1, one 
uses the splitting in Proposition 14.41 we have that e~*^ is the identity on "H^ and 
that A on BqH is similar to A on "H, so \\v\\x ~ ll'^^olh follows from Theorem 110. 1[ 
The modifications when vq G EqL2 are straightforward. □ 

11. The ode in integral form 

Following [3], for radially dependent coefficients we solve fll7p for / by rewriting 
it as 

dtf + [DBo + aN)f = DSf, where St := Bo - Bf 

Recall that solutions ft belong to "H, where "H splits into EqH and E^H by Lemma 
17. 5[ with E^ = x^(-Do) on H. Applying Eq, integrating formally each subequation 
and subtracting the obtained equations we obtain 

(32) ft = e"*^Eo+/o + / e-^'-'^'^E+DEsfsds - / e-^'-'^^'E^ DEJds, 

Jo Jt 

provided limt^o ft = fo and lim^^oo /t = in appropriate sense. We first study 
proper definition, boundedness of the integral operators in ( l32i) on appropriate spaces 
and their limits. The justification of ( 132|) is done in Section [T2l 

Lemma 11.1. If f E 4°"(R+;?/) satisfies dtf + {DB + aN)f = m R+ x 
distributional sense, then 

- f dsvt{t,s)e-^'~'^''E^fJs= [\tit,s)e~^'~'^''E^DSJJs, 
Jo Jo 

/OO POO 
9,r/7(t,s)e-(^-*)^i?o-/.rfs = r]:{t,s)e-^-'-'^''E^D£JJs, 

for all t > 0. The bump functions rjf are constructed as follows. Let rf{t) to 
be the piecewise linear continuous function with support [l,oo), which equals 1 on 
(2,00) and is linear on (1,2). Then let ri^{t) ■.= rf{t/e){l—7f{2et)) andrif{t,s) := 
V%±{t-s)/eMt)r],{s). 

Proof Follow ^ Prop. 4.4]. □ 
Define for / G 4°"(R+; L2{S''; V)), 

pt POO 

Slft:= / vtit,s)e-^'-'^''E+D£sfsds- / r]7 (t, s)e-^'-'^''E^ DSJ^ds. 
Jo Jt 
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In fact, this formula makes sense by extension thanks to the following algebraic 
relations. 

Lemma 11.2. We have S^ft = S\ft — crS^ft = DS'xft, where 
Jo Jt 

pt poo 
JO Jt 

pt POO 

JO Jt 
Here := E^B^^P^^, E^ := E^NB^^P^^^, with P^^^ as m Proposition^. 

Proof. Here, Bq^ denotes the inverse of the isomorphism Bq : l-L ^ BqH. Since 
H{D) = we have 

E^D = E^DPh^ = E^{{DBo + aN) - aN)B^'Pj,^^ = D,E^ - aE^, 

Using that e~"^ and e~"^A extend to bounded operators on "H, this also shows that 
g-MA^+^ extend to bounded operators on L2 for u > 0. We now readily obtain 
S'X = S'X — crS'X. The identity 5*^^ = DS% is a consequence of the intertwining relation 

b{Do)D = Db{Do) 

between the two functional calculi. □ 
Theorem 11.3. Assume < 00. We have bounded operators 

with norms < uniformly for e > 0. In the space X there is a limit operator 

G C{X; X) such that 

lim \\S\f - SAf\\L2{aML2) = 0, for any f e X,0 < a <b < 00. 

The same bounds and limits hold for S\ and S"^ on X . 
In the space y, there is a limit operator G jC{y; y) such that 

lim\\S^J-Syf\\y = 0, for any fey. 

The same bounds and limits hold for S\ and S\ on y. 

Let Sa '■= linie-s-o 5"^, Sa '■= linie_>.o 5"^ and Sa '■= linie^o denote the limit 
operators on y from Theorem I1L3I Since X is densely embedded in 3^, these limit 
operators restricts to the corresponding limit operators on X from Theorem I1L3[ 

One sees that Sa = Sa — cfSa holds, and that 

SAft = lim ^ e-^'-'^^'E+DSsfsds - ^ e-^'-'^^'E^ DSsfsds 
with convergence in L2{a, b; L2) for any < a < 6 < 00, both on y and X . 
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Proof. The proof is essentially an application of [H Sec. 6], where the results were 
proved abstractly. Given Theorems 17.11 and IIU.I^ these results from ^ apply. In 
particular, this makes use of the holomorphic operational calculus of Dq, where 
more general operator- valued holomorphic functions are applied to Dq. It is straight- 
forward, given Theorem 17.11 . to adapt the results in [H Sec. 6-7] and construct this 
S^^ operational calculus of Dq, and we omit the details. 

(i) Consider the operators S\ : X ^ X. Here [H Thm. 6.5] shows that 5*^ : 
L2(R+,dt; L2) — !■ L2(R+, (it; L2) are uniformly bounded, with norm < ||i^||oo5 and 
converge strongly in C{L2(R+,dt; L2)) as e — > 0. Moreover, [H Thm. 6.8] applies 
and shows that 



S:,ft = Z^{£f)t + Ve{t)e 



-tA 



where : L2(R+, dt/t; L2) — )■ L2(R+, dt/t\ L2) are uniformly bounded and converge 
strongly as e — )■ 0. These estimates build on the square function estimates and make 
use of the operational calculus for Dq. On the other hand, using Theorem IIU. II and 
Theorem 17. Ij the last term has estimates 



r],{s)Ae-'''E^Ssfsds 



< 



X 



7],{s)Ae~'''E+Ssfsds 



sup 

||/l||2 = l 



i.sA*e-^^'h,Us)E^£sfs)- 



< WVeSfh^ < \\S\ 



X, 



and is seen to converge strongly in C{X, ^2(0., b; L2)) for any < a < 6 < cx), as in 
in lem. 6.9]. Piecing these estimates together, we obtain 



\SAfh< 



\Z'{£f)\\ 



L-zidt/t-M) 



-\\S\f-Z\£f)\U 

\£\\oo\\f\\L2{dt;L2) + ll^l 

with strong convergence in C{X , ^2(0, b; L2)). 

(ii) For the operators S*^ : A* — )■ A", we note that the estimates for 5^ go through 
when replacing E^ by E^. Since S\ = A"^^^ (with Eq replaced by Eq) and 
: L2{dt/t;'H) — )■ L2{dt/t;'H) is bounded, it only remains to estimate the term 
rj^{t)e~*'^ rj^{s)e~^^EQ Esfsds. But again using the boundedness of gives 



-tA 



r],[s)e -'^E^Ssfsds 



< 



X 



ri^e-'^'E^SJsds 

) 

A / r],{s)e-'^E^SJsds 



< 



2 

sA ; 



and the rest of the estimates go though as for 5*^. Altogether, this proves the stated 
bounds and convergence for S\: X ^ X. 

(iii) Next consider the operators S\:y^y. We have 

ll^y lb < ll%X*<l/)b + II^A(X*>l/)b 

L2{tdt;L2) ~r 2{dt;L2) 
^ \\^\\*\\Xt<lf\\L2{tdt;L2) + \\^\\oo\\Xt>lf\\L2(dt;L2) ^ 11^11*11/11^5 
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where the L2{tdt] L2) estimate follows from ^ Prop. 7.1] and the L2{dt] L2) estimate 
from in Prop. 6.5], along with convergence. This immediately gives the estimates for 
S\:y^y since A"^ : L2{tdt]'H) L2{tdt;'H) and : L2{dt]'H) L2{dt]'H) 
are bounded. □ 

Denote by C(a, h] L2) the space of continuous functions (a, h) 3 t ^ Vt & -^2(5'"; V). 

Theorem 11.4. Assume < 00. If n > 2, then S\f G C(0, oo;L2) for any 

f & There are bounds ||5'^/t||2 ^ 11^11*11/11^; uniformly for all f G 3^, t, e > 0, and 
for each f Ey there is a limit function SAf G C(0, 00; L2) such that lime^^o ll'S'A/t ~ 
SAft\\2 = locally uniformly for t > 0. We have the expression 

pt poo 

(33) SAft = / e-^'-'^^E^£JJs - / e-^-'-'^'^^E^SJ^ds, 

Jo Jt 

where the integrals are weakly convergent in L2 for all f & y and t > 0. Finally, 

SAf = DSAf holds in R_|_ x S"' distributional sense for each f & y. 

If n = 1, then the above results hold ify is replaced by ys, for any fixed 6 > 0. 

Proof, (i) Consider first the case n > 2. The proof is a adaption of the proof of [H 
Prop. 7.2], which we refer to for further details. We split the (0, t)-integral 

/\+(t, s)e-(*-^)^(/ - e-^'^)E^SJ.4s + e"*^ /\+(t, s)e-'^E+SJsds, 
Jo Jo 

The same duality estimate of the second term as in [H Prop. 7.2], given Theo- 
rem 110.11 and Lemma 14.21 goes through here. For the first term, we note the es- 
timates ||e-(*-'')^(/ - e-2''^)|| < min(s/t, 1, l/(t - s)). For t < 2, this yields the 
bound ll^^lloo jlis/t) \\fs\\2ds < ll^'llooH/Hy- On the other hand, for t > 2 we have the 
estimate 

mioo 1^1^' '-Wfshds + II ' -^Jf^W.ds + I' Jfshds^ <|l^||oo||/b. 
The (t, cx))-integral is estimated similarly, by splitting it 

POO POO 

/ C(t, s)e-(^"*)^(J - e-2*^)Eo"^J.c?s + e-*^ / ^(t, s)e-'^E^£J,ds, 
Jt Jt 

The second term is estimated as before, and for the first term we note the estimates 
||g-(s-t)A^j _ g-2tA-j|| < j-Qin^^/g^ i^(^g _ f^^^ which give the bound 

ll^lloo ^-Wfshds + Qj^^Wfshds^ < ll^llooll/b. 

(ii) Consider next the case n = 1. Since e~^^ = J on "H^ and E^ = on "H-*-, 
we also need to estimate the L2-norm of 

uniformly for t > 0, where Pg^ is projection onto "H-*" from Proposition 14.41 So it is 
enough to obtain the bound 



00 

" 



-^Bo^sfs\ds 



< ll^'IU 
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On one hand, we obtain from Proposition 19.71 the estimate 



< ll^l 



oo 

oo ; ll/slhc^-S 
1 



< ll<£^l 



1/2 



OollJ 113^5- 



On the other hand, note that A, hence Bq^, is pointwise strictly accretive by 
Lemma 15.11 and by the exphcit expression in Lemma 15.51 (expressed in other co- 
ordinates), P^^ maps into constant functions and |-P^gM| < f^i \u{x)\dx. Thus 



^ / \Ss{x)\\fsix)\dxds. 

2 JO Js^ 



< 



Pick h : R+ x ^ V such that = 1 and \Ss{x)hs{x)\ = \Ss{x)\ when 

s < 1, and hs{x) = when s > 1. Cauchy-Schwarz inequahty yields 

' / \S,{x)U{x)\dsdx < \\Sh\\y4f\\y < \\SUhU\f\\y < \\SUf\\y. 

Js^ 

This completes the proof of the estimate of || 5*^/411 2- 

(iii) As in the proof of [U Prop. 7.2], replacing rjf by rjf — t]^ in the estimates 
shows convergence of S"^ and yield the expression for the limit operator. The relation 
Sa = DSa follows at the limit from the relation in Lemma [1L2[ □ 

We turn to boundary behavior of the integral operators at t = 0. 

Lemma 11.5. Assume \\£\\* < oo. 

(i) Let f E X (or f Ey) and define /° := SaJ ■ Then and f satisfy 

{dt + D,)f = D8f 

in R+ X distributional sense. If f E X , then there are limits 



limt"^ 



2t 



\SAfs-h \\lds = 0, 



where h := — e '^^Eq DEgfsds E Eq1-L has hounds \\h ||2 < 
(ii) Let n > 2. If f E y and v := SAf, then 

{dt + Do)v = Sf 

in R+ X S*" distributional sense, and there are limits 

lim \\SAft -h^\\2 = 0, 

where h~ := — c'^^^Eq E^fgds E Eq L2 has bounds \\h~\\2 < 
results for SAf hold when replacing y by ys, for any fixed 6 > 



Ifn = 1, these 



Proof, (i) By the convergence properties of 5*^ from Theorem 1 11. 3 [ it suffices to show 
that for (j) E C^(R+ x 5"; C*^^"*""^"^) there is convergence 



J ({~d4t + BlD + aN)(t>uft)dt-^ J {D<Ps^Ssfs)ds, 



e 0, 
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where := 5^/*. For the term (0, t)-integral, Fubini's theorem and integration by 
parts gives 

/ / vtit,s){{-dt + A*)(l)ue-^'-'^''E+D£sfs)dsdt 
Jo Jo 

poo / POO \ 

= -J^ U vtit,s)DiE^ydtie-^'-'^'''(l)t)dt,£sfs]ds 

f*00 / /'OO \ POO 

= J n idtvt)it,s)DiE^ye-^'-'^''*Mt,£sfs]ds^ J {D{E^r(l>s,£sfs)ds. 

Adding the corresponding hmit for the (t, oo)-integral, we obtain in total the hmit 
j^{Dct)s.£J,)ds, since D{{E^y + {E^)*) = {{E^ + E^)DY = D* = D. 

To prove the hmit of SAft for / G A", we note from the proof of Theorem 111.31 
that 



SAft = ZAft + e-'^' / e-'''E,DSJJs, 
Jo 

where ZAf G y*- When taking hmits e — 0, we have used [H Thm. 6.8, Lem. 6.9]. 
This proves the stated hmit. 

(ii) To prove {dt + Dq)v = Sf, we let t G (a, b) and differentiate 5*^/ to get 

dtS^t = \ j'^ e-'\E^£t^sft^s + E^£t+sft+s)ds - 5o(51/i), 

for small e. The first term on the right is seen to converge to £f in L2(a, 6; L2) as 
e — 0, with an argument as in |31 Thm. 8.2]. Note that this uses Eq +Eq = I, which 
holds also when n = 1 by Definition 17. 4[ Letting e — ?■ 0, we obtain dfV = Sf — Dqv 
in distributional sense, since (a, h) was arbitrary. 

The limit for SAft when f & y (or when n = 1) is proved as in [U Prop. 7.2, 
Lem. 6.9]. Note in particular this uses an identity 

SAft = ZAft + e-*^ / e^^'^^E^SJ^ds, 
Jo 

with ZAf G C(0, 00; L2) and lim^^^o ZAft = in L2. □ 

12. Representation and traces of solutions 

We now come to the heart of the matter. The natural classes of solutions for the 
Dirichlet and Neumann problems, with L2 boundary data, use the spaces y° ^ y 
and X° ^ X from Definition 19. 2[ 



Definition 12.1. (i) By a y"-solution to the divergence form equation, with coeffi- 
cients A, we mean a weak solution u of diVxAVw = in O^"^" with HVx^Hyo < 00. 

(ii) By an X°-solution to the divergence form equation, with coefficients A, we 
mean the gradient g := Vx^ of a weak solution u of diVxAVu = in O^'*'" with 
\\g\\xo < 00. 

Note the slight abuse of notation when referring to the gradient VxM rather than 
u as an A:'°-solution. The reason for this convention, here as well as in [1], is that 
the Neumann and regularity problems are BVPs for g (and not for the potential u), 
and ^'''-solutions is the natural class of solutions for these problems. This point of 
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view is the one that lead us to our representations. However, when more convenient 
we call the potential u itself an Af^-solution. 

Remark 12.2. (i) No boundary trace is assumed in our definitions, but will be 
deduced. 

(ii) The semi-norm || Vx^||yo on 3^°-solutions is modulo constants, which is unusual 
for Dirichlet problems. Once we have shown that 3^°-solutions have boundary traces, 
we will be able to put constants back in the norm in a natural way. 

(iii) For any A'^-solution g, the potential u has a boundary trace in appropriate 
sense (replacing pointwise values by averages) and the trace belongs to 1^2^(5'"; C™). 
This is essentially in [30]. We also recover this from our representations. See Sec- 
tion [131 

Here and subsequently, we use the notation e~*^g to denote the function (t, x) H- 
{e~^^g){x). Similarly for e~^^g. 

12.1. A'°-solutions. We begin with representation and boundary trace for solutions 
of the corresponding ODE. 

Theorem 12.3. Assume that \\S\\^ < oo. Let f E X. Then f E L^°'^(R+;'H) 
satisfies dtf + {DB + ^^^N)f = in R+ x S*" distributional sense if and only if f 
satisfies the equation 

(34) ft = e-'^h+ + SAft, for some h+ E E^H. 
In this case, f has limit 

(35) ||/,-/o||^rfs = 0, 

where /o := h"^ + and := — c^^^Eq DSsfgds E Eq'H, with estimates 

max{\\h+h,\\h~h)^\\foh<\\fh. 
If furthermore I — Sa is invertible on X , then 

(36) / = (/ - 5^)"^e-*^/i+ 
^nd ||/|U< ll/i+lh. 

Proof. The proof is an adaption of [H Thm. 8.2] to which we refer for details. Here 
is a quick summary. 

We show that / satisfies (|T71) if and only if / satisfies (1341) . Assume (1171) and 
apply Lemma ril.li Letting e — )■ and applying Theorem 111.31 we obtain the stated 
equation for /, with as a certain weak limit as in part (i) of the proof of [H, 
Thm. 8.2], with A = |Do| here. 

Conversely, if f E X satisfies (IMl) . then we apply Lemma [11.51 with f° '■= f — 
e-*^/i+. Since {dt + Do)e~^^h+ = and e-*^/i+ E X hj Theorem [JOlIl it follows 
that / satisfies ffTTj) . 

Lemma [11.51 also shows existence of the limit /q. The stated estimates follow as 
in part (iii) of the proof of jH, Thm. 8.2]. 

If / — S*^ is invertible, the equation ( !36l) follows immediatly from ([34l) . and the 
estimate \\f\\x ^ II ^^11 2 follows again from Theorem llO.li □ 
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Theorem 12.4. Assume that < oo. Then g is an X°-solution to the divergence 
form equation with coefficients A if and only if the corresponding conormal gradient 
f E X satisfies the equation 

(37) ft = e-'^h+ + SAft, for some h+ G E+U. 

In this case, g has limit 

\im-^— I |5((x) - 5(i(x)prfx = 0, 

r^ll-r A<|x|<(l+r)/2 

where gi := {Bofo)j_n + (/o),, and \\gi\\2 < WgWx" holds. If furthermore I - Sa is 
invertible on X , then ||^'''||2 ~ llfi'ilh ~ 

Proof. The equivalence follows from Corollary 13.41 and Theorem 112.31 The limit and 
the estimates follow on applying the conormal gradient-to-gradient map of Proposi- 
tion [213] from the ones satisfied by /. □ 

It is worth specifying the previous theorem in the case of radially independent 
coefficients. 

Corollary 12.5. Assume A is radially independent. Then any X°-solution has 
corresponding conormal gradient given by f = e~*^/i^ for a unique h'^ G E^T-L. 

Remark 12.6. A careful examination of the proof of Theorem 112.31 in the case of 
radially independent coefficients, shows in fact that for / G L2°'^(R_|_; "H) the weaker 
condition supQ<i<^/2 i ^t* WfsWi^s < oo is sufficient to obtain this corollary, as in this 
case S'a = 0. 

12.2. 3^°-solutions. We now turn to representations and boundary behavior per- 
taining to 3^°-solutions. 



Theorem 12.7. Assume that < oo and f Ey . 

(i) Then f G L^°=(R+; U) satisfies dtf + {DB + ^N)f = R+ x 5" distri- 
butional sense if and only if f satisfies the equation 

(38) ft = De'''^h+ + SAfu for some h+ G E+L2. 

Here h'^ is unique modulo Eq'H'^ and \\h'^\\i^/y^± < \\f\\y, and if furthermore I — Sa 
is invertible on y then 

(39) / = (/ - SA)-'De-H^ 
^^th \\f\\y<W\\ 

(11) If ^ holds, letvt := e-'~^h+ + SAff Then f = Dv anddtV+{BD-^N)v = 
0, and Vt has L2 limit 

(40) \im\\vt - vo\\2 = 0, 



where Vq := h'^+h and h := — e ^^Eq Sgfsds G Eq L2, with estimates \\h II2 ^ 

ll/b (^nd 

(41) ll^t||2<||/^+||2 + ||/b, forallt>0. 
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Proof. The proof is an adaption, with some modifications, of ||4i Thm. 9.2] to which 
we refer for omitted details. 

(i) Assume f lT7|) . We apply Lemma 111.11 to /. Letting e — and applying 
Theorem lll.3[ we obtain for / the equation 

ft = ft + SaJu 

with the limit ft := \im^_^Q J^'^ e~^^~^^^EQ fgds. From here, one can proceed as 
in in Thm. 9.2] to represent ft as DQe~^^h^ for some G EqTI, or use a simpler 
argument (owing to the boundedness of the boundary here): since Dq : EqU — )■ 
EqU is surjective, there exists ht G EqU such that ft = D^ht. From there and 
ftQ+t = e~*^/to, we conclude as in |1] that the weak L2-limit /i"*" := \\mt^oht exists 
and that ft = Doe-*^h+. 

To write Doe~*^/i+ as De^^^h^ for some G EqL2, we use Lemma [4. 61 Indeed, 
there is an isomorphism M : "H — )■ L2/'H'^ with Dq = D o M on D(Z)o). It is 
easy to see that the restriction of M to EqT-L maps onto Eq L2/ EqI-L-^ . Now, on 
□ (Do), Dqc-^^ = c-^^Dq = e-*^D o M = De"*^ o M. By density and boundedness, 
the left and right terms agree on "H. Thus, h'^ = Mh^ G Eq L2/ EqI-L-^ satisfies 
Doe-'^h+ = De-'~^h+. 

We conclude that ft = De~^^h~^ + 5*^/*, with estimates 

(42) fh^h^/H- ^ \\h^h ^ \\Doe~'^h^\\y = \\f - < \\f\\y. 

The middle equivalence uses Theorem 110. 1[ 

(i') Conversely, if / G 3^ satisfies ( !38l) for some G Eq L2, then we apply 
Lemma [mS] with f° = f - De^^^h+ = f - Doe~^^h+ with h+ G E^^n given by the 
isomorphism above. Since {dt + -Do)-Doe~*^/i^ = 0, it follows that / satisfies (fTTT) . 
For the estimate of \\f\\y when J — S'a is invertible on 3^, use that the last estimate 
in 0421) in this case is ~. 

(ii) Lemma [11.51 and Theorem 1 1 1 . 41 show the ODE satisfied by v, existence of the 
limit vq and the estimates of 11^^(112 and ||/?.~||2- This completes the proof. □ 



Corollary 12.8. Assume that \\S\\^ < 00. With the notation from Theorem \1 2. 7 
the following holds. 

(i) Any y° -solution u to the divergence form equation has representation Ur = 



n-l 



r 2 (vt)± with r = e , for some v as in Theorem 12. 7 , boundary trace in the sense 
lim^^i \\Ur — M1II2 = 0, and there are estimates 



— 2— i 

Ft 2 < r 2 VxM 



ui{x)dx 



r G (0,1). 



(ii) The map taking y° -solutions u to boundary functions = EqVq G EqL2 is 
well-defined and bounded in the sense that 



'5" 



Ui{x)dx 



(Hi) If furthermore I — Sa is invertible on y, then this map is an isomorphism 
and its inverse EqL2 3 h'^ ^ u & {y° -solutions} is given by 

(43) Ur := r-"^ ((/ + Sa{I - SaY^ D)e-'~^h+ 
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with estimates \\V:x.u\\yo + | ~ ||/i''"||2- 

Proof, (i) Let / be the conormal gradient of u and define and v applying Theo- 
rem 112.71 As in the proof of Proposition 13.5^ it follows that 

for some c G C™, where r = e~* G (0, 1) and cr = 

Recall that by f lHSj) . /i+ is uniquely defined in EqL2 modulo E^l-L-^ and we now 
use this freedom to choose it in EqL2 such that c = 0. Indeed, by Lemma I7.5[ 
E:^U^ = N-n^ = {[c 0]* ; c G C"'} and since A = al on U^, we have 

(44) e-*^([c 0]*)=e-"*[c O]*, c G C™. 

(The superscript t of the brackets denotes transpose.) Replacing h'^ by li^ — \c O]*, 

then ft remains unchanged, e~^^h~^ is replaced by e~*^h'^ — e~"^ \c O]*, and {vt)i_ 
by {yt)i_ — e~'^*c. Thus we may assume c = 0. 

As Vt has an L2(S'"; C^^"*""^™) limit t^o when t — )■ 0, one can set Ui := (fo)± and 
converges in ^2(5""; C") to Ui. For the estimate on ||u,r||2 it suffices to prove 



n-l 



\Wr-m\\^<r 2 llVx^ilb", rG (0,1). 

with m the mean value of ui on 5". We may assume that m = as by (jHj) this 
amounts to modifying modulo N~'H'^ without changing the conormal gradient 
/ of M. We have 

lk.||2<r--||^^*||2<r-'^(||/^+||2 + ||/b). 
By orthogonal projection onto N^Ti-^, it follows ||/i'''||2 ~ ||^''^||l2/w-l + I fsn{h^)±dx\ 
since G EqL2. We can now conclude since ^ \\f\\y ^^^y since m = 0, 



gn 



< \\h~\\, < 



2 ^ \\J \\y 



{ui — {h )j){x)dx 

IS" 

(ii) The argument using (jHj) shows that given a 3^°-solution u and its conormal 
gradient /, there exists G EqL2 such that = r~°'(e~*^/i"'" + SAft)±- Moreover, 
= EqVq by construction and the estimate ||/i~^||2 ^ llVx^Hyo + \Jgn Uiix)d: 



x\ 



follows from the above argument. To define the map and prove its boundedness, 
it suffices to show uniqueness of such G EqL2. So assume Ur = {e~^^h~^ + 
SAft)± = r~"{e~^^hi + SAft)± with / the conormal gradient of u and h'^,hi G 
E^L2. Thisjmplies that ft = De'^'^h^ ^^^ft = De'^^hf + SAft so we know that 
h+ -h+ e E+n^ by Theorem (HITl As E+U^ = N-R^, write h+ - h+ = [c O]* 
with c G C"*. We have from (j44j), = r~''{e-^^{h+ - hf))j_ = c. 

(ni) Given h+ G E+L2, define ft := (/ - SA)-^De-'^h+, Vt := e-'^h+ + SAft 
and Mr := r~'^{vt)±. By Theorem 110.31 and Lemma 111.5^ v satisfies the equation 
dtv + Dqv = 0, and by Proposition 13.51 u extends to a 3^°-solution and / is the 
conormal gradient of u. For the continuity estimate HVx^iHyo + \JgnUi{x)dx\ < 
\\h+\\2, Theorem (HITj implies \\f\\y < \\h+\\2 and \Jg„Ui{x)dx\ < \\ui\\2 < H^^olh < 
II II 2 + ||/||y ^ II ^^11 2- This map clearly inverts the map in (ii). This completes 
the proof. □ 
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It is worth specifying the Corollary 112.81 in the case of radially independent coef- 
ficients. 

Corollary 12.9. Assume A is radially independent. Then any y°-solution is given 
hyu = r ~ {e~^^}i^)^ for a unique G EqL2 with ||/i~^||2 ~ II Vx^^Hyo + j/^n Uidx\ . 

12.3. Conclusion. It is clear from ( l36l) that provided J — S'a is invertible on X, the 

ansatz 

E+V, X : h+ ^ ft = {I - SA)~^e-'^h+ 
maps onto all conormal gradients of A'^-solutions to the divergence form equation 
with coefficients A. 

Similarly, (H3l) implies that provided / — 5*71 is invertible on y, the ansatz 

^y-.h+^Ur := (^(/ + Sa{I - SAy^D)e~H+^ , 

maps onto all J^^-solutions to the divergence form equation with coefficients A. 

Thus we have a way of constructing solutions and our two main goals towards 
well-posedness results are the following. 

First understand when invertibility of I—Sa holds. This will be done in Section [T6l 

Secondly, introduce the boundary maps that connect the traces of solutions to the 
data for the BVPs and invert them. This is the object of Section [T71 

Before we do this, we continue with different a priori representations of solutions 
in the next section. This will be useful to prove non-tangential maximal estimates 
and obtain convergence of Fatou type at the boundary. 

13. Conjugate systems 

The results in the preceding section allow to represent A'°-solutions in terms of 
the conormal gradient /. Actually, if one is interested in u itself, one can try to 
further describe the corresponding potential vector v. Similarly, representation of 
3^°-solutions is embedded into a potential vector v but it could be interesting to 
describe the properties of the conormal gradient /. Both are related by the rule 
Dv = f. This leads us to the following notion. 

Definition 13.1. A pair of conjugate systems to the divergence equation with co- 
efficients A is a pair {v, f) e 4°'=(R+; ^2(5'^; V)) x 4°'(R+; MS""; V)) with 

(i) Vt G D(D) for almost every t and WDvtW^dt < oo, 

(ii) V is an R"*" x S'"-distributional solution of (HM . 

(iii) ft = Dvt for almost every t > 0, 

(iv) / is a valued R"*" x S'"-distributional solution of (fTTl) . 

By Proposition 13.51 and its proof, a pair of conjugate systems is completely deter- 
mined by V satisfying (i-ii). That is, / defined by (in) automatically satisfies (iv). 
Moreover, the function 

(45) u,:=r-^^-^^/\vt),, r = e-*G(0,l), 

extends to a weak solution of diVxAVxW = in O^"^"' and / must be the conormal 
gradient of u. We say that a weak solution u and a pair of conjugate systems (f , /) 
to the divergence form equation for which ( H5i) holds are associated. 

It is our goal to give a description of the pair (not only / or v) in each case. Recall 
that in integrating Dv = /, Vt is only determined by ft modulo Ji^ so there is a 
choice to make. 
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Theorem 13.2. Assume < oo. Letu be an X°- or y° -solution. Thenu has an 
L2{S"'] C"^) trace ui at the boundary and there exists an associated pair of conjugate 
systems given by 



(46) 



e ^^vq + Wt 



with the following properties. 

(i) If u IS an X°-solution, then ui G {S"" ; C"") , {vq, fo) e D{D) x n with 
Dvq = fo, IIV5M1II2 < II/0II2 < IIVx^IIa-", ll^^olb < ||Vxm|U° + \ jg„uidx\, 
Dwt = wte y*, vt G C(R+; L2) and \\vt - vo\\2 + Ndb = 0{t) for t > 0. 

(ii) Ifms a y -solution, then G ^2(5"; C™) , {vq, fo) e x W^^S""; V) with 
Dvq = fo, \\ui\\2 < |bo||2 + ||/o|lvy-i < ||VxMbo + |/5„Mi(ix|, Dwt = WtEy, 
Vt G CilV'] L2) and \\vt — Vq\\2 + \\wt\\2 = 0(1) for t > and o(l) for t — )■ 0. 

Besides and 3^°-solutions to the divergence form equation, we shall in the 
following sections also consider the following classical class of variational solutions. 

Definition 13.3. By a variational solution to the divergence form equation, with 
coefficients A, we mean a weak solution of diVxAVu = in O^"'"'^ with || Vx^^lb < 00. 

It is illuminating to see how the representation for variational solutions lies in 
between the ones for X°- and 3^°-solutions, independently of solvability issues which 
are well-known for variational solutions. We state this result without proof as it 
is not used in this paper. Note that, as compared to Theorem 113. 2[ the Carleson 
condition \\S\\^ < 00 is not needed in the following result. 

Proposition 13.4. Let u be a variational solution to the divergence form equa- 
tion with coefficients A. Then u has an L2(5'";C™) trace ui at the boundary and 
there exists an associated pair of conjugate systems given by fHUj) with the following 
properties: 

u, G Wl'\S^^-C^), (voJo) e D(|D|V2) X W-'^'iSW) with Dvo = fo, Hh < 
||Vxm||2+ \JgnUidx\, \\ui\\^i/2 < ||/o||,^-i/2 < llVx^lla, Dwt = Wt E L2(R+;L2), 
Vt G C(R+; L2) and \\vt - Vo\\2 + Wmh = 0{t'/^) for t > 0. 

Here W2 is equipped with homogeneous norm and W2 is its dual. 

Proof of Theorem \13.2[ (i) From Theorem I12.3[ we have 

f, = e-'^h+ + SAft = e-'^'fo + Wt, Wt := SAft - e^^^/i". 

with fo = h+ + h- e H, II/0II2 < ||Vxn|Uo and = - j^^ e-'^'Eo DSsfsds. ^ 

We define vo, h~^, h~ and v as follows: is the unique element in Eq L2/ Eq'H'^ 
such that Dh+ = h+{= Do{Do^h+)), := - e-''^Eo Ssfsds, vo = + and 

:= e-H+ + SAft = e-'Ko + Wt := SAft - e-^- . 

Clearly, h+ G D{D). Next, Ah~ G L2 because £f G 3^*, so h~ G D{D) = D(A) 
and Dh~ = h~ . So vo G D(D) and Dvo = fo- 

The estimate on ||e~*^?;o — "yolb follows from Vo G 0{D). 

Next, Dwt = Wt by construction and Wt G y* from the proof of Lemma [11.51 (In 
fact, Wt is nothing but Z^/i defined in that proof.) 
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The estimate on lltytlU follows from 



-{s-t)A 



using SfEy*, the uniform boundedness of the semigroup and its decay at infinity. 
Details are left to the reader. 

Eventually, as in Corollary 112.81 one can adjust by adding an element in 
N~T-i-^ such that u and v satisfy fH3|) . In particular, u has an L2 trace. It also 
follows that / is the conormal gradient of u with a limit /o when t — )■ by ( !35l) . So 
ui e W^{S^- C™) with IIV5M1II2 < II/0II2. 

(ii) By Corollary 112.81 we have description of 



Vt 



-tA 



Vq + Wt, Wt = SAft - e 



with Vq = + h such that u and v satisfy ( H5l) and of trace and growth estimates 
for ||e~*^fo — "Wolh + ll'W^tlh- It remains to consider the representation of /. We have 
by Theorem 112.71 

ft = De^H^ + SAft = De~'\ + Wt, Wt = SAft - De^H" = Dwf 
Define /o := Dvq in distribution sense, so that /o G W2^{S'^]V) and H/oHvi/^ 
IK'olh- We obtain 



1 < 



ft = e-'^f, + wt 

and here, the action of e~*^ is extended to 1^2"^ (S*"; V) by extending the interwining 



formula De 



-tA ^ (,-tAjj 

14. NON-TANGENTIAL MAXIMAL ESTIMATES 



□ 



Theorem 14.1. Assume \\£\\cr\L^ < 00. Then any y° -solution to the divergence 
form equation with coefficients A satisfies 

2 



\Mi<\mu)\\i< 



|Vx^^P(l - |x|)dx 



Ui{x)dx 



When n = 1, the conjugate u of a y°-solution u also satisfies the estimates 



\ni\\i<\\N:mi< 



|Vx^^|^(l - |x|)(ix 



ui{x)dx 



The proof follows the strategy of [1] with a slight modification in view of preparing 
the proof of almost everywhere non-tangential convergence. 

Proof. The estimate ||iV^(ti)||2 > [[o follows from Lemma l9.3l and Corollary ll2.8( i). 
For the upper bound, we proceed as follows. From the representation Ur = r~'^{vt)± 
with Vt = e~*^fo + Wt in Theorem 113.21 it is enough to bound ||A''*((e~*^t;o)x)||2 and 



A^*(ii'x) lb- Theorem 110. 3[ and Lemma [14.21 below, show 

\y 



(47) \\N:iu)h < \\voh + \\f\\y < fh^h + ll^^^lh + 



< \\h' 



htdx 
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and ||/i+||L2/w-L ^ \\f\\ y, \I sn h+dx\ = \ /^^(mi - h^)dx\ < \ J^^uidx] + \\f\\y as in 
the proof of Corollary 112.81 

When n = 1, replacing A by the conjugate coefficients A defined in Section[5]in the 
above argument, and using |Vxw| ~ |Vxm|, proves the estimates of ||A^°('u)||2. □ 

Lemma 14.2. Assume \\S\\cr\L^ < oo. Then we have for each p < 2, 

\\N!{w)h + \mw^)h<\\S\\cnLj\f\\y. 

Here is defined similarly to N^^, replacing L2 averages by Lp averages. When 
n = 1, we also have 

\mwii)h<\mcnLj\f\\y. 

Furthermore, these estimates hold with w replaced by the truncation XktW, and 
11/11^ replaced by \\ft\\lmm{t,T)dt, for any r < 1. 

Proof. The proof will follow closely the strategy of |H Lem. 10.2] on R^"^". We 
remark that < A^* pointwise. Thus we will work with N^, and indicate when 
we need to consider or the normal component. Recall that A^* estimates the 
truncation of the function to J: < 1. 

(i) From Wt = SaIi — e~^^h~ and the definition of h~, 



pt _ roc _ ^ poo _ ^ 

Jo Jt Jo 

rt _ poo 



+ e-*^ / e-'^£sfsds = 1-11 + III. 
Jo 

Note that Eq + Eq = I (also in dimension n = 1) is used in getting III. For the 
first two terms, we use Schur estimates as follows. Since ||e~^*~^^^(/ — e~^^^)|| < s/t, 
we have as in [U Lem. 10.2] 

immi < (^f^ st~'m,ds^ dt/t < iix*<i/ii^. 

Similarly, as ||e~*^*~*^'^(/ — e~^*^)|| < t/s, we have 

\mii)\\i < ts~'\\f,hds^ dt/t 



< 



00 / /■min(s,l) 

2 r„ _ II p||2 



\J0 



tdt/t]\\fs\\ids=\\fry. 



Note that the estimates so far hold for all w, not only for its normal component. 
By inspection, the stated estimates of the truncated maximal function hold for these 
terms. 
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(ii) It remains to consider III = e *^ e ^^Sgfsds. To make use of off-diagonal 



estimates in Lemma 16.21 we need to replace e *^ by the resolvents (/ 
this end, define ijjt{z) := e'*'^' — (1 + itz)^^ and split the integral 



itD^) 



To 



-sA 



Jo Jt 

+ / {I + itDor\e-''' - I)Ssfsds + [I + itDo)-' [ S.Jsds. 
Jo Jo 



For the first term, square function estimates show that iptiDa 
continuous, and Theorem 111.41 shows || e~'^^£sfsds\\2 < 



y 



L2^y* d X IS 

(or < ll/b, when 



n 



1, but 



< 



y for conormal gradients of solutions by Proposition 19.71) . 
For the second and third terms, we proceed as above for I and II by Schur estimates 



using \\ilJt{Do)e' 



-sAu < 



t/s, and II (/ + itDo) 



/)ll <='^A- 



(iii) It remains to estimate (/ + it-Do) ^ /J £sfsds, and this is where we use ||i^||c- 
Consider first A^^. Fix a Whitney box Wq = W(tQ,xo). We proceed by a duality 
argument in the spirit of Corollary 16.31 and bound \\{I + itDo)~^ Jo £sfsds\\Lp{Wo) by 
testing against h G Lq(Wo;V), l/p + 1/q = 1. As in step (iii) of the proof of [H 
Lem. 10.2], this leads to a pointwise estimate implying 



N^[iI + itD, 



Ssfsds 



< 



\S\ 



c 



y- 



Since the proof here is essentially the same as there, but replacing R" by S*", using 
area and maximal functions on S"" instead, we omit the details. The main ingredients 
are the Lp off-diagonal estimates for (/ + itD^)'^ from Lemma f6.2( i) and the tent 
space estimate [201 Thm. 1(a)] of Coifman, Meyer and Stein. 

To estimate N^,{{I +itDo)~^ Ssfsds)j^), we proceed by duality as above. We now 
instead test against h G L2{Wo]V) with /ly = and use the L2 — )■ Lg off-diagonal 
estimates for (/ + itD^)'^ from Lemma F6.2( ii) to obtain 



(/ + ttDo 



,-1 



Ssfsds 



< 



c 



y- 



It remains to see that, when n = 1, the A'^* estimate also applies to the tangential 

part W\\. Consider the transformed conjugate coefficients B = A and Bq = Ax from 

the proof of Proposition l5.4l and let E := Bq—B. Then / := J*/ solves {dt+DB)f = 

0, which yields the estimate of || A^*(u;||) II2 since (5'a/)ii = {J*SAf)± = (S^f)^. This 
completes the proof. □ 

Remark 14.3. Note that the proof also shows a priori estimates for the opera- 
tors Sa when / is not supposed to be a conormal gradient of a solution. Assume 
ll^llcnLoo < 00. If n > 2, then we have for each p < 2, 



When n 



\\N^{SAf)h + \m{SAf).)h<\\S\ 
1, we have for each 6 > 0, 



cnLc 



.11/1 



y^ 



fey. 



mSAf)h<mcnLj\f\\ys, 



feys. 
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15. Almost everywhere non-tangential convergence 

Since solutions are not defined in a pointwise sense, the classical notion of non- 
tangential convergence at a boundary point x is replaced here by 

lim |H^°(rx)|~^ / h{y)dy exists, 
Jw°{rx) 

which we call convergence of Whitney averages at x because the region W°{rx) is 
a Whitney ball. Note that since the Whitney balls at x cover a truncated cone 
with vertex x, it really amounts to a non-tangential convergence. Besides, a slight 
modification of the proofs below yields limits of averages on Whitney regions W"{z) 
for z in a fixed cone with vertex at Xq, as |z| — )■ 1. The exact choice of the Whitney 
balls does not matter. 

Definition 15.1. Let h he a. function in 0^+" with range in the bundle V in the 
sense that h{rx) G for all r > and x G S*". Let Xq G 5"" and 1 < p < oo. 
We say that the Whitney averages of h converge at xq in Lp sense to c G Vx^ if for 
any/some section c^^ G C°°{S^] V) with c^oixo) = c, 

hm \W°{rxo)\-' [ \h{y) - Cx,{y)\Pdy = 0. 

Jw°irxo) 

Here W°{x.) denotes a Whitney ball in O^"'"'^ centered at x. We say that the Whitney 
averages of h converge in Lp sense almost everywhere to Hq with respect to surface 
measure if this happens with c = hQ{xo) for almost every point xq G 5". For 
functions with values in a trivial bundle, the sections c^q are just constant functions. 

Note that the limit does not depend on the choice of the section Cxq, so this 
explains the "any/some" and it suffices to prove the existence of the limit for one 
chosen section. Clearly this notion entails convergence of Whitney averages. 

Theorem 15.2. Let A be coefficients with \\S\\cnLoo < oo. Let u be a y° -solution to 
the divergence form equation with coefficients A and let Ui be the boundary trace of 
u given by Corollaru \12.8{ . Then Whitney averages of u converge in L2 sense almost 
everywhere to Ui. In particular, Whitney averages of u converge almost everywhere 
to Ui. 

The result also holds for the R^"^" setup of with almost identical proof. 
Proof. As in the proof of Theorem 113.21 we can write 

u{^)=e^\e~'\o + Wt)^{x), 

where x = e~*x, a = Vq G L2 with ||fo||2 ^ II Vx^illyo + l/gn Uidx\ and Ui = (fo)±. 

Let p < 2 as in the third inequality of Corollary 16.31 Let Xq be a point on S*", 
and let B{xq, t) be the surface ball centered at Xq with radius t. Adapting the usual 
Lebesgue point argument for p = 1, it is seen that for almost all points xq 

lim |E(xo,t)r^ / \vo{x) - Vxoix)\^dx = 

JB{xo,t) 

for any section Vx^ G C°°(S'"; V) with Vxo{xo) = vo{xo) and one can further assume 
Dvxo ~ 0' which in particular implies that its normal component is the constant 
scalar function {vo{xo))± = ui{xo). The key point is the identity 

(48) m(x) - Mi(xo) = (e'^*e-*^(t;o - Vx,))^^) + e'^' (wt) ^{x) , 
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which follows since Dqv^q = —aNvxQ, and hence Av^^^ = av^Q and e°"*e ^^v^q = v^q- 
From Theorem 114. 11 ||A''^,(x(<T-tt'x)||2 — as r — > 0. Thus we can assume that the 
Whitney averages of Wj^ converge to in L2 sense at Xq. It remains to show, with 

lim\W{to,Xo)\~^ [ |(e"*e"*^/i^J^(x)prftrfx = 0. 

As in Ch. VII, Thm. 4], the rest of the argument consists in using the maximal 
estimates in Theorem 110.31 with some adaptation. As we we do not have pointwise 
bounds on the operators that substitute the Poisson kernel we also have to handle 
more technicalities. Let < coto < t with to, t < 1 to be chosen and Cq ^to < t < coto- 
In the L2 average, write 

(e'^*e-*^/i,Jx = ((l+«t(T)(/ + zt5o)-^/i.Jx + (e"*e-*^/i,o-(l + zta)(/ + 2t5o)-^/i.Jx. 

For the first term, we use fl26p . Fixing t and taking only the L2 average in x, this 
gives us a bound 



J>2 

This is controlled by 



J] 2-^- (\B{xo,2H)\-' [ \K,{x)\Pdx) 

.■^9 V Jb{xo,2H) J 



where M is the Hardy-Littlewood maximal operator over surface balls on S'"', M^iK) := 
M(|/i|^)^/P, and the subscript r means that we restrict the maximal operator to balls 
having radii less than r. This control is obtained by truncating the sum at 2-' ?^ r/t 
and using that t ^ to. The average in t now yields the same bound. 

For the second term, we note that (e'^^e"*^ — (1 + ita){I + itDo)~^)v.j.^^ = 0. Thus 
we may replace h^g by Vq in this term, and write it 

(e'^V(i5o)t^o)x + {e"" - (1 + iaty'){{I + itDoyh'o)^. 
with ijj{X) := e"''*'' — (1 + The first term has estimates 

~ ~ r ~ dt 

\\N,{Xt<MtDo)vo)\\l< ||V'(t/^o)^^o||2- ^0, r^O, 

by Lemma 19.31 and square function estimates. Thus we can assume that Whitney 
averages of {e'^^ip(tDQ)vo)± converge to in L2 sense at xq. By Theorem 110.31 the 
second is controlled by 

tMP{vo){xo). 

Thus it remains to show convergence to zero of 

MP{K,){xo) + {to/T)MP{h,,){xo) + tMP{vo){xo). 

Since M^(fo) G ^2(6'") as p < 2, we can further assume for xq that Mp{vo){xo) < 00. 
For such fixed xq it follows that Mp(/i^J(xo) < Mp{vo){xo) + Mp{vx,,){xo) < 00. We 
now make MP{hxo){xo) + tM^ {vo){xo) small by choosing r small. Then choose 
to < T to make (tQ/T)MP{hxg){xo) small. All the constraints on xq are met almost 
everywhere and this completes the proof. □ 
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Remark 15.3. The proof of almost everywhere convergence for averages apphes to 
V (with NP, p < 2, if n > 2). The starting point is 

replacing fl48l) and the rest of the proof is as above. The only needed modification of 
the argument is that we now use fl^^j) instead of fl2Ul) . We obtain almost everywhere 
convergence of Whitney averages of e'^^v in Lp sense to vo for p < 2. Of course, the 
term e'^* can easily be removed in the end. This factor was needed in order to have 
g^tg-aA = / on N(L)). 

Corollary 15.4. Assume that A satisfies ||i^||cnLoo < ^ '^^^ ^■^ such that all weak 
solutions u to the divergence form equation with coefficients A, for some fixed con- 
stant c > 1, satisfy the local boundedness property 

sup |m(x)| < C (\cB\~^ [ \u{y)\'^dy 
xeB V JcB 

with a constant C independent of u and of closed balls B with cB C O"'^'''". Then 
any y°-solution to the divergence form equation with coefficients A converges non- 
tangentially almost everywhere to its boundary trace. 

The local boundedness property is a classical consequence of local Holder regular- 
ity for weak solutions. For real equations (m = 1), the latter follows from [371 [23]. 
For small complex L^o perturbations of real equations, this is from [2]. For two di- 
mensional systems (n = 1), local regularity follows immediately from reverse Holder 
inequalities described in Theorem 19.81 and Sobolev embeddings. For any dimension 
and system (m > 1,?t, > 1), with continuous in 0^+" or vmo coefficients, this is 
explicitly done in [TO] . 

Proof. Applying the local boundedness property to m — Ui{xq) on Whitney balls 
yields the desired convergence for almost every xq from Theorem 115.21 □ 

We know describe new almost everywhere convergence results for ^'''-solutions. 

Theorem 15.5. Let A be coefficients with ||^^||cnLoo < Let g be an X°-solution 
with potential u to the divergence form equation with coefficients A. Then for any 
p < 2, Whitney averages of g^ = dtu, and of {Ag)\\ = (AVxm)ii, converge in Lp sense 
almost everywhere to {gi)± and {Aigi)^^ respectively, where gi is the boundary trace 
of g given by Theorem \12.4 

Furthermore, if we have pointwise ellipticity conditions on A, then the Whitney 
averages ofV^u and d^^u converge in Lp sense almost everywhere to gi and {Aigi)^ 
respectively. 

Finally, in all cases, Whitney averages of the potential u converge almost every- 
where in L2 sense to Ui. 

Recall that pointwise ellipticity holds when m = 1 (equations) or n = 1 (two 
dimensional systems). If A is continuous in 0^+", then pointwise accretivity can be 
deduced from the strict accretivity in the sense of ([2]), for any m,n. See, e.g. [26] . 
We do not know if this convergence of VxM and d^j^u holds when m > 2 and n >2 
in general. 
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Proof. We begin with the convergence for u. It is a straightforward consequence 
of the growth — folh = 0{t) for t > in Theorem 113.21 and m(x) — Ui{x) = 
{e'^'^Vt — vo)j^{x). Let us turn to the gradient. 

By Theorem 113.21 we have ft = e~*^/o + Wt for some /o G "H and w E y*. From 
the correspondence between g and / in Proposition 13.31 it follows that, modulo a 
rescahng, {g)^n+ {Ag)\\ equals Bf. Thus we need to prove convergence of Whitney 
averages of 

BJt = e~'\Bofo) + (i?oe"*^ - e"*^i?o)/o - StC^'^'fo + BtZUf 

It is clear that any y* element has Whitney averages converging almost everywhere 
to in L2 sense. This applies to the last three terms. Indeed, we have H^Hy* < 00, 
and hence H^wHy. < 00. Also ||^^te~*^/o||3;* < II^IMk^^^/ollA" < 00. Furthermore, 
using Bo{I + UDBqY^ = {I + itB^Dy^BQ, we write 

(i?oe-*^ - e-*^i?o)/o = 5o(e-*l''^"^'^^l - (/ + ^t{DB, + aN)r')fo 

+ fio((/ + ^t{DBo + aN))-^ - (/ + ztDBor')fo 

+ ((/ + itBoD)-' - (/ + it{BoD - aN))-')Bofo 

+ ((/ + zt{BoD - aN)r' - e-*l^°^-'^^l)i?o/o. 

Square function (that is, 3^*) estimates hold for the first and fourth terms, whereas 
the second and third terms have L2 norms bounded by Ct. Hence Xt<i(-Boe~*^ — 
e-'^Bo) e 3^*. 

For the term e^*^(i?o/o) we proceed as in the proof of Theorem 115.21 modified as 
in Remark 115.31 

To complete the proof, we now assume that A is pointwise elliptic. Up to rescaling, 
we have to prove convergence of Whitney averages of the conormal gradient / of u. 
To see this, write / = BQ^i^B^f) using that Bq is now invertible in Loo(5'"; £(V)), 
seen as radial coefficients on O^"^*^. Now the same argument as above replacing Bt 
by Bq shows that the Whitney averages of B^f converge in Lp sense to -Bq/o almost 
everywhere for any p < 2. We claim that the notion of convergence in Lp-sense of 
Whitney averages is stable when p < 2 under multiplication by bounded radially 
independent coefficients. Assume that h has such a convergence property and let 
M G LooiS"^; C(y)). Select smooth sections h^^ and M^^ with /^^^(xo) = h{xQ) and 
M^oixo) = M{xo). Then take Lp(W(to,Xo) average of 

M{y)h{y) - M,,{y)KM = {M{y) - M,,{y))h{y) + M.„(y)(/i(y) - K^{y)) 

with y = e~^y G W{tQ,Xo)/ For the second term, one uses the assumption on h 
and that M^^ is bounded. For the first term, use Holder inequality with exponents 
1/p = 1/r + l/g andp < r < 2. The exponent q falls on M{y)—Mscg{y) and Lebesgue 
convergence theorem applies (this is a further almost everywhere constraint on xq). 
The exponent r falls on h which has uniform control by assumption. □ 

16. FREDHOLM THEORY FOR (/ - 5'^)"^ 

We saw in Section [12] that the invertibility oi I — Sa on X (resp. 3^) allows to 
represent X" (resp. y°) solutions through Cauchy type extensions 

/ = (J-5^)"V*XVo 
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(resp. f = {I — SaJ^^Dg^^^EqVq)). Working in the space X or y, it is clear from 
Theorem 1 1 1 . 3 1 1 hat / — 5^ is invertible provided \\S\\^, is small enough. In this section, 
we use Fredholm operator theory to relax this condition and show that it suffices to 
assume this smallness only near the boundary t = 0. Our discussion in this section 
is limited to the specific but relevant case where a = 

Theorem 16.1. Assume that \\£\\* < oo, so that Sa is bounded on X and y . There 
exists e > such that if £ satisfies the small Carleson condition 

(49) lim ||xt<T^^||* < e, 

then I — Sa is invertible on X and y . 

We remark that f HUj) is equivalent to the small Carleson condition f fTUj) . The proof 
of Theorem 116.11 requires the following lemmas. 

Lemma 16.2. Assume < oo. Then I — Sa is infective on X . 

Proof. Assume that f E X satisfies / = 5*^/. Lemma [11.51 shows that / has trace 
h~ G EqU. As X C L2(R+;L2) and / is valued in "H, we have / G L2(R+]'H). 
Extend / to G L2(R;H), letting 



ft, t > 0, 

t < 0. 



To verify that satisfies dtf^ + {DB^ + aN)f^ = in R x 5" distributional 
sense, where Bl := 5^ for t > and B^ = Bq for t < 0, consider a test function 
<i) G C^{R X C(i+'^)'") and let i,{t) := 1 - r/°(|t|/e), where r/° is the function 
from Lemma 111.11 Then 



H-dt + iB^D + aN^, f')dt 

ii-dt + (B'TD + aiV)((l - ee)0), f) + ii-dt + {B'yD + aiV)(e.0), f'))dt 
+ / Ui-dt + {B'TD + aN)(P, f')dt + e'^ / (0^, f^)dt - e'' [ (0^, fl)dt 

Jr Je J-2e 

-^O + (0o,/i")-(0o,/^") = 0, 



R 



R 



with (f)o{x) := 0(0, x), using that the equation holds both in R+ and R„. Hence 
dtf^ + {DB^ + aN)f^ = in all R X Since a = extending Proposition O 
from O^"*"" to all R^"*"" (see Remark 13.61) . we see that corresponds to a function 
G L2(R^+"; C(i+")"^) solving div^{A^g^) = 0,curlx^i = in all Ri+", with A^ 
corresponding to B^. To verify that this forces g^, and therefore and /, to 
vanish, note that for any fixed i? > we can find u such that g^ = V^u, where 
i|x|<2_R I^P*^^ ^ i|x|<2i? Ifi'^P'^^ t)y Poincare's inequality and the implicit constant 
is independent of R. Take a test function rj G C^(|x| < 2R) with 77 = 1 on |x| < 
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|x|<R 
< 



with |Vx?7| ^ R ^, and use that divy^{A^g^) = in the distributional sense to get 
/ rfx < Re / {A'g\ V^u)7]d^ = - Re / (A'g', W^ri)ud^ 

J\x\<R J J 

([ \g'\'d^Y'([ \g'?dS'\(l \g'\^d^"\\g%. 

\J R<\yL\<2R J \J\:k\<2R J \J R<\x\<2R J 

Letting R oo this shows that g^ = 0, which proves the lemma. □ 
Lemma 16.3. Assume < oo and fix t > 0. Then there are lower bounds 

||/||L2(r,oo;W) < || (-^ " 'S'a) / 1| L2{t/2,oo;-H) 5 

where the implicit constant depends on t, for all f G L2(R+;'H) such that ft = 
for t < T. 

Proof By Lemma [HSl / and f^:={I- 5^)/ satisfy {dt + DBq + (jN)f^ = {dt + 
DB + aN)f. As in Proposition 13.31 combined with Proposition 12.11 this can be 
translated to 

'diVx(Ai^O) = diY^{Ag), 
curlxfif" = curlx5^, 

in O^"'""' distributional sense, where g^ = r~*^"+^)/^((5o/°)x?^ + and gr = 

r~^^+^)l'^{{Bft)^n + {ft)\\). Write 01+^^ := {|x| < e-"}, so that 0^+" C 0^+"- 
particular, the last equation implies that there is a potential u : O^^g" -> C™ such 
that 

g-g' = V^u inO^+T, 
and we may choose M so that ||ii||i2(o^+") ~ llfi'~fi'°llL2(o^+")- Fix ^ ^ C^(0"'^"^") such 
that ?7|oi+" = 1 and supp?7 C O]^^- Using the first equation and suppf? C 0\'^'^ 
gives 

Re j {Ag, g - g'')dx = / Vx(^7m))c?x = Re ^ V^{r]u))dx 

= Re / {Aig^,r]{g - g°) + (y^r])u^d^ < ||^°||i2(oi+") 11^ " ^°llL2(oi+"r 

Note that {gr)n = r-("+^)/2(/t)|| € R(V5), so that gr G lii. The accretivity dHj) of 
Aj,, for each fixed r G (0, 1), and integration for < r < e""^ imply that 

ll^llL(oi+") ~ R^ / {^9,9)dx<Re {Ag,g-g°)dx 

+ II5'IIl2(oi+")II5'°IIl2(oi+") ~ ll5'llL2(oi+")llfi'°llL2(oi+") + ll5'°lli2(oi+")' 

and hence that H^? 11^2(0?-+") ^ II^^°IIl2(o^+")- the isomorphism ffTSjl . this translates 

to ||/|U2(r,oo;W) < ||/°||L2(r/2,oo;W) aud provcs the lemma. □ 

Lemma 16.4. Assume \\S\\^ < oo. Let t] : R+ H be a Lipschitz function, that is 
\T]{t) — r]{s)\ < C\t — s\ for all t,s > 0. Then the commutator 

[r], Sa] = vSa - SaV 
is a compact operator on L2(R+, dt; L2). 
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Proof. Write Sa = Sa — ctSa as in Theorem lll.3[ Since S'a = A ^Sa, except that 
are replaced by Eq, it is enough to show compactness of [rjo, Sa]- It suffices to 
verify that 

(50) F(A) : ft ^ !\r^{t) - r/(s))Ae-(*-^)^/,rfs, 

is a compact operator on L2(R+,dt;'H). (The proof below only depends on the 
fact that A has compact resolvents.) Indeed, by duality this implies that also ft ^ 
J^{r]{t) — r]{s))Ae~^'^~^'^^ fsds is compact, upon changing A to A*. Since E^£ are 

bounded L2(R+; L2) — )■ L2(R+; V.) and commute with rj, we conclude that [77, Sa] is 
compact. 

Consider the symbol 



F(A) -.ft^ ivit) - v{s))\e-^'^''^'fsds. 



To estimate the norm of this integral operator, acting in L2(R+; C) for fixed A G 
S^^_^_, we apply Schur estimates as in [H Lem. 6.6]. We need to estimate 

sup / |(r/(t) -r/(s))Ae-(*~^)^|rfs + sup / \{7]{t) - 7]{s))Xe-^'-'^^\dt. 

t>0 Jo s>0 Js 

Using Lipschitz regularity, the first integral has estimate 

[ (t-s)Aie^(*-^)^^dt = Ar^ / \e-''dx<\~\ 
Jo Jo 

where Ai := ReA ~ |A| for A G S°^_^, and a similar estimate for the second integral 
gives the bound 

||-^(A)||l2(R+;C)->L2(R+;C) < A"^ 

It is also clear that E{X) defines a compact operator on L2(R+;C) (for example 
truncate the kernel and show from the Schur estimates that F{X) is a uniform limit 
of Hilbert-Schmidt operators). 

Consider now the Dunford integral 

F(A) = — / F(A)(A- A)~MA, uj<9<u. 

From the compactness of F(A) : I^2(R+; C) -> iv2(R+; C), and of (A - A)~^ : H 
Ti by Proposition 14. 3[ we deduce compactess of F(A)(A — A)~^ : L2(R+;'H) — )■ 
L2(R+; "H) (for example by approximating (A — A)~^ uniformly by finite rank oper- 
ators). Since ||F(A)(A — A)~^|| < A~^, the Dunford integral converges in norm, at 
least when a > 0, and we conclude that F{A) is a compact operator on L2(R+;'H) 
(for example, approximate with Riemann sums, using norm continuity of A H- 
F(A)(A — A)~^). In dimension n = 1, i.e. cr = 0, note that A = does not be- 
long to the spectrum of Do on'H. Hence it is not needed to integrate through A = 
in the Dunford integral, in which case the Dunford integral converges in norm also 
here. This proves the lemma. □ 

Lemma 16.5. Assume \\S\\^ < 00. Let < a < b < 00 and write xo '■= X{o,a) one? 
Xoo '■= X(6,oo) for the characteristic functions of these intervals. Then 

XoSaXoo : X X, and XooSaXo ■ y y 
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are compact operators. 

Proof. As in the proof of Lemma I16.4[ we may replace Sa by Sa as straightforward 
modifications of the proof below give the result for Sa- (i) We claim that the integral 
operator 



is a Hilbert-Schmidt (hence compact) operator F{\) : L2{0,a; sds) — L2 {b, oo; dt). 
Indeed, a straightforward calculation shows that 



00 pa 



\\e-^'-'^^\^sdsdt < -^e-^(^-^\ 

As in the proof of Lemma 116.41 it follows by operational calculus that 

L2(0, a; sds; H) ^2(6, oo;H): ft^ f Ae-^'-'^^fsds 

Jo 

is compact. Since EqS is bounded on ^2(0, a; srfs; "H), this proves that XoqSaXo '■ 
y ^ y is compact. 

(ii) To prove that XqSaXoo : A' — )■ A' is compact, it suffices to show that 

POO 

(51) L2(6, oo;n)^X:ft^ Xoit) / Ae^^^^'^^/.ds 

Jb 

is compact, since E^S is bounded on ^2(6, 00; Ti). To prove this, we write for t < a, 

roo roo POO 

/ Ae-^'-'^^fsds= / Ae-(^+*)^/,rfs + / (J - e-2*^)Ae-(^-*)^/,c/s 

Jb Jb Jb 



/ /--p-atAX roo 

+ (Vte"("-*)-^ J- j e"^^ A3/2e-(^-"-^)^/,ds = / + //, 



where 5 > is small enough. Cauchy-Schwarz' inequality shows that the integral 
expressions in both / and // define bounded operators L-^ih, 00;%) — )■ whereas 
e~^^ = Dq^^Dqc'^^^^^) is compact on T-L by Proposition 14.31 For /, the factor 
e~*^ : H — )■ A* is bounded by Theorem 110.11 Since y* C X, boundedness of the 
first factor in // follows from boundedness of v^e"*^""*-*^ for t E (0,a), and square 
function estimates for A since ip{X) = (1 - e~2^)/yA G ^(5°+). This completes the 
proof. □ 

Proof of Theorem \16.1\ (i) Consider first invertibility in the space M. By Theo- 
rem [HISl we have ||S'yi||A'->A' ^ W^W*^ for any perturbation of coefficients £. Thus, 
for any r > 

II'S'a/IU < C'||Xi<ri^ll*ll/IU, whenever /t = for t > r, 

with C independent of r. This follows upon writing £f = {xt<T£)f ■ Under the 
hypothesis, we can choose r > such that C||xt<ri^||* < 1/2- We obtain 

ll(/ - SA)f\\x > ll/IU - ill/IU = lll/IU, whenever /, = for t > r. 



58 



PASCAL AUSCHER AND ANDREAS ROSEN ^ 



Next consider an arbitrary f E X. Pick r/o ^ C°°(Ii^) such supp?7o C [0, r] and 
770 = 1 for t < r/2. Write r/i := 1 - Vo- Then ||(/ - 5A)(r/o/)IU > |lho/IU, and 
Lemma [16.31 shows that ||(/ — SA){vif)\\x ^ ||'7i/||a'- This gives 

il/IU < WvoflU + ll^i/IU < lia - SA){vof)h + SA){vimx 

< \\rjo{I - SA)fh + II [r/o, SaUU + ||r/i(/ - 5^)/|U + II [^1, 5^]/IU 

<||(/-5^)/|U + ||[r/o,5A]/|U. 

To show that [770, Sa] : X ^ X is compact, we write 

[?7o, Sa] = Xo[vo, Sa] + (1 - Xo)[Vo, Sa] = Xo5'a(1 - Vo) + (1 - Xo)[vo, Sa], 

where Xo '■= X(o,t/4)- Hence, compactness of the first term is granted from Lemma [16.5[ 
Next, as the X and L2 norms are the same away from the boundary. Lemma [16.41 im- 
phes that the second term is compact from X ^ X. This shows that I — Sa'-X—^X 
is a semi-Fredholm operator. 

To see that it is a Fredholm operator with index 0, note that the lower estimate 
on I — Sa above goes through with £ replaced by a£, a G [0, 1]. Apply the method 
of continuity. Since / — S'a is injective on X by Lemma 116.21 it follows that it is 
invertible. 

(ii) Consider now invertibility in the space y. That / — S'a : 3^ — 3^ is a Fredholm 
operator with index follows as in (i), provided we show that [?7o, Sa\ : y ^ y is 
compact. Here we write 

[^70, Sa] = [vo, Sa]xo + [Vo, 'S'aKI - Xo) = ivo - ^)SaXo + [Vo, Sa]{1 - Xo), 

and Lemmas 116.51 and 116.41 are applied in the same way. 

To verify bijectivity, we note that A* C 3^ is a dense continuous inclusion, where 
I — Sa '■ X ^ X is an isomorphism. This implies that I — Sa '■ y ^ y has dense 
range, hence is an isomorphism since its index is 0. □ 



17. Solvability of EVPs 

17.1. Characterization of well-posedness. For A such that / — S'y^ is invert- 
ible, we introduce boundary maps and characterize well-posedness in terms of their 
invertibility. 

Definition 17.1. For coefficients A such that \\£\\* < 00 and I — Sa '■ X X is 
invertible, define the perturbed Hardy projection 

POO 

E\h := E^h - / e-'^DSsfsds, h e L2(^"; V), 
Jo 

where f := {I - SAY^e'^^E^h. Write E^ := I - E\. Here, E^ denote the Hardy 
projections associated to the corresponding radially independent coefficients Ai. 

Proposition 17.2. The operators E^ : L2{S''; V) L2(^"; V) are bounded projec- 
tions and the range E^Ti C "H consists of all traces /o of conormal gradients f of 
X° -solutions to the divergence form equation with coefficients A in O"'^"'"". 

Proof. That E^ are bounded follows from their construction. The projection prop- 
erty (E^)'^ = E^ follows from E^Eq = 0. Next, the statement about the range 
follows from Theorem 112. 3[ □ 
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Definition 17.3. For coefficients A such that \\S\\^ < oo and I — Sa '■ y ^ y is 
invertible, define the perturbed Hardy projection 

E\h := Efh - E^ / e-'^Ssfsds, h G L^iS""; V), 
Jo 

where / := {I - SA)~^De~^^E+h. Write E^ := I - E+. Here, E^ denote the Hardy 
projections associated to the corresponding radially independent coefficients Ai. 

Proposition 17.4. The operators E^ : L2(5";V) L2(5'";V) are bounded pro- 
jections and {{E^h~^)j^ ; /i"*" G -E(fL2} consists of all traces of y° -solutions to the 
divergence form equation with coefficients A in O"'^"''". 

Proof. That E^ are bounded follows from their construction. The projection prop- 
erty (E^)'^ = E^ follows from EqEq = 0. Next, the statement about the trace 
space follows from Corollary 112. 8( ii). □ 

We remark that, unlike the case of r-independent coefficients, the complementary 
projections E^ and E^ are in general not related to solutions of a divergence form 
equation in the complementary domain R^"^" \ O^"*"". 

Proposition 17.5. For coefficients A such that I — Sa is invertible on X for (i) 
and (a), or I — Sa is invertible on y for (Hi), the following hold. 

(i) The Neumann problem (with coefficients A) is well-posed in the sense of Def- 
inition if and only if 

(52) E^n ^H^:h+ ^ {E\h+)^ 
is an isomorphism. 

(a) The regularity problem (with coefficients A) is well-posed in the sense of Def- 
inition if and only if 

(53) E+n ^%\:h+ ^ (^A^^)ll 
is an isomorphism. 

(Hi) The Dirichlet problem (with coefficients A) is well-posed in the sense of Def- 
inition if and only if 

(54) KL2iS^; V) L2(5"; C"^) : ^ {E\h+), 
is an isomorphism. 

Proof, (i) The ansatz flHUl) in Theorem 112.31 gives is a one-to-one correspondence 
between /i"*" G E^l-L and conormal gradients f = (I — S a)~^ e~^^h'^ of A'°-solutions 
to the divergence form equation. Moreover, /o = E^h^ by Proposition 117.21 Under 
this correspondence, invertibility of h'^ i— >■ {E^h'^)j^ translates to well-posedness of 
the Neumann problem. The proof of (ii) is similar. 

(iii) The ansatz (1431) from Corollary I12.8( iii) gives a one-to-one correspondence 
between G EqL2 and 3^°-solutions u to the divergence form equation. Moreover, 
{E\h^)^ = Ui by Proposition 117. 4[ Under this correspondence, invertibility of 
I— {E\h^)^ translates to well-posedness of the Dirichlet problem. □ 
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17.2. Equivalence between Dirichlet and Regularity problems. We show 
that the Dirichlet and regularity problems are the same up to taking adjoints. 

Proposition 17.6. Assume that A are coefficients such that I — Sa is invertible on 
X and I — Sa* is invertible on y. Then the regularity problem with coefficients A is 
well-posed if and only if the Dirichlet problem with coefficients A* is well-posed. 

It is not clear to us whether invertibility oi I — Sa on X implies or is implied by 
invertibility of J — Sa* on 3^. Thus we assume both. We need three lemmas, the 
first being useful reformulations of invertibility of the Dirichlet boundary map, the 
second an identity between Hardy projections and the third an abstract principle. 

Lemma 17.7. The maps 
and 

E+{L2{S''; V)/H^) ^ L2(5"; C'")/C'" : h+ ^ {E+h+)^ 
are simultaneous isomorphisms. 

Proof. This amounts to mod out H^. We recall that is preserved by A and 
E^, and annihilated by so from the definition E^h~^ = E^h^ E "H"*" for E 
H^. By Lemma Ipl (^ol^+)x = {h^)± for h+ E H^, so E^{L2{S'';V)/n^) 
-^2(6'"; C™)/C™ : H- {E^h'^)^ is a well defined map. That the two maps simulta- 
neously are isomorphisms can now be verified from {{E^h'^)j^; E "H"*"} = C™. □ 

Lemma 17.8. On ^2(5'"; V) we have the duality relation 

(55) (E^r = NE\,N. 

Proof. The proof of this duality builds on the formula 

{DaJ = -NDa'N 

on L2{S^] V) from Lemma [4.21 with Ai equal to the boundary trace of A and where 
we used the notation at the end of Definition 14. 1[ Using this observation and short 
hand notation E^ = E^^, A = I-Dyij, E^ = E^, and A = |-Da*|, it follows that we 
have 

{E^y = NE^N, A* = NAN. 

Note that when n = 1, these identities can be also checked from the extensions of 
the projections in Definition 17.41 This implies that 

/•OO /"OO 

/ iNft,£tiSAf)t)dt= / iNiSA*f)s,£sfs)ds, fEyjEX, 
Jo Jo 

which follows from Fubini's theorem and the formula defining S"^ from Lemma [11. 21 
and then letting e — ?■ using boundedness on X and y. Details are left to the reader. 
Note that Sa* is defined using the coefficients £t := A\ — A* , while £t = Ai — A. 
This duality relation between Sa and Sa* clearly extends to their resolvents. 
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For h,h & L2, using the isomorphism assumption on J — S'a and I — Sa*, we let 
f = {I - SAy^e-'^E+h and f := (J - SA^y^De-'^E+h e y and calculate 

{Nh, E\h) = {Nh, E^h) - / {Nh, E^e"'^D£sfs)ds 

Jo 

= [NE^h, h)+ {NDe-'^E+h, £^{{1 - SA)'^e-'^E+h),)ds 
Jo 

= {NE^h, h)+ {N{{I - SA*)-^De-'~^Efh)u £te-'^E+h)ds 
Jo 

= (A^^o"^. ^) + / {NE^e~''^8ju h)dt = (NE^^h, h). 
Jo 

This completes the proof. □ 

Lemma 17.9. Assume that N"^ and E"^ are two pairs of complementary projections 
in a Hubert space H, i.e. {N"^)"^ = N"^ and + = /, and similarly for E^ . 
Then the adjoint operators (N^)* and [E"^)* are also two pair of complementary pro- 
jections on H* , and the restricted projection : E'^'H — )■ N^l-i is an isomorphism 
if and only if {N~)* : {E'yH* — )■ {N~)*T-i* is an isomorphism. 

Proof. This is [3 Prop. 2.52] □ 

Proof of Proposition \17.6[ We apply the abstract result as follows. Here "H is the 
Hilbert space R(-D) C L2(S'";V) = L2 and we realize its dual "H* as L2/'H^. The 

operators A^^ are A^"*" : / 1— )■ 



and N~ : f H- 



from Definition 13.11 As 



both preserve "H, their adjoints induce operators on T-L*. We choose E~^ = E\ and 
E~ = E^. By Proposition 117. 5( ii) and reformulating (l53i) using A^"*", well-posedness 
of the regularity problem for A* is equivalent to A""*" : E\M — )■ N^l-L being an 
isomorphism. By Lemma [US] this is equivalent to (A^^)* : {E^,)*T-L* — > {N~)*T-i* 
being an isomorphism. By ( l55i) with the roles of A and A* reversed, and written 
as an identity on "H* since both terms preserve "H"*", this translates into {N~)* : 
E^K* — )■ {N~)*T-i* is an isomorphism. Using the definition of E\, {N~)* = N~ and 
n* = L2/n^, this amounts to E+{L2/H^) ^ L2(5"; C'")/C"* : h+ ^ {E\h+)^ is 
an isomorphism. Using Lemma [17.71 and Proposition I17.5( iii). this means that the 
Dirichlet problem for A is well-posed. □ 

17.3. Perturbations results. Proposition 117.61 shows that it suffices to consider 
the Neumann and regularity problems and to study invertibility of the maps 
and fl53p . Note that for r-independent coefficients A = Ai, we have E\ = Eq and 
therefore {E;^h+)^ = h+ and = h+. 

Lemma 17.10. Assume that A are coefficients such that I — Sa is invertible on X . 
Then the maps 15^) and 153]) are injective. 

Proof. Assume that G Eq'H is such that {E^h'^)± = 0. As in Theorem 112.31 
let f E X be such that /o = E\h'^ , so that we are assuming (/o)± = 0. For the 
corresponding ^'"-solution g = VxM to diVxA^f = 0, Green's formula shows that 



/ {Ag,g)d:si= / {Aigi)^uidx, 
Jo^+" Js" 
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where g e X" ^ ^2(0^+"; C(i+'^)'"), (Ai^i)x = (/o)x e 1.2(5"; C") and u G 
iJ^(0^'^'"; C"). The accretivity of A then shows that (7 = 0. Hence / = and 
/i+ = Eo'^/o = 0. 

The proof that the map t— {E\h'^)^^ is injective is similar. In this case, we use 
that Ui is constant, and /o G H so that Jg„{fo)±dx = 0. □ 

We can now derive two perturbations results. Our first result is about L^o pertur- 
bation within the class of radially independent coefficients. We need two preliminary 
lemmas. 

Lemma 17.11. Let Pt be bounded projections in a Hilbert space "H which depend 
continuously on a parameter t G {—5,5), and let S : Ti ^ }C be a bounded operator 
into a Hilbert space JC. If S : PqH K is an isomorphism, then there exists 
Q < e < 5, such that S : PtH IC is an isomorphism when \t\ < e. If each 
S : PtH JC is a semi-Fredholm operators with index it, then all indices it are 
equal. 

Proof. The first conclusion is in [TU Lem. 4.3] and the second one is proved similarly 
using in addition the continuity method. □ 

Proposition 17.12. The operators x^{DBo + o'N) G C{l-L), defined for strictly 
accretive coefficients Ai G Loo{S"'] CiV)) and a G R, depend continuously on Ai 
and a . 

Proof. This is a corollary of Theorem 17. II and ^ Prop. 2.42]. □ 

Here, note that for fixed a we called this operator E^. Only its action on 1-i mat- 
ters for well-posedness issues. In particular, this does not depend on the extension 
defined in Definition 17.41 when a = 0. 

Theorem 17.13. Assume that Ax are r -independent coefficients for which the Neu- 
mann problem is well-posed. Then there exists e > such that the Neumann problem 
is well-posed for any r -independent coefficients A[ such that \\Ai — A[\\oo < e. The 
corresponding results for the regularity and Dirichlet problems hold. 

Proof. Lemma [17.111 and Proposition 117.12] give the result for Regularity and Neu- 
mann problems as in [5] . For the Dirichlet problem, apply Proposition 117.61 □ 

The second result is perturbation from radially independent to radially dependent 
coefficients. 

Theorem 17.14. Assume that Ai are r -independent coefficients for which the Neu- 
mann problem is well-posed. Then there exists e > such that the Neumann problem 
is well-posed for any r-dependent coefficients A such that limT—^o \\Xt<T^t\\* < e- Th^ 
corresponding results for the regularity and Dirichlet problems hold. 

Proof. The condition on the coefficients implies that I — Sa is invertible on X and 
/ — Sa* invertible on y by Theorem 116.11 
We write the map ( l52i) as 
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for h"^ G EqH, where ^ ll^^lh by Theorem I12.3[ By assumption the map 

E^n ^ : h+ ^ is invertible. By [H Lem. 6.9], the norm of E^n 
Ti^ : t-j- (/ii)x is ^ ||Xt<r^^t||*- Fix r small enough so that E^Ti — ?■ : 
{h~^ + is invertible. For the last term, we then have estimates 

roc l*QQ 

\\h2h < / ||e-(^-^/2)^D||2->2||^||oo||/.||2rfs < II^IU / s-'\\fshds 




Here we used the estimate ||e^(^~^/2)A^^||^ < \\Ae-is~r/2)A(^jj^_^pj^-Q^ip^^^^gy < 

{is-T/2r' + ais-T/2)-')\\gh. 

It follows that Ei^n ^ : h+ ^ (e^(^/2)A/^^^^ 

is a compact operator since 
g-(T/2)A jg compact as a consequence of Proposition 14.31 We conclude that E^H — )■ 
'H±_ : I— is a Fredholm operator with index 0. Lemma [17. 101 shows that 
it is injective, hence an isomorphism. 

Replacing normal components by tangential parts (■)|| in the proof above 
shows the result for the regularity problem. Proposition 117.61 then gives the result 
for the Dirichlet problem. □ 



17.4. Positive results. We now give examples of radially dependent coefficients 
for which one has well-posedness. Given Theorems 117.131 and 117.141 this induces 
results for perturbed coefficients. 

Proposition 17.15. If A are r -independent coefficients, and if A is a block matrix, 
i.e. Aj^ii = = Ayx, then the Neumann, regularity and Dirichlet problems with 
coefficients A are well-posed. 

Proof. By Proposition 117. 6[ it suffices to consider the Neumann and regularity prob- 
lems. Consider the projections E^ = E^. As the maps fl3^ and fl3^ act on 
EqH C "H, it suffices to consider their action on H. throughout this proof. In this 
case, we have Eq := sgn{DBQ + aN) = Eq — Eq . Consider also the "H preserving 
projections N"^ from Definition 13. 1[ Define the anti-commutator 

C:= l{EoN + NEo). 

Since Bq is a block matrix, commutes with Bq, which shows that NEqN = 
Nsgn{DBo + aN)N = sgn{N{DBo + aN)N) = -sgn{DBo - aN), using ND = 
-DN. Hence, 

C={Eo + NEoN)N/2 = (sgn(D5o + aN) - sgn{DBo - aN))N/2 
= {{DBo)'' + aY'^\{DBo + aN) - (DBo - aN))N/2 = a{{DB^f + a^)~^'\ 

and it follows from Proposition 14.31 that C is a compact operator on "H. 
We claim that 

{2E+)N+\^^^ = I + N+{2E^)U+n= I + CU^u. 

{2E^)N-\^^^ = I - N-{2E^)U-n= I - CU-n. 
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The first identity follows from tlie computation 

(2Eo+)iV+/i+ = E^{I + N)h^ = h^ + 1(1 + Eo)Nh+ 

= h+ + l{Nh+ + 2Ch+ - NEoh+) = h+ + Ch+, for all h+ G E+n, 

and the other three identities are proved similarly. This proves that the maps 
E^Ti -7- 1-Lj_ : ^ and EqH — )■ : ^ are Fredholm operators for any 
0" G R, and for a = it follows that they are isomorphisms. By Lemma 117.111 the 
indices of these operators are zero for any a G R, and Lemma 117.101 implies that in 
fact the operators are isomorphisms for a = (n — l)/2. □ 

Proposition 17.16. If A are r-independent coefficients, and if A is Hermitean, i.e. 
A* = A, then the Neumann, regularity and Dirichlet problems with coefficients A 
are well-posed. 

Proof. By Proposition 117.61 it suffices to consider the Neumann and regularity 
problems. Let G EqT-L and define ft := e~^^h^ . By Theorem 110.11 we have 
dtft + Doft = 0, lim^^o ft = h'^ and rapid decay of as t — oo. We calculate 

POO 

{Nh+, Boh+) = - / dt{Nft, Boft)dt 
Jo 

= ^ [{NDoft, Bo ft) + {Nft, BMt))dt 
= [{{NDBo + DB*,N)ft, Bo ft) + {Bo + iV5oiV)/i))rft 

POO 

= a / ift,iBo + B;)ft)dt. 
Jo 

On the last line, we used that A* = A, or equivalently Bq = NBoN, so that 
NDBo + DB^N = 0. This gives the estimate 

POO 

I - (Boh-'),) + (5o/i+)ii)| < ^ / WftWldt. 

Jo 

From this we deduce the estimate 

\\hm < Re{h^,Boh^) < \{ht, W+)x)| + ||/|lL(R,;«) < \\hlh\\h^2M\f\\Un,,ny 

This shows that the map ( 152|) is a semi-Fredholm map, if we prove that the map 
Ti —J- L2(R+; "H) : /i H- {e~^^h)t>o is compact. To see this, note that square function 
estimates for Do give the estimate 

r midt = r mDo){A-'^'f)\\i^ < \\A-'^'f\\i 

Jo Jo 

where ipt{z) ■= \/t\z\e~^^^^ , and A~^/^ can be seen to be a compact operator on "H by 
Proposition 1131 Taking = x^{.DB' + aN) in Lemma [nUl where B\ s G [0, 1], 
denotes the straight line in L^i^S"^] '^C^')) from I to So, shows that the index of the 
map (152|) is 0. By Lemma 117.101 this map is in fact an isomorphism. 

The proof for the regularity problem is similar, using instead the estimate 

\\hm < UK, {Boh^),)\ + ll/llL(R,;«) < ll/^llhll/^+b + ll/llL(R,;K)- 

□ 
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Proposition 17.17. If A is a Holder regular C^^'^~^'^{S^] C{C^^~^^^"^)), r -independent 
coefficients, for some £ > 0, then the Neumann, regularity and Dirichlet problems 
with coefficients A are well-posed. 

For the proof, we need the following lemmas. 

Lemma 17.18. Let Bq G C^/^+^(5'"; C{V)) be the matrix associated to A. Then for 
allf^gen, \i[\D\y',Bo]f,g)\<\\f\\2\\g\\2. 

Lemma 17.19. Under the same assumptions, D{\D\^^^)n'H = D{\Do\^/'^) nU with 
equivalent graph domain norms. 

Proof of Proposition | j 7. i 7[ Consider first the Neumann and regularity problems. 
Let h+ e E^n and define ft := e-*^/i+. By Theorem [mi we have dtft + D^ft = 
and lim^^o ft = and lim^^oo ft = ^ with rapid decay. We begin with the observa- 
tion that (sgn(D)/i+,/i+) = Re(V5(-div5V5)-i/2/^t,/i+). Thus |(sgn(D)/i+, < 
II II 2 II ^^11 2- Now, we calculate for fixed T > 

(sgn(D)/i+,/i+) - (sgn(Z})/r,/r) = - r dt{sgn{D)ftJt)dt 

Jo 

\sgn{D){DBo + aN)ft, ft) + {ft, sgn{D){DBo + aiV)/0)c^t 

= 2 Re r {\D\BoftJt)dt 
Jo 

= 2Re^^ [{Bo\D\'/'ft, \D\'/' ft)dt + {[\D\'/^ Bo]ft, ft))dt, 

using that sgn{D)D = \D\ and sgn(D)A^ + Nsgn{D) = in the third equality and 
Lemma [nH in the last since ft E dI\Do\^/'^) n K C D{\D\^/'^). Accretivity of Bq 
and Lemma 117.181 lead to the estimate 

r \\\D\'/'ft\\ldt < WhthWh^h + \{sgniD) frjT)] + T ||/t II2II l/^r/'/^lhc^t, 
Jo Jo 

and by absorption, to the same estimate but with last term equal \\ft\\2dt- Due 
the rapid decay of ||/t||2 when t — 00, we conclude that 



/•OO /"OO 

/ \\\D\'/'ft\\ldt<\\hth\\hth+ / II/. 
Jo Jo 



twldt. 



Since || |-Do|"'^/^/t||2 ^ || |-D|^/^/t II2 + H/ilh from Lemma [17. 19^ we may replace D by 
Dq in the left hand side. Since square function estimates for Dq give 



this implies 



r iiiDor/v.ii^c^t = r ii(tiDoi)^/v*i^«i/.+ii^^ ^ \\h 

Jo Jo 

11^+112 <l|/^^l|2||/^;||2+ / WftWldt. 

Jo 



Well-posedness of the Neumann and regularity problems now follows as in the proof 
of Proposition II 7. 161 Proposition [TTSlthen gives the result for the Dirichlet problem. 

□ 
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Proof of Lemma \n.l^ Note that agrees on "H with the (positive) Hodge-Laplace 
operator 

divs'Vs 

Vsdivs' — curl^curlg 

where curls : L2(5'^; (Tc^")™) ^ ^2(5"; a2(Tc^")"*) is the tangential curl/exterior 
derivative on S"". Since /, (7 G "H, we have 



\D\'/',Bo]f, 



{[Ay\Bo]f,9) 



and it suffices to prove that [A^^^,i?o] is bounded on L2. Since the action of Bq 
mixes functions and vector fields, some care has to be taken. 

(i) First, by functional calculus we can replace Ashy Tq = As + X for any A G R"*", 
to be chosen large later, as — (A5 + A)^/'^ is bounded. 

(ii) Next, the commutator estimate is a local problem and by a partition of unity 
argument and rotational invariance of the assumptions, we can assume that / is 
supported in the lower hemisphere and it is enough to show that ||C[7c}^'^, -Bo]/||2 < 
II /II 2 when the smooth scalar function is 1 a neighborhood of the support of /. 
Indeed (1 — ()[T^^'^, Bolf = —[[(, tJ''^], -Bo]/; where the inner commutator is seen to 
be bounded on L2. 

(iii) Now using rescaled puUback p* to R" from the proof of Theorem 17.11 yields 

p*(To/) = Ti(p*/) with 



Ti:=- 



divRnd^ "Vr^ 




d"" 

VR..(i"divRnd2-n _ rf'^-^curl^^ci^-'^curlR 



A/ 



in L2(R";C(^+")'"), with d{y) = (jyp + l)/2 inside \y\ < 1 and extended to a 
smooth function on R", with d{y) = 2 for \y\ > 2 and 1/2 < d{y) < 2 for all y. Any 
extension would do since p*f is supported in \y\ < 1. (The proof of this equality 
builds on the fundamental differential geometric fact that the standard puUback 
operation intertwines V on S"" and R", as well as curl, and the adjoint results for 
div and curl*. Note that the rescaled pullback p* from Theorem 17.11 equals the 
standard pullback on vectors, but is times the standard pullback on scalars.) A 
further calculation shows that Ti = — divRnd^VR" + R + XI, where i? is a first order 
differential operator with smooth coefficients and divRnrf^VRn acts componentwise 
on C''^"'""^™-valued functions. Note that the coefficients of R must vanish outside 
|?/| < 2 by construction. We now choose A large enough to guarantee the accretivity 
condition Re{Tig,g) > 5||c/||^i with 5 > and all g G W^{W;C^^+''^'^). Consider 
K, 7] and g as in the proof of Theorem 17.11 and = {p*)^^f] and / = {p*)~^g. We 
claim that ||^tJ''^/ — {p*)~^7fTl^^g\\2 < \\g\\2 ~ ll/lh- For both operators Tj, we use 
the identity 

poo POO ^ 

(56) T//^ = c / s'/'Ul + s'rr'ds = c (/-(/ + s'T,)-')^. 

Jo Jo ^ 

The part with s > 1 gives rise to a bounded operator for each Tj. For the s < 1 
integral of the difference, we use the identity obtained as in Theorem 17.11 



C(/ + s'To)-'f - (p*)- V(/ + s'T^r'g 



~\pT's'hmi+s'n] 
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SO that 

iic(/ + s'nr'f - ipTWii + s'T.r'gh < shh, 

using that the commutator [77, Ti] is a first order operator. 

(iii) We are reduced to showing that [tI^^,Bo] is bounded on L2(R"; C-""^^"^™) 
with Bq := p* Bo{p*)^^ of Bq on \y\ < 1 extended to a bounded matrix function of 
class (7^/^+^ on R". We now eliminate the R part of Ti. Set T2 := — divRnd^VRn + 1 
acting componentwise in L2(R'^; C^^"'"'^-''"). The chosen extension of d insures that T2 
is accretive (in fact self-adjoint) as Ti. We claim that tI^'^ — T^^^ is bounded. We use 
again ( l56l) for each Tj . The part with s > 1 gives rise to a bounded operator for each 
Ti. For the s < 1 integral of the difference, we use || (/ + s^Ti)^^ ~ (-^ + ■s^^2)~^|| ^ s 
by the resolvent formula, because Tj have same second order term. This proves the 
claim. ^ 

(iv) Hence, it remains to estimate the commutator C = [T2^'^,Bq\. Since T2 acts 
componentwise, so does T2^^ and the commutator consists of a matrix of commu- 
tators with each component of Bq. Thus it suffices to estimate C = [T2^^, b] in 
L2(R"; C), with b scalar- valued. To see this, we use the the different representation 
for T2^^^ to obtain 

ds 



POO 

C = c [s'T,e-^'^\b] 
Jo 



3/2 ■ 



The s > 1 integral is trivially bounded, using boundedness of b and s^T2e~'^ For 
s < 1, we have ||[s^T2e~'^ '^'^ ,b]\\L^^L^ < s^^^"*"^ using pointwise decay and regularity 
for the kernel of s^T2e~^ and regularity of b. See, for example, [2] where it is proved 
that under continuity of the coefficients (here d"^), the kernel of the semigroup e~^^', 
s < 1, has Gaussian estimates (this is in fact due to Aronson for real measurable 
coefficients) and Holder regularity in each variable with any exponent in (0,1), in 
particular larger that 1/2 + e. From here, the same estimates hold for sT2e~''^^ = 
—sdse~^'^^ by analyticity of the semigroup. This takes care of the s < 1 integral. 
Further details are left to the reader. □ 



Proof of Lemma \17. 19[ Recall that Dq = DB^ + aN. As before, by a representation 
formula it is easy to prove that \DBo + aN]^^"^ — \DBq\^/'^ is bounded on L2. Hence 
we may replace -Dq by DBq. We remark that "H is invariant for both D and DBq. 

As P-hBo is an isomorphism of 7/, for / G "H, / G D{\DBq\) if and only if Py^Bof G 
D(|D|) and in this case 

\\\DBo\fh ~ WDBofh ^ \\D(PnBof)\\2 ^ || |D|(Pw5o/)||2. 

Complex interpolation for sectorial operators (see [8]) shows that for / G 7/, / G 
D{\DBo\^/'^) if and only if PnBof G D{\D\^/'^) and 

\\\DBo\'/'fh^\\\D\'^\PHBof)\\2. 
Next, for / G -H n D{\D\^/^), we have \D\^/^f e n so that 

\\\D\'/'fh^\\PnBo\D\'/'fh. 

Thus it suffices to show that for / G / G D{\D\^/^) if and only if PnBof G 

we use the regularity of Bq to yield \\\D\^/^ (PnBof) - 
PnBo\D\^''^ f\\2 < II/II2 when / G "H as a direct consequence of Lemma [17.181 and 
the fact that D and Pn commute. □ 
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Remark 17.20. Using the Tl theorem, the commutator C of the proof of Lemma 
nrm is bounded on L2 when (-A + 1)^^% e BMO (and h e L^). The converse is 
also true. This can be shown to be a regularity condition between C^/^ and C^/'^^^ . 
So well-posedness holds under this condition (expressed in local coordinates on the 
coefficients of Bq). This is probably the best conclusion we can draw from this 
method. However, we suspect that should be enough in general. 

18. Uniqueness 

The following is the class of solutions in Definition 11.81 

Definition 18.1. By a T>° -solution to the divergence form equation, with coefficients 
A, we mean a weak solution of diVxAVw = in O"*^^" with ||A^°(m)||2 < 00. 

Note that unlike the previous classes, ©"-solutions are defined through an estimate 
on u itself, not on the gradient VxM. 

Under the Carleson control on the discrepancy, we know that 3^°-solutions are 
©"-solutions. We would like to know the converse. At this stage we need assump- 
tion of well-posedness in the sense of Definition II. 2[ It goes via identification with 
variational solutions for smooth data which will be also useful later. 

Lemma 18.2. Let A he coefficients such that < 00 and I — Sa is invertible on 
y and on X , and assume that the regularity problem and the Dirichlet problem in the 
sense of Definition M.^ both are well posed. Let ip G ^2(6'"; C™) be Dirichlet datum 
such that V G -^2(5""; (TcS*")™). Then the solution u to the Dirichlet problem in 
the sense of Definition coincides with the variational solution with datum if. 

Proof. By Proposition [T731 there is a unique /i"*" e E^Ti such that {E^h'^)^^ = Vsf, 
since the regularity problem is well-posed. From Lemma 17.5^ we know that D : 
EqL2 — )■ EqI-L is surjective. Let h'^ G Eq L2 be such that Dh'^ = h~^. Consider now 
(f := {E^h'^)j_. We claim that Vs'f = ^sf- Indeed, this follows from taking the 
tangential part in the intertwining formula 

DEi = E\D, 

which is readily verified from Lemma [4.21 and definitions of E\, E\. Thus (p — if is 
constant. As in the proof of Corollary 112. 8^ by adding a normal constant in EqH^ 
to h'^, we may assume that if = ip. 

Given this h'^, the solution u to the Dirichlet problem with datum f is given by 
the normal component of 

as in Corollary 1 1 2 . 8( iii) . Next, we have 

and / is the conormal gradient to the solution to the regularity problem with datum 
Vsf- In particular f e X C L2(R+ x S""; V) by Lemma [931 

Translated to O^"*"", this shows that the solution u to the Dirichlet problem with 
datum f has V^u G ^2(0^+"; C(^+")™). This shows that M is a variational solution. 
Uniqueness of the Dirichlet problem in this class completes the proof. □ 
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Remark 18.3. Note that since X C L2(R+ x S"^; V), solutions to the regularity and 
Neumann problem always coincide with the variational solutions, by the uniqueness 
of such. In the setting of the half-space, as in [71 H] , it was shown in [TT] that this 
uniqueness result does not hold. As pointed out in [TJ Rem. 5.6], the problem occurs 
at infinity for the regularity and Neumann problems, which explains why uniqueness 
holds for the bounded ball. Although the analogue of [H] for the Dirichlet problem 
on the ball is not properly understood at the moment. Theorem 119.41 below shows 
that uniqueness of solutions essentially holds also for the Dirichlet problem on the 
unit ball. 

Proposition 18.4. Let A be radially independent coefficients and assume that the 
regularity problem and the Dirichlet problem in the sense of Definition M.^ are both 
well-posed. Then all ^"-solutions are given by u = e~"^{e~^^h'^)^ for a unique 
G EqL2. In particular, the class of V- solutions is the same as the class of 
y°-solutions, and the estimate 

\\N:{u)\\1^ [ |Vxn|2(l-|x|)dx + 

holds for all weak solutions. 

Proof Let m be a ©"-solution. For almost every p G (0, 1), V^Wp G ^2(5"; (Tc5")'") 
and Up G L2{S"'; C"^). Fix such p. As in the proof of Proposition 118.21 we can find 
h+ G E+n, h+ G E+L2 with Dh+ = h+, = VsMp and = Up on 

Using radial independence, the function Up{rx) := e'^*(e~*'^/ip (here, p is fixed 

and e~* = r G (0, 1)) thus extends to a solution of the divergence form equation with 
coefficients A, and it is a variational solution by Proposition 118.21 Since x 1— )■ u(px) 
is also a variational solution and agrees with Up on S"", we conclude by uniqueness 
that u{pr-) = e'"*(e-*^/i+)x as ^^(S"; C™)-functions for all e"* = r G (0,1], and 
almost every p G (0, 1). 

From this representation, we see that the right hand side is continuous in t, with 
range in L2, so the left hand side is continuous in r. We also have ||itpr||2 ^ II II 2 ~ 
||Mp||2 for every r G (1/2, 1] and almost every p G (0, 1). The last equivalence comes 
from the well-posedness of the Dirichlet problem, and the implicit constants are 

independent of p. As sup (1 — p)~^ J ^ ||?^s||2'^'5 ^ ||A^°(^)||2 < 00, we conclude 

l/2<p<l 

that II Zip II 2 is bounded for 1/2 < p < 1. Consider a weak limit /i^ G EqL2 of a 
subsequence h^^ with p„ — 1. Reversing the roles of p and r, for almost every r < 1, 

Up^r converges in L2(5'"; C") to u^., so that u^. = e"^{e~^^h^)^. Extending to all r, 
the representation is proved. 

In particular, this shows that the classes of 3^°-solutions and of P^-solutions of 
Lu = coincide under our assumptions. □ 

Note that the full force of ||A^°(u)||2 < 00 is not used and the condition 

sup / ||Mp||2(ip < 00 

l/2<r<l Jl-2r 



ui{x)dx 



suffices in the proof of Proposition 118. 4[ 
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Remark 18.5. If A is not r- independent, we need to know that A{p-) satisfies the 
large Carleson condition for all 1/2 < p < 1 to run the argument. This is not clear 
if we just assume this for A. However, if we assume that A is continuous on 0^+" 
and satisfies the square Dini condition of Theorem II.IH then this can be checked. 

Proof of Theorem \1.7\ We consider Ai G Loo(«S'"; £(C^^+"-''")), radially independent 
coefficients which are strictly accretive in the sense of (Ej). Assume that the Dirichlet 
problem with coefficients Ai is well-posed. By Corollary 112.91 we have VrUi = 
r~'^(e~*'^Uo)x with r = e~* and Vq given by the inverse of the well-posedness map 
(l54l) from applied to mi. The assumed uniqueness of the solution u allows us to prove 
the product rule of Vr by considering VrUi as another boundary data. The existence 
of the generator with domain contained in 1^2^(5'"'; C"*) is as in [3] in the setting of 
the upper-half space. There, the if direction was deduced using the duality principle 
between Dirichlet and regularity. An examination of the argument there reveals that 
the only if direction was implicit. We can repeat the same duality argument using 
Proposition 117.61 □ 

Proof of Theorem \1.9[ By Proposition 118.41 we know that the two classes of V°- and 
3^°-solutions are the same. Thus the assumed well-posedness for 3^°-solutions carries 
over to ©"-solutions. This completes the proof. □ 



We now specialize to the case of equations (m = 1) with real coefficients, and make 
this assumption for the coefficients A throughout this section unless mentioned oth- 
erwise. For such equations the theory of solvability for the Dirichlet problem using 
non-tangential maximal control is rather complete for real symmetric equations, but 
not so much for non symmetric equations. In [21], the extensions of the tools for 
real non symmetric equations are discussed and we refer there for details. 

We have developed a strategy using square functions rather than non-tangential 
maximal functions and our goal here is to tie this up. It is convenient to introduce 
the square function 



[Tx denoting a truncated cone with vertex x and axis the line (0, x)) and the diver- 
gence form operator L := — diVxAVx- We note that a weak solution to Lu = is in 
y° if and only if S{u) G L2{S"'), the measure being the surface measure. We have 
so far studied 3^°-solutions and well-posedness in this class, which is convenient to 
denote here by well-posedness in y°. (This was called "in the sense of Definition II. 21 ' 
in the introduction.) 

Recall that by a 1)° -solution of Lu = 0, we mean a weak solution with ||A^°(u)||2 < 
oo. As said in the introduction, we may replace N°{u) by the usual pointwise non- 
tangential maximal function. For the Dirichlet problem, we shorten well-posedness 
in the sense of Dahlberg in Definition 11.81 to well-posedness in V°. 

On regular domains such as O^"^", there is always a unique variational solution 
u G VF2^(0^"*""), which is in addition continuous in 0^+*^, to the Dirichlet problem 
with data (p G (7(5*") by results of Liftman, Stampacchia and Weinberger [33] which 
extend to real non-symmetric equations (see [29]). Thus, it is natural to ask whether 



19. New well-posedness results for real equations 
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this solution satisfies ||A'^^(m)||2 ^ C'll'/'lh with C depending on the Lipschitz charac- 
ter of S"'. By a density argument, it suffices to do this for smooth tp, say ip G C^(S'"). 
If this is the case, then the Dirichlet problem (-0)2 is said to be solvable. 

From the maximum principle and Harnack's inequalities, one can study the L- 
elliptic measure w, say at 0, which is the probability measure C{S"') 3 (p ^ u{0) 
with u the above solution. The question whether u is absolutely continuous with 
respect to surface measure is central. 

The result, somehow folklore but we have not seen it stated explicitely in the 
literature, summarizing the state of the art is the following. 

Theorem 19.1. Let L = — diVxAVx be a real elliptic operator in O"^"'"", n > 1. 
Then the following statements are equivalent 

(i) The Dirichlet problem is well-posed in T>° . 

(ii) {0)2 is solvable. 

(iii) The L-elliptic measure w is absolutely continuous with respect to surface 
measure and its Radon-Nikodym derivative k satisfies the reverse Holder B2 
condition, i.e. there is a constant C < 00 such that for all surface balls B 



Proof. The proof that (ii) is equivalent to (iii) is stated for real non-symmetric 
operators in [221 p. 241]. The proof that (i) implies (ii) is trivial. For ip G (7(5""), 
the variational solution is bounded, hence satisfies ||A^°(n)||2 < 00 since O^"'"'^ is 
bounded. By uniqueness in (i), it is the unique solution and the continuity estimate 
that follows from well-posedness shows ||A^°(ti)||2 < C||</2||2- So (ii) holds. It remains 
to see (ii) implies (i). Existence and continuity estimate are granted from (-0)2- 
Uniqueness follows the argument in [2H p. 125-126], using the equivalent assumption 
(iii) instead of (ii). The extension to non-symmetric real operators is allowed from 
the details in □ 

Theorem 19.2. Let L be an elliptic operator with real coefficients. Then all weak 

solutions to Lu = satisfy ||>S'(n)||i2((i^) ^ ||^* (^) IU2(dM) /'^'^ measure fi with 

respect to L*- elliptic measure. 

Proof. This is the result of [22] where this is proved when L = L*. Aside from 
properties of solutions that are valid for all real operators, the proof relies on the 
use of [221 eq. (7) p. 101]. Since Lu = and there is an integration by parts, it has 
to be modified with g{X, Y) and u being respectively the Green's function and the 
elliptic measure of the adjoint on the domain Q there. This is why the Aoo property 
with respect to L*-elliptic measure intervenes in the hypotheses. Further details are 
in [22]. □ 

Corollary 19.3. Let L be an elliptic operator with real coefficients. Assume that 
the Dirichlet problem is well-posed in V° for L* . Then all weak solutions to Lu = 
satisfy ||S'(m)||2 < ||A''°(m)||2. In particular, V" -solutions of Lu = Q are y° -solutions 
of Lu = under this assumption. 

Proof. By Theorem 119.11 L*-elliptic measure is Aoo with respect to surface measure. 



on S' 




and vice- versa by [T9|. So ||5'(m)||3;> 



•o < \\N°{u)\\2 follows from Lemma [19.2[ 



□ 
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Note that Corollary 119.31 and Proposition 118.41 are close but incomparable. First, 
Proposition 118.41 applies to systems of equations, whereas Corollary 119.31 applies to 
radially dependent coefficient. Secondly, the well-posedness assumptions are differ- 
ent. The next results reconciles the two approaches. 

Theorem 19.4. Let L = — diVxAVx be a real elliptic operator in O^"*"", n > 1. 
Assume further that L has coefficients with lim^-^o llxt<r^t llcnLoo sufficiently small. 
Then the following are equivalent. 

(i) The Dirichlet problem is well-posed in V" for L and L* . 

(ii) The Dirichlet problem is well-posed in y° for L and L* . 

Moreover, in this case the solutions for L (resp. L*) from a same datum are the 
same. 

Proof. It suffices to prove the conclusion for L in each the assumptions are 

invariant under taking adjoint s. 

Assume (i). Uniqueness in is immediate since the class of ©"-solutions a priori 
contains the class J^^-solutions when ||i^||cnLoo < oo. Next, for the existence, there 
is by assumption a unique ©"-solution with given boundary datum if G L2{S"'). 
Since the Dirichlet problem is well-posed in V" for L*, Corollary 119.31 shows that 
this solution is in fact a 3^°-solution. 

Conversely, assume (ii). By Theorem ll9.lt it suffices to show that (-0)2 is solvable 
for L. To this end, it suffices to consider ip G C^(S'") and the associated variational 
solution u. By Proposition 118.21 which applies because of Theorem 116.11 (/ — Sa is 
invertible on X and on y) and Proposition 117.61 u coincides with the solution in the 
sense of Definition 11.21 that is, it is a 3^°-solution. Now Theorem 114.11 provides the 
non-tangential maximal estimate that shows that (-0)2 is solvable for L. □ 

Remark 19.5. In the case of radially independent coefficients (or more generally 
for continuous, Dini square coefficients) Proposition 118.41 (or the remark that follows 
it) proves the converse also for systems. 

We can generalize results of [30] to non-symmetric perturbations of r-independent 
real symmetric operators. 

Corollary 19.6. In 0^+'", the Dirichlet problem is well-posed in V° for all real 
operators L with coefficients A such that limT-_>,o ||xt<T^^t||cnLoo small enough and 
its boundary trace Ai real symmetric. 

Proof. Let Li be the second order operator with r-independent coefficients Ai. By 
Proposition 117.16^ we know that that the Dirichlet problem for Li = L\ is well- 
posed in y°. Thus, by Theorem 117. 14[ it is well-posed in y° for L and L* . Thus, 
we conclude with Theorem 119.41 □ 

We continue with generalizations of results in |2l] , where well-posedness for Dirich- 
let was obtained for real symmetric coefficients. Well-posedness for regularity (which 
we denote here by well-posedness in X°) is new. 

Theorem 19.7. Assume that A are coefficients with limT-^o II Xt<r^^i 1 1 cnLoo small 
enough and boundary trace Ai which is real and continuous. Then the Dirichlet 
problem is well-posed in 7)° and in y° , and the regularity problem in X" is well- 
posed. In particular, this holds for real continuous coefficients in 0^+" satisfying 
the Dini square condition J^w\{t)j- < oo, where WA{t) = sup{\A{rx) — A{x)\]x G 
^",1 - r < t}. 
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Proof. Let Li be the operator with coefficients Ai. Recall that under smallness 
of liniT-^o II Xi<T^^t II cnLoo 5 suffices to prove the result for Li by Proposition 117.141 
Next, by Proposition 117.6^ the regularity problem (in X°) for Li is well-posed if and 
only if the Dirichlet problem for is well-posed in y°. On applying Theorem I19.4[ 
it suffices to prove that the Dirichlet problem with coefficients Ai is well-posed in 
V", as the same would then hold for Al by symmetry of the assumptions. To do 
this, we prove that Li-harmonic measure satisfies the property (iii) in Theorem 119. 1[ 
The argument is inspired by the one of [21], p. 139-140. 

Assume ffist we work on some boundary region of O^"*"". For r small, set Qr = 
{py G (0,1) X S""; 1 — r < p < l,y E B{xo,r)} where B{xo,r) is a surface ball 
of radius r, with real radially independent coefficients Ai being the restriction of 
some matrix defined on 0^~^" that we still denote by Ai and which is close in L^o 
to the constant matrix Ai{xo). Let g he a. non-negative function supported 
on the part of the boundary of Qr/2 in S"'. Let v be the variational solution to 
the Dirichlet problem Liv = in Qr/2 and v = g on the boundary of Qr/2 in S"^ 
and = on the part of the boundary that is contained in O^"*""". Recall that 
V G W2{Qr/2) ^C{Qr/2)- By Theorem 117. 131 because Ai is Loo close to a (constant) 
matrix for which one knows well-posedness by Proposition 117.171 one can construct 
the unique solution u in O^"*""" to the Dirichlet problem in y° with u = g on S*", 
that is Liu = with \V^u\^{1 - |x|)rfx < C||^||^. As 5- G (7^(5"), we know 

on applying Lemma [18.21 that the solution u is variational, i.e. m G 14^2^(0^^"). We 
can apply Stampacchia's minimum principle to obtain ffist that m > in O^"*"", and 
next the maximum principle in Qr/2 to conclude that v < u. From there, it remains 
to repeat the argument in [21], to obtain that v{py) < C(l — p)~'^''^||5'||2 for all 
py G Qr/ii which in turn, yields an L2 estimate on the Radon- NikodjTii derivative 
of the Li-elliptic measure. 

The localisation argument as in [21], and using the continuity of Ai to cover a layer 
of the boundary with a finite number of such small Qr/21 allows us to conclude. □ 

Corollary 19.8. With the same assumption as above andn = 1, then the Neumann 
problem with coefficients A is well-posed in X° . 

Proof. By the results in Section [5l it follows that the Neumann problem for coef- 
ficients A is well-posed in X° if and only if the regularity problem for conjugate 
coefficients A is well-posed in X° . The latter follows from the previous result since 
A satisfies the same assumption as A. □ 

Remark 19.9. As in [21], the Dini square condition in the normal direction can be 
replaced by a Dini square condition in a transverse direction to the sphere. It 
suffices to perform locally changes of variables that transform the transverse direction 
to normal ones. 
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